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1.  Reference  letter,  0AL0-TSP-C11 ,  11  May  1983,  SAB. 

2.  The  Deputy  Chief  of  Staff  for  Logistics  requested  that  the  US  Army 
Concepts  Analysis  Agency  study  US  Army  transportation  workload  forecasting 
and  develop  procedures  to  Improve  the  system. 

3.  This  report  describes  the  study  approach,  the  current  forecasting 
system,  and  several  alternative  systems  and  methods  that  should  result  in 
improved  forecasts  of  over-ocean  cargo  transportation  requirements. 
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THE  PRINCIPAL  FINDINGS  of  the  work  reported  herein  are  as  follows: 

(1)  The  transportation  workload  forecasting  system  has  produced 
inaccurate  forecasts  resulting  in  inefficient  Military  Sealift  Command 
(MSC)  industrial  fund  operations. 

(2)  Accurate  forecasting  of  cargo  transportation  requirements  can  be 
accomplished  by  forecasting  at  a  single  activity. 

(3)  Either  HQ  Military  Traffic  Management  Command  (MTMC)  or  HQ,  US 
Army  Materiel  Development  and  Readiness  Command  (DARCOM)  is  a  suitable 
location  for  a  single  point  forecasting  activity. 

(4)  The  Box-Jenkins  and  Winters  Forecasting  Models  can  provide 
accurate  forecasts  when  used  in  conjunction  with  program  information. 

(5)  Changes  to  the  allocation  of  transportation  account  codes  and 
requirements  for  forecasting  shipping  mode  are  also  required  to  improve 
forecasting  accuracy. 

THE  MAIN  ASSUMPTION  on  which  the  work  reported  herein  rests  is  that  trans¬ 
portation  workload  forecasting  requirements,  contained  in  JCS  Publication 
15,  would  not  be  changed. 

THE  PRINCIPAL  LIMITATIONS  of  this  work  which  may  affect  the  findings  are  as 
follows: 

(1)  Only  the  forecasting  of  peacetime  over-ocean  surface  cargo  trans¬ 
portation  requirements  was  evaluated. 

(2)  Historical  lift  data  was  extracted  exclusively  from  MSC  records 
and  could  not  be  validated  from  Army  sources. 


THE  SCOPE  OF  THE  STUDY  was  taken  to  include  an  analysis  of  the  Army's  long- 
range  cargo  transportation  requirements  forecasting  process  and  its  impact 
on  budgets  and  transportation  costs. 


THE  STUDY  OBJECTIVE  was  to  develop  cost  effective  systems  and  methods  for 
improving  the  forecasting  of  Army  over-ocean  surface  cargo  transportation 
requirements. 


THE  BASIC  APPROACH  followed  in  this  study  can  be  defined  as:  research  was 
conducted  into  the  nature  and  extent  of  the  forecasting  problem,  to  identify 
its  impact,  and  its  systemic  and  methodological  causes.  Several  alternative 
systems  were  evaluated  based  on  their  relative  costs  and  efficiency.  Then 
a  series  of  mathematical  techniques  was  evaluated  for  suitability  as  fore¬ 
casting  tools.  Two  of  the  techniques,  the  Box-Jenkins  and  Winters  models, 
were  used  to  forecast  the  1982  cargo  transportation  requirements  based  on 
1977  to  1981  MSC  cargo  lift  data. 


REASONS  FOR  PERFORMING  THE  STUDY  are  mainly  as  follows:  recent  forecasts 
of  Army  over-ocean  surface  cargo  transportation  requirements  have  been  in¬ 
accurate.  As  a  consequence  MSC  industrial  funds  have  incurred  significant 
losses  and  the  MSC  controlled  fleet  was  not  efficiently  utilized  for  cargo 
transport.  This  study  was  directed  to  develop  methods  to  improve  the  fore¬ 
casts. 


THE  STUDY  SPONSOR  was  the  Deputy  Chief  of  Staff  for  Logistics,  who  also 
established  the  objectives  and  monitored  the  study  activities. 


THE  STUDY  EFFORT  was  directed  by  LTC  James  N.  Keenan,  Strategy,  Concepts 
and  Plans  Directorate. 


COMMENTS  AND  QUESTIONS  may  be  directed  to  CAA,  ATTN:  Assistant  Director 
for  Strategy,  Concepts  and  Plans. 


(Tear-out  copies  of  this  synopsis  are  at  back-cover.) 
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TRANSPORTATION  WORKLOAD  FORECASTING  (TWF)  STUDY 


CHAPTER  1 
INTRODUCTION 


1-1.  INTRODUCTION.  The  Department  of  Defense  transports  approximately  7.5 
million  tons  of  cargo  annualy  via  the  Defense  Transportation  System  (DTS). 
In  excess  of  50  percent  of  this  cargo  is  generated  by  the  Army.  Planning 
for  and  use  of  military  and  commercial  shipping  is  dependent  on  the 
accurate  forecasting  of  the  services'  requirement  for  movement  of  cargo. 

The  current  Army  transportation  workload  system  does  not  produce  accurate 
forecasts  of  Army  cargo  movement  requirements.  As  a  consequence,  the 
Deputy  Chief  of  Staff  for  Logistics  (DCSLOG),  Headquarters,  Department  of 
the  Army  (HQDA),  tasked  the  US  Army  Concepts  Analysis  Agency  (CAA)  to  study 
the  transportation  workload  forcasting  issue  and  to  report  the  study 
findings  in  January  1984.  This  report  discusses  the  study  approach,  the 
forecasting  system,  other  studies  of  the  issue,  other  service  solutions  to 
the  problem,  and  approaches  to  the  Army's  forecasting  problem. 

1-2.  BACKGROUND.  The  Army  is  required  to  submit  periodic  forecasts  of  its 
over-ocean  transportation  to  Military  Sealift  Command  (MSC).  Recent  Army 
forecasts  have  shown  considerable  variance  from  actual  cargo  lift.  This 
discrepancy  between  forecasted  and  actual  lift  requirements  impacts 
adversely  on  the  operation  of  the  MSC  and  Military  Traffic  Management 
Command  (MTMC)  industrial  funds,  on  the  MSC  controlled  fleet  sizing 
process,  and  on  Army  transportation  budget  preparation  and  execution. 

1-3.  PURPOSE  AND  SCOPE.  The  purpose  of  this  study  was  to  develop 
procedures  to  improve  the  US  Army  transportation  workload  forecasting 
system.  The  study  focus  was  primarily  on  the  long-range,  over-ocean 
surface  transportation  requirement  forecasting  process. 

1-4.  OBJECTIVES.  Objectives  of  the  study  were  to  determine  the  nature  and 
extent  of  the  transportation  workload  forecasting  problem  and  to  explore 
and  evaluate  alternative  solutions. 

1-5.  LIMITATION.  Historical  lift  data,  used  in  model  development,  was 
extracted  exclusively  from  Military  Sealift  Command  records  and  could  not 
be  validated  from  Army  sources. 

1-6.  ELEMENTS  OF  ANALYSIS.  The  elements  of  analysis  of  the  study  were 
designed  to  explore  those  facets  of  the  transportation  workload  forecasting 
environment  which  wduld  help  identify  the  extent  of  the  problem, 
contributing  factors,  and  potential  solutions.  These  elements  were  as 
follows: 


a.  What  are  the  recorded  variances  between  long-range  transportation 
workload  forecasts  and  actual  utilization  of  cargo  shipping? 

b.  Do  systemic  conditions  exist  which  contribute  to  unrealistic 
forecasting?  If  so,  what  are  they? 

c.  What  is  the  economic  and  operational  impact  of  long-range  forecasts 
which  are  at  variance  with  actual  utilization? 

d.  Do  short-range  forecasts  impact  on  long-range  forecasting?  If  so, 
how? 


e.  What  is  the  impact  of  the  current  separation  of  responsibility  for 
long-  and  short-range  forecasting? 

f.  What  methodologies  exist  in  the  other  services  which  could  be 
applied  to  the  resolution  of  the  Army  problem? 

g.  What  DARCOM  activities  affect  major  items  of  equipment  planned  for 
overseas  distribution? 

h.  What  is  the  impact  of  the  Total  Army  Distribution  Program  (TAEDP)  on 
forecasting  of  requirements? 

i.  What  is  the  commodity  impact  of  Army  and  Air  Force  Exchange  Service 
(AAFES)  forecasting  procedures? 

j.  What  are  the  feasible  and  cost  effective  methods  for  improving 
forecasting  accuracy? 

k.  What,  if  any,  unique  coronodities  are  marked  as  general  or  special 
cargo? 

1-7.  STUDY  TASKS.  To  fulfill  the  study  objectives  and  to  answer  the 
elements  of  analysis,  five  principal  tasks  were  identified  and  were  the 
basis  of  the  study  methodology.  These  tasks  were: 

a.  Identification  of  the  nature  and  extent  of  the  variances  between 
forecasts  and  actual  lift. 

b.  Determination  of  the  Impact  of  erroneous  forecasts. 

c.  Development  and  examination  of  forecasting  systems  and  methods. 

d.  Evaluation  of  alternative  locations  for  forecasting  responsibility. 

e.  Documenting  recommended  methodology. 


1-8.  STUDY  METH000L06Y.  The  methodology  used  in  the  study  is  depicted  in 
Figure  1-1. 


Figure  1-1.  Study  Methodology 


a.  Activities  in  the  determination  of  the  nature  and  extent  of  the 
problem  consisted  of: 

(1)  Research  into  completed  and  ongoing  studies  of  the  transportation 
workload  forecasting  problem. 

(2)  Research  into  Joint  Chiefs  of  Staff  (JCS)  publications,  MSC, 
Military  Airlift  Command  (MAC),  and  MTMC  directives,  and  Army,  Navy,  and 
Air  Force  regulations  which  were  related  to  transportation  workload 
forecasting  and  operations. 

(3)  Analysis  and  evaluation  of  the  procedures  used  by  Army 
transportation  workload  forecasting  and  budgeting  agencies. 

(4)  Analysis  of  MSC,  MTMC,  and  MAC  operations  and  the  relationship 
between  these  commands  and  the  Army  transportation  forecasting  system. 

(5)  Collection  and  analysis  of  4  years  of  forecasts  and  6  years  of 
cargo  lift  data. 


b.  Activities  in  the  determination  of  the  impact  of  long-range 
transportation  workload  forecasts  on  the  Army  budgets  and  on  the  industrial 
funds  of  MSC  and  MTMC  included: 

(1)  Examination  of  the  interfaces  between  the  forecasting  system  and 
budgetary  process  and  quantification  of  the  effects  of  incorrect  forecasts. 

(2)  Review  of  HQDA  budget  documents  and  interviews  with  personnel 
involved  in  the  budget  process. 

(3)  Examination  of  the  development  and  operation  of  MSC  industrial 
funds. 

(4)  Investigation  of  the  MSC  controlled  fleet  sizing  process. 

c.  Actions  in  the  development  and  examination  of  alternative  methods 
for  increasing  forecasting  accuracy  consisted  of: 

(1)  Analysis  and  evaluation  of  previously  collected  performance  data, 
focusing  on  systemic  contributions  to  erroneous  forecasts. 

(2)  Identification  of  the  forecasting  systems  data  requirements  and 
the  most  accurate  data  sources. 

(3)  Formulation  of  alternative  forecasting  systems  and  evaluation  of 
costs  in  terms  of  personnel  and  facilities. 

(4)  Investigation  of  Navy  and  Air  Force  forecasting  systems. 

(5)  Evaluation  of  mathematical  forecasting  techniques  to  determine 
their  suitability  as  tools  in  transportation  workload  forecasting. 

(6)  Application  of  several  mathematical  techniques  to  historical 
cargo  lift  data  to  determine  the  most  appropriate  technique  and  to  identify 
the  parameters  of  the  forecasting  model. 

d.  Examination  and  evaluation  of  alternative  locations  for  forecasting 
responsibility  consisted  of  the  following  actions: 

(1)  Evaluation  of  the  availability  and  accuracy  of  the  information 
required  to  manage  the  transportation  workload  forecasting  system. 

(2)  Identification  and  evaluation  of  the  availability  of  data 
processing  support  which  is  required  by  proposed  systems. 

e.  Documenting  the  recommended  methodology  consisted  of  documentation 
of  the  Box-Jenkins  and  Winters  models  and  the  actions  required  to  improve 
the  systemic  and  external  aspects  of  the  forecasting  system. 


CHAPTER  2 

THE  CURRENT  ARMY  TRANSPORTATION  WORKLOAD  SYSTEM 


2-1.  GENERAL.  The  Army's  transportation  workload  forecasting  system  is 
governed  by  Joint  Chiefs  of  Staff  (JCS)  Publication  15,  AR  55-23,  AR  55-30, 
and  several  MSC,  MAC,  and  MTMC  directives.  This  chapter  defines  the  regu¬ 
latory  basis  for  the  system,  observations  on  how  the  actual  system  works, 
the  performance  of  the  system  in  terms  of  forecasting  accuracy,  the  effects 
of  inaccurate  forecasts,  and  some  external  factors  which  may  contribute  to 
the  inaccuracies. 

2-2.  THE  REGULATORY  FORECASTING  SYSTEM 

a.  JCS  Pub  15,  Mobility  Systems  Policies,  Procedures  and  Considerations, 
2  June  1974,  contains  approved  joint  transportation  procedures  applicable 
to  the  submission  of  common  user  movement  requirements.  Specifically  Chap¬ 
ter  4  (Transportation  Requirements,  Allocations,  and  Priorities)  addresses, 
inter  alia,  shipper  service  forecasted  cargo  movement  requirements.  It 
requires  that  each  military  service  and  the  Defense  Logistics  Agency  (DLA) 
submit  four  specific  forecasts  of  sealift  requirements. 

(1)  On  1  May  a  preliminary  annual  forecast  (MSC-9)  is  submitted  which 
states  the  worldwide  MSC  surface  movement  requirements  for  the  fiscal  year 
which  begins  17  months  later  (e.g.,  1  May  84  for  fiscal  year  1986). 

(2)  An  annual  forecast  (MSC-10)  is  submitted  on  1  March  for  the  sub¬ 
sequent  fiscal  year  (e.g.,  1  Mar  85  for  FY  86).  This  forecast  refines  the 
preliminary  forecasts. 

(3)  A  sealift  cargo  requirement  (short  range)  is  submitted  by  the 
fifteenth  day  of  each  month  for  the  succeeding  3  months.  Each  of  the  re¬ 
ports  states  the  monthly  sealift  cargo  requirements  in  measurement  tons  for 
each  traffic  route,  program,  commodity,  and  type  of  shipment  or  mode. 

(4)  Change  reports  are  required  when  significant  changes  to  the  above 
forecasts  are  anticipated. 

b.  Each  military  service  is  also  responsible  for  the  collection  and 
submission  of  movement  requirements  of  government  agencies  outside  DOD  for 
which  the  service  has  sponsorship  responsibility  and  whose  requirements 
have  been  approved  by  competent  authority  as  eligible  to  be  handled  within 
the  DOD  transportation  system. 

c.  AR  55-23,  Military  Sealift,  implements  JCS  Pub  15  within  the  Army. 

It  identifies  57  numbered  traffic  areas  and  their  associated  geographic 
areas.  These  areas  are  the  terminals  of  the  traffic  channels  for  which 
forecasts  are  submitted.  Additionally,  AR  55-23  identifies  sponsor  codes, 
budget  programs,  cargo  classes/commodities,  types  of  shipment,  and  formats 


for  reports  submitted  to  the  Military  Sealift  Command.  These  are  discussed 
below. 

(1)  Budget  programs  to  be  used  in  the  forecasts  are  Troop  Support, 
Military  Construction,  Military  Aid,  and  Civilian  Aid. 

(2)  Cargo  classes/commodities  to  be  forecasted  in  the  reports  are: 

Reefer-chill  (CHL) 

Reefer-freeze  (FRZ) 

Coal  and  Coke  (COK) 

Bulk -Other  (BLK) 

Privately  Owned  Vehicles  (POV) 

Household  Goods  (HHG) 

Ammunition  and  Explosives  (AMO) 

General  Cargo  (GEN) 

Special  Cargo  (SPC) 

Assembled  Aircraft  (AAC) 

Empty  Conex  (CNX) 

Cargo  Carrying  Trailers  (CCT) 

(3)  Type  shipments  or  modes  to  be  forecasted  in  the  reports  are  break- 
bulk,  container,  and  MILVAN. 

d.  AR  55-30,  Space  Requirements  and  Performance  Reports  for  Transporta¬ 
tion  Movements,  prescribes  procedures  for  the  preparation  and  submission  of 
cargo  requirements  and  performance  reports  and  defines  responsibilities  for 
report  submission. 

(1)  Responsibilities  defined  in  AR  55-30  are  as  follows: 

(a)  The  Deputy  Chief  of  Staff  for  Logistics,  HQDA,  is  responsible 
for  developing  longe-range  cargo  movement  requirements  (preliminary  and 
annual  forecast  reports)  and  for  programing  and  budgeting  for  transporta¬ 
tion  services. 

(b)  DARCOM  Logistics  Control  Activity  (LCA)  has  DA  responsibility 
for  developing  and  programing  short-range  movement  requirements. 

(c)  The  following  commanders  and  agency  heads  are  responsible  for 
the  provision  of  inputs  to  both  long-range  and  short-range  forecasts: 

Military  Postal  Service  Agency 

Army  and  Air  Force  Exchange  Service 

Armed  Forces  Courier  Service 

US  Army  Intelligence  and  Security  Command 

Chief  of  Engineers 

National  Security  Agency 

Ballistic  Missile  Defense  Systems  Command 

US  Army  Communications  Command 


US  Army,  DARCOM,  Logistics  Control  Activity 

US  Army,  Europe 

US  Army,  Japan 

Eighth  US  Army 

Western  Command 

US  Army  Forces  Command 

193d  Inf  Bde  (Panama) 

172d  Inf  Bde  (Alaska) 

Deputy  Chief  of  Staff  Logistics,  HQDA 

(2)  The  commands  and  agencies  listed  in  paragraph  2-2d(l)(c)  above, 
are  required  to  submit  their  long-range  reports  to  US  Army  Management  Sys¬ 
tems  Support  Agency  (USAMSSA).  USAMSSA  provides  a  consolidated  report  to 
the  Director  for  Transportation,  Energy  and  Troop  Support,  ODCSCLOG,  who 
analyzes  and  adjusts  the  stated  requirements.  The  adjusted  data  is  then 
provided  to  USAMSSA  for  preparation  and  submission  to  MSC  and  MTMC.  As  an 
exception,  the  US  Army  Installation  Support  Activity,  Europe,  Enerqy  Center, 
Rheinau  is  required  to  submit  long-range  solid  fuel  (coal  and  cokej  require¬ 
ments  semiannually  to  HQDA  ODCSLOG. 

(3)  Short-range  requirements  for  surface  cargo  movement  are  to  be 
submitted  monthly  to  DARCOM  LCA.  LCA  is  required  to  consolidate  the  reports 
and  forward  the  Army's  statement  of  requirements  to  MSC  and  MTMC. 

(4)  Change  reports  are  to  be  submitted  when  there  is  a  600-measurement 
ton-change  over  a  traffic  area  (e.g.,  Gulf  Coast  to  Europe). 

e.  Military  Sealift  Command  uses  the  long-range  forecasts  to  plan  the 
use  of  its  nucleus  fleet  and,  when  required,  plans  for  augmentation  by  com¬ 
mercial  or  National  Defense  Reserve  Fleet  resources.  (The  total  of  these 
assets  is  the  MSC  controlled  fleet.)  The  long-range  forecasts  are  also 
used  by  MSC  to  formulate  the  shipping  rates  to  be  charged  to  the  services 
for  cargo  shipped  and  by  the  services  for  budget  preparation.  The  stated 
use  of  short-range  forecasts  by  MSC  and  MTMC  is  the  scheduling  of  ship  and 
port  workloads. 

f.  JCS  Pub  15  directs  that  monthly  utilization  reports  be  provided  to 
the  Army  comparing  final  forecasted  requirements  with  actual  cargo  for  a 
particular  month. 

g.  The  Army's  forecasting  system  is  portrayed  at  Figure  2-1,  and  the 
report  submission  sequence  is  shown  at  Figure  2-2. 
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2-2.  OBSERVATIONS  OF  THE  ARMY  FORECASTING  SYSTEM 

a.  Long-range  Forecasts.  The  long-range  forecasts,  as  required,  are 
submitted  by  the  forecasting  agencies  identified-  in  AR  55-30.  However, 
methods  and  procedures  for  preparation  of  these  reports  are  inconsistent 
and  vary  extensively  in  their  accuracy.  A  synopsis  of  forecasting  proce¬ 
dures  by  the  various  commands  follows. 

(1)  DARCOM  forecasts  for  major  items  and  ammunition  are  prepared  at 
the  DARCOM  major  subordinate  command  inventory  control  points  (ICP).  The 
reports  are,  in  theory,  based  on  The  Army  Equipment  Distribution  Plan 
(TAEDP),  Materiel  Movement  Reports  (MMR),  Logistics  Intelligence  File 
(LIF),  and  forecaster  judgment.  Due  in  part  to  forecaster  mistrust  of  the 
TAEDP  accuracy,  inability  to  accurately  correlate  LIF  data  with  actual 
movement,  and  the  absence  of  accurate  timely  feedback,  significant  sub¬ 
jective  forecasting  is  made.  Additionally,  indications  are  that  there  is 
inadequate  information  exchange  between  the  item  managers  and  the  personnel 
making  the  forecasts.  Consequently  the  program  information  available  to 
the  transportation  workload  forecaster  may  not  be  the  most  current. 

(2)  DARCOM  forecasts  for  secondary  items  are  prepared  by  the  Logistic 
Control  Activity  using  the  previous  year's  actual  lift  data  extracted  from 
MSC  billing  records. 

(3)  Army /Air  Force  Exchange  System  forecasts  are  prepared  using  the 
arithmetic  mean  of  three  previous  years  of  lift  history.  This  is  adjusted 
by  anticipated  sales  growth;  inflation  factors,  changes  in  troop  strengths; 
opening,  closing,  or  consolidation  of  facilities;  and  changes  in  mode  of 
transportation. 

(4)  The  Military  Postal  Service  Agency  (MPSA)  uses  historical  data  on 
the  volume  of  consolidated  Army  and  Air  Force  mail  moved  over  an  MSC  channel. 
Adjustments  are  made  based  on  knowledge  of  strength  increases,  new  missions, 
and  changes  to  traffic  routes.  History  is  developed  from  MSC  billing  tapes 
which  are  provided  directly  by  MSC  to  the  MPSA.  An  automated  management 
information  system,  the  Military  Automated  Mail  Accounting  System  (MAMAS), 

is  used  by  MPSA  in  their  forecast  formulation. 

(5)  The  Chief  of  Engineers  forecasts  are  submitted  only  by  the  US 
Army  Engineer  District,  New  York,  and  apply  only  to  cargo  originating  at 
Bayonne,  NJ  and  destined  for  Thule,  Bremerhafen,  Rotterdam  and  the  Azores. 

The  forecasts  are  based  on  previous  forecasts  and  have  not  changed  for 
several  years.  Forecasts  are  not  made  by,  or  for,  other  engineer  commands 
or  agencies. 
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(6)  The  Ballistic  Missile  Defense  Systems  Command  forecasts  are  pre¬ 
pared  by  the  Logistics  Support  Contractor  (LSC)  at  Kwajalein  Missile  Range 
from  inputs  submitted  by  six  range  contractors  who  perform  range  services. 
Forecasts  are  based  on  3  years  of  historical  data  as  determined  by  the  con¬ 
tractor  records  of  items  shipped.  The  quantities  obtained  from  the  contrac¬ 
tors'  historical  data  have  not  agreed  with  the  cargo  movement  feedback 
reports  provided  by  LCA. 

(7)  The  US  Army  Communications  Command  reports  that  their  forecasts 
are  "strictly  guesswork." 

(8)  US  Army  Europe  has  not  submitted  the  required  surface  cargo  move¬ 
ment  forecasts  for  at  least  2  years. 

(9)  Eighth  US  Army  forecasts  are  derived  from  data  provided  by  nine 
subordinate  elements  based  on  forecast  history. 

(10)  US  Army  Japan  (USARJ)  forecasts  are  based  on  history  and  projected 
requirements.  USARJ  does  receive  monthly  feedback  reports  from  LCA. 

(11)  Western  Command  forecasts  are  straight  line  projections  from  his¬ 
torical  data  adjusted  by  experience. 

(12)  193d  Infantry  8rigade  (Alaska)  uses  monthly  tonnage  reports  com¬ 
piled  locally  as  the  historic  basis  for  their  forecasts. 

(13)  172d  Infantry  Brigade  (Panama)  forecasts  are  based  on  historical 
cargo  lift  data  provided  by  the  local  MTMC  facility. 

(14)  ODCSLOG,  HQDA  consolidates  and  adjusts  the  aggregated  inputs  of 
the  Army  forecasting  commands  and  agencies.  It  also  prepares  the  long- 
range  CONUS  outbound  household  goods  and  POV  cargo  space  requirements  fore¬ 
casts.  The  household  goods  and  POV  forecasts  are  based  on  the  most  recent 
year  of  complete  cargo  movement  data  available. 

b.  Short-range  Forecasts.  Short-range  forecasts  are  prepared  by  all 
activities  which  submit  long-range  forecasts  and  by  Headquarters,  Forces 
Command . 

(1)  Many  of  the  short-range  forecasts  are  a  straight  line  monthly 
average  of  the  annual  forecast  and  do  not  account  for  monthly  or  seasonal 
variations.  Others,  particularly  DARCOM,  use  program  and  historical  data. 

(2)  Forces  Command  does  not  prepare  a  long-range  forecast.  It  does 
prepare  and  submit  a  short-range  forecast  for  all  CONUS  outbound  household 
goods  and  privately  owned  vehicles,  although  it  does  not  have  knowledge  of 
all  movements.  In  FY  81  and  FY  82,  household  goods  were  approximately  20 
percent  and  11  percent,  respectively,  overforecast  in  the  long-range 
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forecasts  and  approximately  80  percent  and  38  percent,  respectively,  over 
forecast  in  the  short-range  forecasts. 


(3)  LCA  consolidates  the  short-range  forecasts  provided  by  the  fore¬ 
casting  agencies  and  forwards  the  report  to  MSC  and  MTMC.  Additionally, 

LCA  prepares  the  short-range  forecast  of  DARCOM  secondary  items  using  a 
modified  Kalman  filter  algorithm  and  data  from  the  LIF.  A  recent  DARCOM 
Inventory  Research  Office  study  concluded  that  the  current  smoothing  al¬ 
gorithm  was  not  optimal  and  recommended  a  replacement. 

c.  Performance  Reports 

(1)  Chapter  3,  AR  55-30,  requires  DARCOM  LCA  to  produce  monthly  feed¬ 
back  reports  for  each  reporting  command  or  agency  which  reflect  the  fore¬ 
casted  and  actual  lift  of  surface  cargo  tonnage.  The  data  provided  are  to 
be  used  as  guidance  in  the  preparation  of  the  long-range  and  short-range 
forecasts.  Most  Army  forecasters  cite  there  is  a  lack  of  information  on 
shipments  and  identify  this  as  a  barrier  to  improved  forecasts. 

(2)  Feedback  reports  produced  by  LCA  are  extracted  from  the  MSC 
billing  tapes  and  the  forecasts  provided  by  the  forecasting  agencies  or 
commands.  LCA  attempts  to  identify  the  cargo  shipper  either  through  a  com¬ 
parison  of  the  shipping  information  to  the  LIF  or  through  the  Transporta¬ 
tion  Account  Code  (TAC).  The  LIF  to  shipping  information  comparison  cannot 
identify  a  large  number  of  transactions  due  primarily  to  the  absence  of  the 
information  within  the  LIF.  Identification  of  the  shipper  through  TACs, 
under  which  all  shipments  are  reported,  is  not  currently  achievable.  This 
is  primarily  a  result  of  the  current  TAC  structure  which  allocates  a  single 
TAC  to  all  cargo  shipped  by  DARCOM  and  several  other  commands  and  agencies. 

d.  Cargo  Transportation  Budgeting.  The  MECHTRAM  (Mechanization  of 
Selected  Transportation  Movement  Reports)  system  produces  long-range  cargo 
forecast  reports  and  provides  them  to  DCSLOG  and  DCSPER  to  be  used  in  bud¬ 
get  preparation.  Review  of  the  budget  formulation  process  indicates  that 
the  MECHTRAM  reports  are  not  used  and  that  other  tools  have  been  developed. 
MILPERCEN  personnel  reports,  MSC/MTMC  provided  shipping  rates  and  histori¬ 
cal  factors  are  the  basis  of  the  transportation  costs  within  Military  Per¬ 
sonnel  Army  (MPA)  budget,  while  history  and  program  knowledge  are  the  basis 
of  the  second  destination  transportation  costs  in  Military  Assistance  Pro¬ 
gram  (MAP),  and  Operations  and  Maintenance  (OMA)  budgets.  None  of  the 
forecasting  commands  or  agencies  have  any  responsibility  for  budget  formu¬ 
lation  and  execution.  Payment  of  shipping  charges  is  made  by  US  Army 
Finance  and  Accounting  Center  based  on  MSC  and  MTMC  billings. 
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2-3.  FORECASTING  SYSTEM  PERFORMANCE 

a.  Forecasts  submitted  by  the  Army  in  recent  years  have  generally  over¬ 
stated  requirements.  While  the  total  forecast  compared  to  actual  shipment 
has  shown  significant  improvement,  forecasts  by  commodity  and  type  of  ship¬ 
ment  have  continued  to  exhibit  significant  variances.  Table  2-1  shows  in 
thousands  of  measurement  tons  (MTON)  the  aggregated  long-range  forecast  and 
the  actual  lift  and  error,  by  commodity  and  type  shipment  for  FY  1981. 

Table  2-2  portrays  the  same  data  for  FY  1982.  It  should  be  noted  that 
overforecasts  for  some  commodities  in  some  modes  are  offset  by  underfore¬ 
casts  in  the  same  commodities  in  other  modes.  For  example,  breakbulk 
general  cargo  was  overforecast  by  115,000  tons  and  seavan  general  cargo  was 
underforecast  by  191,000  tons. 

b.  Forecasting  errors  by  the  Army  in  FY  81  were  much  greater  than  those 
of  other  services.  In  FY  82,  the  differences  were  less  pronounced  as  Army 
forecasts  were  more  accurate.  As  the  Army's  cargo  accounts  for  approximately 
60  percent  of  all  cargo  carried  by  MSC,  the  effects  of  percentage  errors 

are  much  more  severe  than  errors  by  the  other  services.  Tables  2-3  and  2-4 
show  comparisons  of  the  aggregate  1981  and  FY  82  breakbulk  forecasts  and 
shipments  by  the  shipper  services. 


Table  2-1.  FY  81  Forecast  Performance  (000  MTON) 


Col  4ltlos 

Breakbulk 

Sear a n 

MILVAN 

Forecast 

1  L,,t  1 

Error 

Forecast 

r-n 

Error 

Forecast 

1  l'ft  J 

Error 

General 

MS 

170 

-216 

1.736 

1,610 

-126 

11 

21 

-10 

Household 

as 

42 

-3 

65 

36 

-29 

-- 

13 

♦13 

Special 

671 

575 

-92 

19 

20 

♦1 

- 

-- 

-- 

POT* 

350 

786 

-64 

136 

190 

♦62 

1 

-- 

-1 

Reefer 

0 

1 

♦l 

36 

.  41 

♦5 

- 

-- 

- 

^munition 

16$ 

120 

-45 

-- 

3 

♦3 

20 

20 

-8 

CCT 

54 

76 

♦22 

-- 

-- 

- 

- 

" 

- 

cm 

S 

2 

-3 

- 

-- 

- 

- 

” 

- 

Coal/balk 

462 

377 

-as 

-- 

2 

♦2 

-- 

-- 

-- 
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Table  2-4.  FY  82  Comparison  of  Service  Forecasts 


FY  1982  breakbulk  cargo  (000  MTON) 


FY  1982  seavan  cargo  (000  HTOH) 


Component 

[  Forecast 

|  Actual 

|  Error 

|  Percent 

|  Forecast 

|  Actual 

j  Error  | 

Percent 

Army 

1.860 

1,407 

453 

24  over 

1,874 

2,186 

312 

17  under 

Navy 

375 

393 

18 

18  over 

918 

934 

16 

2  under 

Air  Force 

523 

463 

16 

11  over 

557 

684 

127 

23  over 

2-4.  EFFECTS  OF  FORECASTING 

a*  General.  Erroneous  forecasting  has  several  negative  effects.  It 
results  in  inefficient  MSC  shipping  and  MTMC  port  operations,  as  measured 
by  losses  in  MSC  and  MTMC  industrial  funds,  in  ship  utilization,  in  unre¬ 
liable  service  budgets,  and  in  inaccurate  MSC-controlled  fleet  sizing. 

b.  Industrial  Fund  Impact.  The  MSC  industrial  fund  losses  for  FY  82 
attributable  to  erroneous  Army  forecasts  are  shown  in  Table  2-5.  Losses  to 
the  industrial  fund  are  initially  absorbed  by  the  MSC  but  are  recouped  in 
the  following  year  through  increases  in  the  rates  charged  to  the  shipper 
services.  This  would  appear  to  be  an  accounting  problem.  However,  the 
industrial  fund  losses  do  reflect  inefficient  use  of  cargo  carrying  re¬ 
sources  and  expenditures  in  the  subsequent  year  which  would  not  occur  if 
the  forecasts  were  accurate. 


•  V  y->>v-V'X'!v'vA'!vAiv' ' 
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Table  2-5.  Losses  in  MSC  Industrial  Fund 


FY  1982  rates  (developed  on 
preliminary  requirements) 

Actual  FY  1982  Army  lift 

Difference 


Tons 

Miles 

(000) 

(000.00) 

2,076 

8,378  : 

1,470 

5,240 

(  606) 

(3,138)  : 

c.  Budget  Impact.  Budget  preparation  for  second  destination  (over, 
ocean)  costs  of  shipping  cargo  is  impeded  by  the  lack  of  accurate  forecasts 
of  the  cargo  movement  requirements.  As  stated  earlier,  budget  preparers 
have  developed  their  own  tools  to  produce  the  budget.  OMA  (p  7)  is  the 
predominant  source  of  funds  for  over-ocean  cargo  movement.  Table  2-6  shows 
the  FY  82  performance  of  this  budget  element. 


Table  2-6.  Army  Budget  Performance 


FY  82  budget  OMA  (p  7) 


(1)  The  MSC-controlled  fleet  consists  of  a  nucleus  fleet  of 
Navy-owned  cargo  ships  and  time  chartered  commercial  vessels.  These  are 
used  primarily  to  carry  breakbulk  cargo.  The  size  of  this  fleet  is 
calculated  from  the  cargo  movement  requirements  forecasted  by  the  shipper 
services.  MSC  uses  the  percentage  of  cargo  that  historically  was  shipped 
on  the  various  types  of  lift  within  the  controlled  fleet  and  then  allocates 
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that  percentage  of  the  forecasted  cargo  to  each  type  lift.  From  this,  con¬ 
tracts  are  entered  into  between  MSC  and  the  commercial  carriers  for  provi¬ 
sion  of  required  vessels.  If  the  forecasted  cargo  is  not  generated,  more 
shipping  is  available  than  required.  Since  ships  are  contracted  for  a 
minimum  of  1  year,  they  are  not  returned  until  the  end  of  the  contracted 
period  but  are  placed  in  reduced  operating  status.  During  FY  1982  the  MSC 
controlled  fleet  average  utilization  rate  was  40  percent.  Army  breakbulk 
forecast  for  FY  1981  was  555,000  measurement  tons  and  in  FY  1982  was 
453,000  measurement  tons  in  excess  of  actual  shipments.  This  equates  to  a 
peacetime  requirement  excess  of  a  minimum  of  four  breakbulk  cargo  vessels 
within  the  MSC  controlled  fleet. 

(2)  MSC  controlled  fleet  vessels  are  also  considered  by  MSC  when  cal¬ 
culating  its  contingency  shipping  requirements.  Therefore,  it  may  be  ex¬ 
pected  that  the  MSC  fleet  would  not  be  reduced  by  a  number  equal  to  the 
savings  generated  by  accurate  forecasts  and  that  the  cost  of  operating  con¬ 
tingency  vessels  would  be  absorbed  by  a  source  other  than  the  shipper  ser¬ 
vices.  An  excellent  study  on  MSC  controlled  fleet  sizing  is  contained  in 
the  OSD  Logistics  Systems  Analysis  Office  Dry  Cargo  Sealift  Forecast  Study, 
December  1982.1 

2-5.  OTHER  OBSERVATIONS.  During  the  study  of  the  Army's  forecasting  sys¬ 
tem,  several  factors  were  observed  which  have  major  impacts  on  the  prep¬ 
aration  and  the  use  of  the  forecasts  which  were  not  all  under  the  control 
of  the  Army.  These  include  mode  determination,  household  goods  shipment 
mode,  the  use  of  short-range  forecasts,  commodity  classification,  and  the 
MSC  billing  process. 

a.  Mode  Determination 

(1)  JCS  Pub  15  and  AR  55-30  require  that  forecasts  state  the  mode  or 
type  shipment  as  breakbulk,  container,  or  MILVAN.  However,  the  decision  as 
to  whether  cargo  is  transported  by  breakbulk  or  container  vessel  generally 
does  not  rest  with  the  Army  forecaster  or  shipper,  but  with  MTMC.  As  most 
breakbulk  cargo  is  carried  on  MSC  controlled  fleet  vessels  and  most  con¬ 
tainer  cargo  on  commercial  shipping,  this  decision  has  a  significant  impact 
on  forecast  accuracy.  As  stated,  MSC  uses  this  data  for  MSC  controlled 
fleet  sizing  and  in  the  development  of  commercial  cargo  rates.  Examination 
of  FY  82  forecasted  and  lift  data  shows  that  breakbulk  cargo  was  overfore¬ 
cast  by  453,000  MTON  or  24  percent  and  that  container  cargo  was  underforecast 
by  312,000  MTON  or  17  percent.  Table  2-7  shows  those  commodities  in  which 
offsetting  forecasted  versus  shipped  cargo  variances  occurred  in  FY  82. 
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Table  2-7.  Selected  FY  82  Breakbulk  vs  Container  Forecast  Performance 

(000  MTQN) 

Forecast  I  Actual  I  Error 


Breakbulk 


General 

297 

182 

115  over 

Household 

61 

32 

29  over 

POVs 

344 

319 

66  over 

Container 

General 

1,639 

1,830 

91  under 

Household 

39 

72 

13  under 

POVs 

138 

240 

102  under 

Commodity 


(2)  From  this  analysis  of  the  FY  82  forecast  performance,  it  appears 
that  a  significant  portion  of  the  forecasting  error  is  attributable  to  the 
breakbulk  versus  container  shipping  decision.  Continuation  of  the  require¬ 
ment  for  the  shipper  to  forecast  the  type  or  mode  of  shipment  will  most 
likely  result  In  continued  forecasting  problems  and  consequently  in  MSC 
controlled  fleet  sizing  errors.  Routing  of  annual  forecasts  through  MTMC 
to  MSC  would  allow  MTMC  and  MSC  to  determine  the  allocation  of  cargo  be¬ 
tween  the  various  modes. 

b.  Household  Goods  Shipments.  International  household  goods  shipments 
are  classified  into  six  categories,  one  of  which  (door-to-door  container, 
surface  code  5)  is  processed  by  MSC.  Door-to-door  container,  code  4,  is 
processed  totally  by  commercial  carrier  and  does  not  appear  in  MSC  billing 
information.  Currently,  code  5  shipments  are  less  than  3  percent  of  those 
shipped  under  code  4.  Minor  shifts  of  cargo  from  code  4  to  code  5  have  a 
major  Impact  on  the  accuracy  of  the  forecast  of  household  goods  to  be  pro¬ 
cessed  by  MSC.  As  the  forecasting  agencies  do  not  make  the  modal  decision 
for  household  goods  shipments,  the  accuracy  of  their  forecasts  is  not  under 
their  control. 


(1)  JCS  Pub  15  and  AR  55-30  direct  submission  by  the  shipper  services 
of  sealift  cargo  requirements  (short  range)  each  month  for  the  succeeding  3 
months.  Discussions  with  MTMC  personnel  indicate  that  ships  are  not  booked 
until  MTMC  has  been  notified  by  the  shipper  that  a  shipment  is  available 
and  that  short-range  forecasts  are  not  used  in  the  process.  Given  the 
nature  of  shipping  schedules.  It  Is  highly  unlikely  that  schedules  can  be 
firmly  arranged  if  only  the  month  and  area  of  shipment  are  known  and  not 
the  date  and  port. 


(2)  Evaluation  of  short-range  forecasts  indicate  that,  in  general, 
they  are  as  inaccurate  as  the  long-range  forecasts  and  for  some  commodities 
and  modes  they  are  much  less  reliable.  This  may,  in  part,  explain  their 
nonuse.  Table  2-8  shows  the  performance  of  short-range  forecasts  for  FY  81 
and  FY  82.  Long-range  error  is  shown  in  parentheses. 

(3)  Observations  of  most  commodity  shipping  history  reveals  monthly 
and  seasonal  trends  that  should  be  predictable  with  reasonable  accuracy  as 
components  of  the  long-range  forecast  when  progrm  information  is  included 
in  the  assessment.  Discussion  with  MSC  personnel  indicate  that  an  accurate 
long-range  forecast  and  periodic  change  reports  would  be  adequate  for  their 
use. 


(1)  The  special  cargo  commodity  contains  all  military  wheel  and  track 
vehicles  (including  sedans)  and  all  items  measuring  30  feet  in  any  dimension. 
Most  special  cargo  is  carried  by  the  MSC  controlled  fleet.  Privately  owned 
vehicles  are  a  separate  commodity  and  are  transported  by  breakbulk  and  sea- 
van  carrying  vessels.  As  there  is  no  physical  difference  between  POVs  and 
military  sedans  and  other  small  military  vehicles,  the  transportation  re¬ 
quirements  should  be  the  same.  Therefore  military  sedans  and  other  small 
vehicles  should  be  in  the  same  commodity  group. 


(2)  The  general  cargo  category  contains  a  significant  number  of  end 
items,  particularly  communications-electronics  equipment.  Current  Army 
transportation  account  code  (TAC)  structure  does  not  facilitate  the  identi¬ 
fication  of  these  items.  As  a  result,  general  cargo  requirements  contain 
elements  which  are  affected  by  program  decisions  and  for  which  impact  on 
forecasts  cannot  be  readily  assessed. 


(3)  The  empty  CONEX  and  assembled  aircraft  commodities  constitute  a 
small  segment  of  the  total  cargo  and  do  not  vary  significantly.  The  pre¬ 
vious  year's  lift,  extracted  from  the  billing  information,  should  be  ade¬ 
quate  forecasts  for  these  commodities  unless  program  information  indicated 
otherwise. 


e.  MSC  8i 11 inq  Process.  Analysis  of  several  years  of  MSC  billing  in¬ 
format  ionT^TsoisnorrwTth  personnel  involved  in  the  budget  process  and 
forecasting  personnel,  review  of  data  developed  by  US  Navy  Support  Command, 
and  information  received  at  Air  Force  Logistics  Center  indicate  that  the 
monthly  billings  received  from  MSC  usually  do  not  contain  all  bills  for  the 
reported  month  and  usually  contain  erroneous  entries.  Improvements  to  the 
MSC  billing  process  are  needed  if  the  services  are  to  accurately  capture 
monthly  and  seasonal  cargo  shipping  behavior  to  improve  their  forecast. 
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Table  2-8.  Short-range  Forecasting  Error  (000  MTON) 


Comnodlty 

FY  81 

FY  82 

Break bulk 

|  Container  j 

MILVAN 

Breakbulk  i 

Container  f 

MILVAN 

General 

+42( -286 ) 

+H0(  -126 ) 

-38( +10) 

+39( — 115 ) 

+248(+191) 

--(-4) 

Household 

— 3(— 3) 

-89(-29) 

+13(+13) 

-71 (-29) 

-20(+13) 

+5(+5) 

Special 

-♦(-92) 

♦15 (+1) 

— 15( -198) 

+7 ( -18) 

POV 

-26(-54) 

-12(+62) 

-41 (-66) 

+17 (+102) 

+1(+1) 

Reefer 

♦1(+1) 

+l(+5) 

+1(+1) 

♦9 ( +9 ) 

Anno 

+19(-45) 

♦3(+3) 

-14 (-8) 

-23 (-3) 

+3(+3) 

Cargo  trailer 

+36 (+22) 

-31(0) 

+S0( -11 ) 

-IO(-) 

CONEX 

+2C-3) 

-(-3) 

+2 ( +2 ) 

Coal/bulk 

-334(-85) 

♦2(+2) 

-108{ —53 ) 

(  )  Long-range  forecast  error. 

2-6.  SUMIARY  OF  OBSERVATIONS 

a.  General.  Observations  in  the  study  of  transporation  workload  fore¬ 
cast  system  were  categorized  into  three  principal  areas:  systemic— those 
aspects  involving  the  structure  and  interactions  within  it;  methodological- 
-those  aspects  relating  to  the  methodologies  used  by  various  forecasters; 
and  external  environmental— those  aspects  within  the  JCS  directed  forecast¬ 
ing  system  which  the  Army  does  not  control. 

b.  Systemic  Observations 

•  The  Army's  transportation  workload  forecasting  is  not  a  closed  loop 
system  and  can  be  classified  as  many  forecasters  performing  dupli¬ 
cative  functions  at  widely  separated  locations. 

•  The  long-range  and  short-range  forecasting  systems  are  not  indepen¬ 
dent  but  flow  on  separate  reporting  and  aggregation  channels. 

•  The  budget  and  cargo  forecast  systems  are  not  supportive  of  each 
other. 

•  Review  and  adjustment  of  forecaster  input  at  HQDA  and  LCA  is 
restricted  by  aggregation  of  inputs  and  lack  of  performance  data. 
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•  Army  historical  lift  and  forecast  files  are  limited  to  1  year  of 
lift  data. 

•  Present  allocation  of  TAC  does  not  facilitate  forecast  performance 
evaluation  or  identification  of  shippers'  programs,  or  specific 
sources  of  errors. 

•  Some  agencies  are  required  to  submit  cargo  movement  requirements 
for  commodities  and  programs  about  which  they  have  limited  infor¬ 
mation. 

•  The  link  between  the  various  programs  and  statements  of  their 
transportation  requirements  is  not  well  defined. 

•  The  present  operation  processes  of  the  MECHTRAM  system  are  not  ade¬ 
quate  for  accurate  forecast  or  budget  formulation. 

c.  Methodological 

•  There  are  no  standard  forecasting  methodologies  or  tools  in  use 
within  the  system. 

•  Extensive  subjective  forecasting  is  used  throughout  the  system. 

•  An  adequate  data  base  does  not  exist  within  the  Army  on  which  to 
base  forecasts. 

•  MSC  billing  information  is  the  most  complete  source  of  shipping 
data  but  does  not  currently  identify  shipper  and  program  in  suffi¬ 
cient  detail  to  develop  accurate  forecasts  for  all  commodities. 

•  Current  methodologies  generally  do  not  produce  accurate  forecasts. 

•  Those  agencies  using  several  years  of  historical  data  and  program 
knowledge  as  the  basis  for  their  forecasts  have  traditionally  been 
most  accurate. 

d.  External  Environment 

•  JCS  Pub  15  requires  the  Army  to  forecast  the  shipping  mode 
(container  vs  breakbulk),  but  the  decision  on  mode  is  generally 
made  by  MTMC. 

•  Household  goods  shipments  through  MSC  constitute  a  small  percentage 
of  total  household  goods  shipments.  The  accuracy  of  shipper  fore¬ 
casts  can  be  significantly  affected  by  shifts  between  modes.  This 
decision  is  not  made  by  the  shipper  service. 
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•  Indications  are  that  short-range  forecasts  are  not  used  for  port 
and  ship  scheduling  and  are  otherwise  of  questionable  value. 

•  The  MSC  billing  system,  which  is  the  primary  data  source,  is  slow 
and  may  contain  significant  inaccuracies. 

•  The  TAEDP  is  not  currently  an  accurate  source  of  information  for 
forecasting  of  equipment  distribution. 
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CHAPTER  3 

RECENT  STUDIES  OF  THE  TRANSPORTATION  WORKLOAD 
FORECASTING  PROCESS 


3-1.  GENERAL.  During  the  study  of  transportation  workload  forecasting, 
several  related  studies  and  reports,  which  were  conducted  during  the  pre¬ 
vious  10  years,  were  reviewed  to  minimize  duplication  of  effort  and  to 
benefit  from  data  already  available.  Most  recent  of  these  were  the  Dry 
Cargo  Sealift  Movement  Forecasting  and  Sizing  of  the  MSC  Fleet  Study,  con¬ 
ducted  by  the  OSD  Logistics  Systems  Analysis  Office  (LSAO)  in  December  1982 
for  the  Deputy  Assistant  Secretary  of  Defense  (L&MM);  the  Over-ocean  Cargo 
Forecasting  Procedures  Study,  conducted  by  the  DARCOM  Inventory  Research 
Office  (IRO)  in  September  1983;  and  an  MTMC  Forecasting  Methodology  (draft) 
1983. 

3-2.  OSD  LSAO  STUDY.  The  OSD  LSAO  Study  focused  on  breakbulk  cargo  fore¬ 
casting  only  and  its  impact  on  the  MSC-controlled  fleet  sizing  process. 

The  study  was  critical  of  the  Army  forecasting  process  and  highly  rec¬ 
ommended  both  the  current  Air  Force  and  developmental  MTMC  methodologies. 
They  concluded  that  a  wide  range  of  forecasting  methods  with  varying 
degrees  of  accuracy  were  now  in  use,  and  that  the  most  accurate  methods 
(USAF)  use  minimum  judgment,  are  mathematically  simple,  use  program  data, 
have  a  single  computation  point,  and  have  accurate  feedback  systems.  Their 
recommendations  include  a  proposal  that  the  DASD  (L&MM)  direct  the  services 
to  develop  forecasting  methodologies  containing  the  characteristics  as  out¬ 
lined  above  and  contained  in  the  USAF  system.  As  stated,  the  LSAO  Study 
only  addressed  breakbulk  and  not  container  cargo.  This  omitted  consider¬ 
ation  of  the  relationship  between  the  cargo  shipping  modes.  Further,  it 
did  not  address  (as  previously  discussed)  the  function  of  MTMC  in  deter¬ 
mining  the  type  of  shipment  or  mode.  The  study  focused  on  FY  81  data  only. 

3-3.  THE  DARCOM  STUDY.  The  DARCOM  IRO  has  recently  completed  a  study  of 
the  procedures  now  in  use  within  DARCOM.  Their  observations  paralleled 
many  of  those  contained  in  this  study  and  were  developed  in  coordination 
with  the  study  team.  They  recommended  a  new  mathematical  forecasting  model 
(Winters),  an  expansion  of  the  number  of  TACs  to  identify  program  items, 
that  more  effort  be  given  to  producing  an  accurate  long-range  forecast,  and 
that  the  use  of  the  TAEDP  be  evaluated. 


3-4.  MTMC  DEVELOPMENTAL  FORECASTING  METHODOLOGY.  HQ  MTMC  is  developing  a 
methodology  which  would  allow  MTMC  to  prepare  improved  transporation  work¬ 
load  forecasts.  Initial  efforts  were  to  create  and  evaluate  a  4-year  aver¬ 
age.  This  was  later  revised  to  include  a  smoothing  technique.  The  OSD  Dry 
Cargo  Sealift  Movement  Forecasting  Study  recommends  that  the  MTMC  model  be 
used  to  evaluate  the  service  forecasts. 


CHAPTER  4 

OTHER  SERVICE  FORECASTING  SYSTEMS 


4-1.  GENERAL.  USAF  and  USN  transportation  workload  forecasting  systems 
were  investigated  to  determine  the  methodologies  used  and  their  efficacy, 
and  to  identify  elements  which  had  potential  application  to  the  solution  of 
the  Army  forecasting  problem. 

4-2.  USAF  FORECASTING  SYSTEM 

a.  USAF  forecasts  of  transportation  requirement  procedures  are  con¬ 
tained  in  AFR  75-15,  dated  6  June  1979  (currently  being  revised).  Head¬ 
quarters,  USAF  Logistics  Command  (AFLC)  has  responsibility  for  forecasting 
and  budgeting  for  Air  Force  cargo  requirements.  AFLC  automated  data  pro¬ 
cessing  procedures  are  contain  in  AFLC  Regulation  AFL171-125  dated  15 
November  1979.  This  regulation  provides  a  detailed  operational  description 
of  the  system  and  its  associated  environment. 

b.  The  Air  Force  long-range  forecast  is  developed  from  a  bivariate  re¬ 
gression  analysis  of  projected  flying  hours  with  an  8-year  data  base  of 
actual  lift  history  which  has  been  developed  from  MSC  billing  information. 
Additionally,  designated  commands  are  required  to  report  nonrecurring  or 
unusual  requirements  which  are  added  to  the  forecast.  In  those  cases  where 
a  commodity  constitutes  a  small  percentage  of  the  total  forecast,  the  pre¬ 
vious  year's  actual  shipment  is  used  for  the  forecast.  Mode  of  shipment  is 
based  on  the  percentage  of  the  total  shipment  sent  by  the  particular  modes 
during  the  previous  year. 

c.  Short-range  forecasts  are  prepared  by  a  triple  exponential  smoothing 
of  the  8-year  data  base. 

d.  USAF  TACs  are  structured  to  identify  both  the  program  and  the  ship¬ 
ping  customer. 

e.  The  USAF  system  includes  effective  procedures  to  validate  MSC 
billing  information. 

4-3.  US  NAVY  FORECASTING  SYSTEM 

a.  Forecasting  of  US  Navy  sealift  requirements  is  governed  by 
NAVSUPINST  4620. 7D.  Navy  Material  Transport  Office  (NAVMTO),  Norfork, 
Virginia,  collects  and  submits  both  long-range  and  short-range  forecasts 
for  the  Naval  Support  System  Command,  which  is  responsible  for  the  fore¬ 
casting  system. 


b.  Navy  cargo  transportation  requirements  forecasting  is  currently  a 
manual  process  which  uses  historical  data,  major  program  requirements,  and 
forecaster  judgment.  Sealift  requirements,  excluding  personal  property. 
Coast  Guard  items,  construction  material,  and  aircraft,  are  reported  only 
as  an  exception  to  the  previous  year's  lift. 

c.  Navy  over-ocean  cargo  movement  requirements  are  approximately 
one-third  of  those  required  by  the  Army  and,  as  a  consequence,  their  errors 
have  less  impact.  Their  greatest  diffculties  are  determining  the  shipping 
mode  and  the  route  to  be  used.  Additionally,  the  delay  in  receiving  com¬ 
plete  actual  lift  and  billing  data  from  MSC  complicates  both  forecasting 
and  budget  execution. 

d.  The  Navy  has  approved  a  program  to  automate  their  forecasting 
system.  The  objectives  of  the  automated  system  are  to  maintain  historic 
lift  and  forecast  files,  and  generate  all  required  forecasts.  It  is 
expected  that  the  system  will  incorporate  the  essential  features  currently 
used  in  the  USAF  system. 

4-4.  APPLICATION  TO  ARMY  FORECASTING.  The  USAF  system,  with  the  exception 
of  the  forecasting  routines,  is  an  appropriate  paradigm  on  which  to  base  a 
revised  Army  forecasting  system.  The  USAF  system  contains  all  the  com¬ 
ponents  required  for  forecasting,  cost  and  performance  reporting. 
Transportation  Operating  Agency  (TOA)  bill  processing,  and  data  base 
maintenance. 
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CHAPTER  5 

ALTERNATIVE  NETHOOOLOGIES  CONSIDERED  TO  IMPROVE  THE  FORECASTING  SYSTEM 


5-1.  GENERAL.  The  requirements  of  the  forecasting  system,  the  sources  of 
information,  and  the  systemic,  external  environment,  and  methodological 
aspects  of  the  forecasting  system  were  considered  in  the  development  of 
proposals  for  changes  to  the  current  projects.  Systemic  and  external  en¬ 
vironment  factors  are  discussed  in  this  chapter.  Forecasting  methodologies 
are  discussed  in  Chapter  6. 

a.  System  Requirements.  To  satisfy  the  essential  requirement  of  JCS 
Pub  15,  AR  55-23,  AR  55-30,  and  an  accurate  cargo  forecasting  process,  the 
systems  should  provide  the  following  capabilities: 

(1)  Forecasts  of  cargo  movement  requirements  by  commodity,  route, 
mode,  and  program. 

(2)  Forecasts  of  budget  requirements. 

(3)  System  performance  data  by  commodity,  route,  mode,  and  program. 

(4)  A  historical  data  base  of  cargo  movement  requirements. 

b.  Information  Sources.  Cargo  which  generates  peacetime  over-ocean 
movement  requirements  is  generally  in  support  of  personnel  and  equipment 
deployed  overseas  or  in  support  of  specific  programs  such  as  force  moderni¬ 
zation.  As  peacetime  overseas  deployed  personnel  and  equipment  densities 
do  not  vary  significantly,  historical  cargo  movement  requirements  and  cur¬ 
rent  program  information  contain  most  of  the  data  required  for  forecast 
development.  The  MSC  cargo  billing  is  the  most  comprehensive  and  could  be 
the  most  accurate  source  of  historical  lift  data.  Sources  of  program  data 
are  (1)  personnel  related  programs— MILPERCEN  and  ODCSPER,  (2)  materiel 
related  programs— DARCOM  and  DCSRDA,  (3)  fuels  (coal  and  coke)— 0DCSL0G, 
and  (4)  AAFES  and  MPSA  for  their  specific  programs. 

5-2.  ALTERNATIVE  SYSTEM  STRUCTURES 

a.  General .  Several  system  structures  were  examined  considering  the 
requirements  of  the  system,  the  sources  of  the  information,  and  the  costs 
in  personnel  and  facilities.  In  each  case  it  was  assumed  that  (1)  improved 
forecasting  models  or  procedures  were  available,  (2)  TACs  were  restructured 
to  ensure  identification  of  force  modernization  items  and  the  specific  gen¬ 
erators  of  general  and  special  cargo,  and  (3)  that  the  present  forecasting 
requirements  will  continue.  There  are  many  possible  system  structures 
which  can  satisfy  the  forecasting  requirements  to  some  degree.  Of  these, 
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four  predominate:  (1)  the  current  system  with  improved  methodology,  (2) 
the  current  system  with  the  long-range  and  short-range  systems  combined, 

(3)  forecasting  by  a  single  Army  agency  using  historical  program  data  from 
DARCOM,  DCSLOG/DLA  and  MILPERCEN,  (4)  forecasting  by  MTMC  and  other  TOA 
with  program  input  from  the  Army.  The  first  two  represent  decentralized 
options  while  the  latter  two  are  centralized  options. 

b.  Current  Forecasting  System  Retained.  This  alternative,  which  is 
portrayed  in  Figure  5-1,  envisions  the  retention  of  the  present  forecasting 
agencies  and  channels,  correction  of  current  forecasting  anomalies,  and 
improved  forecast  and  budget  coordination.  It  assumes  that,  given  an  im¬ 
proved  methodology  and  restructured  TACs,  improvements  in  forecast  accuracy 
could  be  achieved  within  the  framework  of  the  present  system.  This  option 
requires  that  each  forecasting  agency/command  be  provided  (or  develop)  an 
accurate  mathematical  forecasting  model  capability  an  accurate  data  base, 
and  that  the  forecasters  have  access  to  program  data.  The  TAC  structure  is 
such  that  it  allows  each  forecasting  agency  to  identify  the  amount  of  each 
commodity  that  it  has  shipped  and  the  program  responsible  for  the  shipment. 
The  time  sequencing  and  frequency  of  reports  remain  as  required  by  the  cur¬ 
rent  procedures.  Reports  stating  forecasted  tonnage  by  route,  commodity, 
and  program  are  prepared  using  a  directed  forecasting  methodology  for  all 
commodities  except  for  nonrecurring  or  program  specific  items.  Program 
data  is  derived  from  the  most  appropriate  sources.  These  sources  include 
TAEDP,  Materiel  Management  Reports  (MMR),  the  Committee  on  Ammunition 
Logistics  Support  Report,  and  the  Annual  Solid  Fuels  Requirement  Report. 

The  final  annual  forecast  and  each  short-range  forecast  are  updated  ver¬ 
sions  of  the  previously  submitted  forecast  including  the  effects  of  program 
changes.  Also,  in  this  option,  the  LCA  system  and  MECHTRAM  are  modified  to 
identify  forecasts  which  are  significantly  different  from  expected  values 
and  periodically  compare  performance  reports  with  individual  command/agency 
forecasts.  Significant  variances  are  resolved  between  the  forecasting 
agency  and  HQDA  for  long-range  forecasts  and  with  LCA  for  the  short-range 
forecasts.  The  monthly  performance  reports  are  segregated  at  LCA  and  for¬ 
warded  to  the  forecasting/shipper  agencies  for  comparative  analysis  and 
retention  for  data  base  purposes. 
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Systems  Combined.  This  alternative,  Figure  5-2,  requires  that  the  same 
agencies  as  in  the  previous  alternative  prepare  and  submit  forecasts,  ex¬ 
cept  that  both  the  long-range  and  short-range  forecasts  are  submitted  to 
one  agency  versus  two  for  consolidation  and  adjustment.  This  agency  could 
be  either  HQDA,  HQ  DARCOM,  LCA,  or  any  other  designated  agency.  All  other 
processes  are  as  in  the  previous  alternative.  Consolidation  of  the  two 
procedures  provides  a  more  positive  linkage  between  the  long-range  and 
short-range  forecasts. 


Figure  5-2.  Current  Forecasting  System  Modified 


d.  Single  Agency  Forecasting  with  Program  Input.  The  alternative. 
Figure  5-3,  envisions  a  single  agency  preparing  the  forecasts  for  all  com¬ 
modities,  routes,  and  programs  based  on  forecasting  model  outputs.  These 
are  modified  by  current  information  on  specifically  identified  programs 
such  as  coal  and  coke,  ammunition,  and  force  modernization  items.  Program 
data  is  submitted  by  AAFES,  MPSA,  DARCOM,  and  DCSLOG  annually,  stating  pro¬ 
jected  monthly  shipping  requirements  by  item,  origin,  and  destination.  All 
other  commands  report  projected  or  nonrecurring  requirements.  Negative 
reports  are  required.  Projected  annual  personnel  movements  are  provided 
from  MILPERCEN/DCSPER  to  the  forecasting  agency.  Where  the  personnel  move¬ 
ment  forecast  is  significantly  different  from  the  previous  year's,  the 
model  generated  forecasts  for  POV  and  household  goods  (if  forecast)  ship¬ 
ments  are  modified  to  reflect  the  expected  movements.  Agencies  submitting 
annual  program  information  also  forward  quarterly  reports  indicating 
changes  to  the  previously  submitted  annual  reports.  Other  change  reports 
would  be  submitted  when  shipments  on  a  given  route  are  expected  to  differ 
from  the  forecast  by  600  MTON.  Monthly  performance  reports  are  analyzed  by 
the  forecasting  agency  to  evaluate  the  accuracy  of  the  forecasts.  The  re¬ 
ports  are  used  to  Identify  and  resolve  discrepancies  between  forecasts  and 
actual  lift,  and  to  modify  forecasting  models.  The  forecasting  activity  in 
this  option  requires  an  automated  data  processing  system  which  provides  for 
those  capabilities  currently  in  the  MECHTRAM  system.  It  provides  a  route 
and  system  performance  monitoring  process  such  as  that  employed  by  USAF  and 
incorporates  mathematical  forecasting  capabilities  developed  in  this  study. 
MSC  bills  received  by  the  forecasting  agency  are  validated.  They  are  then 
used  as  source  data  for  forecast  performace  evaluation  and  for  the  data 
base- to  be  used  in'future  forecast  formulation.  Activities  of  the  agency 
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would  consolidate  all  actions  now  performed  in  the  MECHTRAM  system  at  HQDA 
and  the  AOP  support  to  forecasting  at  LCA.  Personnel  in  the  forecasting 
all  should  be  of  military  rank  04/03,  MOS  49,  or  their  civilian  equivalent 
who  possess  a  knowledge  of  the  Defense  Transportation  System  and  computer 
operations. 


Figure  5-3.  Single  Agency  Forecasting 


e.  Forecasting  by  the  TOA  (Objective  System).  In  this  alternative. 
Figure  5-4,  MTMC  prepares  the  annual  forecasts  based  on  cargo  movement 
history  and  program  information  provided  by  the  shipper  services.  HQDA 
consolidates,  program  information  provided  by  DARCOM,  MILPERCEN,  and  the 
Army  Staff  and  forwards  it  to  the  MTMC.  MTMC  computes  cargo  requirements 
based  on  forecasting  model  output  and  program  information,  provides  the 
forecast  to  MSC,  and  from  this,  transportation  costs  are  estimated.  The 
TOA  provide  HQDA  with  forecasted  tonnages  and  costs  for  service  budget 
preparation  purposes.  Transportation  costs  are  paid  as  under  the  present 
system.  Budget  performance  reporting  and  validation  of  billing  information 
functions  are  performed  by  the  MECHTRAM  or  similiar  system.  MECHTRAM  also 
produces  monthly,  and  as  required,  cost  and  performance  reports  for  the 
various  HQDA  staff  elements  (DCSLOG,  DCSPER).  As  the  historical  shipping 
data  used  in  forecasting  is  originated  from  TOA  billing  and  the  TOAs  are 
supported  by  extensive  computing  facilities,  it  is  feasible  and  desirable 
for  the  TOAs  to  develop  the  forecasts.  In  particular,  HQ  MTMC  is  capable 
of  producing  forecasts  of  over-ocean  cargo  movement  requirements  for  both 
MTMC  and  MSC  as  a  result  of  its  own  forecasting  model  development  which  has 
occurred  concurrent  with  this  study.  Model  parameters  developed  by  this 
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study,  and  discussed  in  the  next  section,  can  be  provided  to  MTMC  to  assist 
in  forecast  development. 


Figure  5-4.  Forecasting  by  TOA 


5-4.  EVALUATION  OF  ALTERNATIVE  SYSTEMS.  Table  5-1  portrays  the  identifi¬ 
able  direct  cost  associated  with  the  four  alternatives  previously 
described.  This  section  discusses  the  four  alternatives  in  terms  of 
personnel  and  operational  costs. 


Table  5-1.  Alternative  System  Costs 


(1) 

(2) 

(3) 

(4) 

Current 

system 

Current 

system 

modified 

Single  Army 
forcasting 
activity 

Forecast  by 

TOA 

(objective  system) 

Man-years 

14 

13 

6 

5 

Models 

17 

16 

1 

1 

Data  bases 

17 

16 

1 

1 
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a.  Alternative  1,  the  present  system,  requires  at  least  23  personnel. 

Of  these,  9  are  full-time  forecasters.  The  efforts  of  the  other  14  range 
from  10-50  percent  of  their  time.  This  effort  equates  to  approximately  14 
man-years  annually.  By  eliminating  the  dual  reporting  system  as  in  Alter¬ 
native  2,  an  estimated  one  space  would  be  saved  as  the  result  of  combining 
the  HQDA-LCA  effort.  Largest  savings  of  personnel  occur  in  Alternatives  3 
and  4.  Alternative  3  envisions  a  three-person  forecasting  cell,  a  HQDA 
program  officer,  a  DARCOM  program  coordinator,  and  part-time  effort  by  one 
person  at  AAFES  and  MPSA.  Alternative  4  would  reduce  the  size  and  effort 
of  the  Army  forecasting  cell  and,  if  adapted  by  all  shipper  services,  would 
result  in  similar  personnel  savings  by  each  shipper  service. 

b.  In  addition  to  the  direct  costs  in  personnel,  already  discussed, 
there  are  significant  indirect  costs  associated  with  the  forecasting  ac¬ 
tivities.  These  include  management,  clerical,  and  communications  support. 
Currently,  17  separate  commands  or  agencies  and  HQDA  prepare  forecasts. 
Additionally,  several  of  these  agencies  consolidate  forecasts  prepared  by 
subordinate  commands,  e.g.,  LCA,  a  forecasting  agency,  also  consolidates 
forecasts  prepared  by  five  DARCOM  commands.  This  would  be  reduced  signifi¬ 
cantly  in  Alternatives  3  and  4.  Reduction  of  the  number  of  agencies,  as  in 
Alternatives  3  and  4,  also  reduces  lead  time  required  for  forecasts,  thus 
allowing  more  recent  data  to  be  considered  when  they  are  developed. 

c.  Accurate  forecasting  and  use  of  a  predictive  model  requires  each 
forecasting  agency/command  to  maintain  an  accurate  historical  record  of 
cargo  forecasted  and  cargo  shipped.  Most  forecasting  models  require  a  min¬ 
imum  of  50  data  points.  This  results  in  a  data  base  at  least  50  months  of 
shipping  data.  Current  estimates  of  MSC  records  of  Army  shipments  are 
30,000  records,  each  containing  140  characters  per  year.  Alternative  1  has 
a  complete  data  base  at  MSC/MTMC,  HQDA,  and  LCA.  Each  forecasting  agency 
in  turn  maintains  the  data  relevant  to  its  operations.  Alternative  2  would 
eliminate  either  the  HQDA  or  LCA  data  base.  Alternative  3  requires  only 
one  Army  data  base  and  Alternative  4  requires  only  the  retention  of  the 
MSC/MTMC  data  base,  which  is  the  source  of  data  for  all  other  data  bases  in 
the  other  alternatives. 

d.  Accurate  forecasting  requires  the  use  of  some  form  of  predictive 
model.  Alternatives  1  and  2  require  the  development  and  maintenance  of  a 
large  number  of  models,  as  it  is  unlikely  that  a  general  model  would  be 
applicable  to  each  commodity  and  each  shipper.  This,  combined  with  the 
maintenance  of  the  data  bases,  requires  a  significantly  greater  effort  than 
Alternatives  3  and  4.  (Specific  forecasting  models  are  discussed  in 
Chapter  6.) 

e.  Forecasting  systems  require  accurate  and  timely  information  on  their 
programs.  All  agencies  which  currently  forecast,  with  the  exception  of  HQ 
FORSCOM  and  District  Engineer,  New  York,  have  access  to  the  program  infor¬ 
mation  which  is  required  to  formulate  the  cargo  movement  requirement  for 
which  they  are  responsible.  However,  aggregated  program  information  is 
available  at  HQDA/HQ  DARCOM  to  determine  most  materiel  associated  cargo 
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movement  requirements  and  at  HQOA/MILPERCEN  to  determine  personnel  asso¬ 
ciated  cargo  movement  requirements.  As  AAFES  and  MPSA  have  other  service 
responsibilites,  these  agencies  alone  have  cognizance  of  the  programs  which 
will  affect  their  requirements  for  cargo  movement.  Historic  cargo  lift 
data  is  available  from  the  MSC  and  MTMC  billing  tapes  at  HQDA,  LCA,  and  to 
a  limited  extent  at  some  of  the  forecasting  agencies.  Options  1  and  2  re¬ 
quire  the  processing  of  program  Information  and,  as  discussed  earlier,  the 
NSC  billing  information  by  all  forecasting  commands/agencies.  Options  3 
and  4  use  consolidated  program  information  and  TOA  bills.  They  convert 
program  and  historic  information  into  forecasts  at  a  single  location.  The 
quality  of  the  program  Information  available  to  forecasters  in  any  of  the 
options  should  be  comparable.  However,  the  translation  of  the  data  into 
cargo  movement  requirements  would  appear  to  be  more  efficient  if  performed 
at  a  single  location  as  in  Options  3  and  4. 

5-4.  COLLOCATION  OF  CARGO  MOVEMENT  REQUIREMENTS  AND  BUDGET  FORECASTING 
ACTIVITIES 

a.  Discussions  during  the  early  stages  of  the  study  surfaced  a  concern 
that  an  underlying  cause  of  current  forecasting  inaccuracies  is  the  sepa¬ 
ration  of  budget  and  cargo  movement  requirements  forecasting.  Investigation 
of  the  forecasting  system  did  not  confirm  this  as  a  condition  or  source  of 
error.  From  a  management  perspective,  if  the  budgets  are  the  estimated 
costs  of  shipping  the  cargo,  then  the  budgets  should  be  developed  from  the 
forecast  of  the  cargo  movement  requirements  and  the  projected  cost  schedule. 
This  does  not  imply  that  the  cargo  movement  and  budget  forecasts  should  be 
performed  by  the  same  office  but  does  imply  that  a  procedure  exists  which 
ensures  coordination  and  Interdependence  of  the  two  processes.  The  docu¬ 
mentation  of  the  MECHTRAM  process  describes  this  function  as  occurring  at 
HQDA;  however,  as  discussed  earlier,  the  current  budget  development  process 
does  not  use  the  MECHTRAM  reports. 

b.  USAF  and  USN  budget  systems  were  examined  to  determine  what  effects, 
if  any,  could  be  attributed  to  the  separation  or  collocation  of  cargo  fore¬ 
casting  and  budget  preparation  processes.  USAF  budget  development  and  exe¬ 
cution  and  cargo  movement  requirements  forecasting  are  all  performed  by  a 
single  office  at  USAF  Logistics  Command.  USN  has  separated  the  functions. 
Budget  processes  are  performed  by  Naval  Supply  Systems  Command,  Arlington, 
Virginia,  while  cargo  movements  forecasting  is  performed  by  the  Navy  Material 
Transport  Office  (NAVMTO),  Norfolk,  Virginia.  The  systems  are  equally  accu¬ 
rate.  No  other  systems  were  identified  wherein  the  separation  of  operations 
forecasting  and  fiscal  forecasting  made  a  difference  to  the  accuracy  of  the 
operations  forecasting,  and  no  alternatives  were  developed  which  varied 
budgeting  responsibility. 
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CHAPTER  6 

METHODOLOGY  DEVELOPMENT 


6-1.  GENERAL.  The  purpose  of  this  section  is  to  describe  the  development 
of  the  forecasting  methodologies  used  to  model  the  transportation  workload 
forecasting  (TWF)  data.  Paragraph  6-2  describes  the  data  gathering  and 
data  manipulation  processes.  Paragraph  6-3  discusses  time  series  analysis 
and  the  reasons  for  using  this  type  of  data  analysis.  Paragraph  6-4 
provides  a  general  description  of  the  Box-Jenkins  modeling  approach  for 
time  series  data.  (A  detailed  description  of  this  approach  is  discussed  in 
Appendix  E.)  Paragraph  6-5  discusses  the  Winters  approach  to  modeling  time 
series  data.  (A  detailed  description  of  the  Winters  method  is  discussed  in 
Appendix  G.)  Paragraphs  6-6  thru  6-14  discuss  the  individual  commodities 
that  were  forecasted,  the  processes  involved  in  developing  the  mathematical 
models,  and  the  forecasting  results.  Figures  6-2  through  6-10,  presented 
later,  compare  the  commodity  forecasts  of  the  Box-Jenkins  and  Winters 
models  with  the  actual  lift  data  for  FY  82.  Appendix  F  details  specific 
analytic  steps  of  Box-Jenkins  model  building  for  each  commodity.  Appendix 
H  contains  similar  information  for  Winters  models.  Paragraph  6-15 
summarizes  the  results  of  the  forecasting  models. 

6-2.  DATA  DESCRIPTION.  This  paragraph  details  the  steps  taken  to  develop 
the  TWF  data  base:  (1)  data  acquisition,  (2)  computer  system  conversion, 
(3)  data  format  and  condition,  and  (4)  data  reduction. 

a.  Data  Acquisition.  Six  years  of  lift  information  (tonnage)  on  Army 
cargo  shipments  was  acquired  from  MSC.  The  data  provided  by  MSC  contained 
lift  information  by  route,  commodity,  and  mode  of  shipment. 

b.  Computer  System  Conversion.  The  following  specifications  were  fol¬ 
lowed  during  tape  creation: 


Tracks 

Code 

File(s) 

Density 

Logical  record  length 

Block  size 

Label 


9 

ASCII 

Sequential,  unformatted 
1,600  BPI 
80 

3,200 

None 


c.  Data  Format  and  Condition.  Each  record  contained  35  columns  of 
data,  formatted  shown  in  Table  6-1. 
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Table  6-1.  Data  Format 


1-2 

3-4 

5-6 

7 

8-10 

11-12 

13-15 

16-17 

18-21 


Fiscal  year 

Sail  year 

Sail  month 

Container  type 

Port  of  embarkation  (POE) 

POE  traffic  area 

Port  of  debarkation  (POD) 

POD  traffic  area 

TAC 


22-23  Commodity 

24-35  Measurement  tons 


Initial  attempts  to  read  the  data  with  a  univariate  statistics  program 
failed,  due  to  the  existence  of  nonnumeric  entries  in  data  fields.  Records 
containing  aberrative  subfields  were  corrected  or  deleted,  as  appropriate. 
Records  indicating  zero  shipment  of  measurement  tons  were  removed.  This 
process  reduced  the  original  data  base  of  184,645  records  to  a  final  count 
of  150,387  records  covering  the  timeframe  October  1977  through  June  1983. 

d.  Data  Reduction.  Several  special  purpose  FORTRAN  programs  were  de¬ 
veloped  in  order  to  transform  the  revised  data  base  into  forms  suitable  for 
analysis. 

(1)  In  preparation  for  preliminary  time  series  analysis,  the  data 
base  was  broken  out  into  12  general  commodity  files,  each  of  which  con¬ 
tained  records  of  shipment  tonnage  versus  time. 

(2)  The  12  (raw)  commodity  files  were  segregated  into  various  route 
subfiles  and  7,202  "credits"  discovered  in  the  original  data  base  were  re¬ 
moved.  The  credits  had  been  entered  into  the  original  data  base  in  order 
to  offset  certain  debits. 

(3)  Time  series  data  were  developed  for  424  unique  routes,  and  the 
routes  were  rank -ordered  according  to  decreasing  levels  of  shipment 
activity. 

(4)  Finally,  tonnage  was  segregated  on  a  monthly  basis  and  according 
to  chronological  shipment  date.  One  program  was  designed  to  process  routes 
individually.  The  other  was  designed  to  form  composites  of  all  routes  for 
each  commodity. 
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e.  Time  Series.  Aggregated  time  series  data  for  each  commodity  are 
contained  in  Figures  G-l  through  G-12,  Appendix  G.  Of  the  81  months  of 
data,  60  were  deemed  necessary  for  model  construction  and  adequacy 
verification.  It  is  generally  recommended  that  at  least  50  observations  be 
used,  if  possible,  when  developing  Box-Jenkins  time  series.  Specifically, 
data  from  October  1977  through  September  1981  were  used  to  build  both  the 
Box-Jenkins  and  Winters  models. 


f.  Forecast  Interval.  The  12  observations  from  the  sixth  year  of  data 
(FY  82)  were  used  as  the  basis  for  testing  the  forecasting  accuracy  of  each 
model.  Two  measures  of  forecasting  accuracy  were  selected  from  among 
several  potentially  acceptable  alternatives: 

(1)  Annual  Percentage  Error.  The  annual  percentage  error  between 
forecasted  and  actual  cargo  shipments  was  used  as  a  simple  aggregate 
measure  of  relative  accuracy  from  a  practical,  point  of  view. 

(2)  Mean  Square  Error  (MSE).  The  mean  squared  error  of  the  12 
monthly  observations  was  used  to  capture  the  absolute  accuracy  from  a 
statistical  standpoint. 

6-3.  TINE  SERIES  ANALYSIS 

a.  Background.  Time  series  analysis  is  the  application  of  analytical 
techniques  to  model  historical  data.  It  differs  from  most  types  of  histor¬ 
ical  data  analysis  in  that  no  attempt  is  made  to  determine  the  factors  af¬ 
fecting  the  behavior  of  the  historical  trend.  Another  difference  is  that 
time  series  analysis  is  concerned  with  historical  data  that  are  serially 
correlated  and  not  independent.  Time  series  analysis  attempts  to  model 
behavioral  trends  and  relationships  between  consecutive  observations. 

b.  Analysis  Considerations.  Time  series  analysis  techniques  were  con¬ 
sidered  the  most  apropriate  analytic  tools  to  forecast  the  US  Army 
transportation  workload  requirements  after  considering  three  important 
factors:  (1)  data  pattern,  (2)  time  and  efficiency,  and  (3)  accuracy  of 
forecasts. 

(1)  Time  series  analysis  is  usually  the  analysis  of  a  single  variable 
over  time,  i.e.,  univariate.  In  many  cases,  the  univariate  series  exhibits 
a  behavioral  pattern  over  time.  The  pattern  may  be  seasonal  or  a  persis¬ 
tent  trend  may  exist  in  the  data.  This  behavior  is  usually  made  up  of  de¬ 
pendent  observations  and  unlike  most  analytic  techniques,  time  series  anal¬ 
ysis  does  not  rest  on  the  assumption  of  independent  observations.  Instead, 
time  series  analysis  capitalizes  on  the  dependency  of  the  observations, 
determines  the  patterns  in  the  data  and  then  models  the  patterns.  The  data 
contained  in  the  TWF  study  was  highly  dependent,  therefore  time  series 
analysis  was  selected  as  the  most  appropriate  statistical  modeling  tool. 

(2)  Due  to  the  complexity  of  events  that  produced  the  transportation 
lift  history,  it  was  virtually  impossible  to  determine  all  the  factors  that 
affected  the  historical  process.  In  many  cases  the  factors  affecting  the 


process  could  not  be  defined  quantitatively  or  a  causal  relationship  be¬ 
tween  the  factors  and  the  process  could  not  be  determined.  Sufficient 
proxies,  such  as  troop  density,  could  have  been  developed  to  resolve  these 
issues,  however,  the  time  and  energy  required  to  obtain  statistically  rele¬ 
vant  factors  would  have  been  immense.  Time  series  analysis  enables  quick 
development  of  forecasting  models  since  factors  affecting  the  process  are 
for  the  most  part  ignored.  Several  studies  which  compared  simple  time 
series  models  to  complex  econometric  models  indicated  that  time  series 
models  can  provide  comparable  forecasting  results  and,  in  some  cases, 
better  forecasting  results  than  the  complex  models. 

(3)  Finally,  as  indicated  earlier,  time  series  models  are  primarily 
concerned  with  identifying  systematic  behavioral  patterns  in  the  data. 

Once  the  systematic  pattern  is  identified,  the  data  can  be  separated  into 
two  components:  (1)  the  data  pattern  and  (2)  the  random  error  of  the  data 
pattern.  This  statement  is  expressed  as: 

DATA  =  PATTERN  +  RANDOMNESS 

If  the  future  observations  of  a  time  series  can  be  predicted,  without  any 
error,  then  a  deterministic  model  is  sufficient  because  the  randomness  of 
the  data  is  negligible.4  If  the  data  exhibit  random  behavior  over  time, 
the  predictions  of  future  observations  are  mainly  concerned  with  the 
identification,  modeling,  and  prediction  of  the  error  term.  Time  series 
analysis  is  the  best  forecasting  methodology  to  employ  when  the  accurate 
modeling  of  the  error  term  is  critical.  TWFS  data  exhibited  a  random 
behavior  over  time;  thus,  stochastic  time  series  analysis  was  determined  to 
be  the  most  appropriate  modeling  approach. 

6-4.  BOX-JENKINS  MODELING  APPROACH 

a.  Background.  Box-Jenkins  models  are  a  unique  set  of  linear  time 
series  models  used  to  model  stochastic  time  series  data.  Box-Jenkins 
models  fall  into  three  classes:  autoregressive  (AR),  moving  average  (MA), 
and  mixed  (ARMA).  Box-Jenkins  models  find  their  origin  in  the  AR  models 
that  were  first  introduced  by  Yule  (1926)  and  later  generalized  by  Walker 
(1931).  MA  models  were  first  developed  by  Slutzky  (1937),  and  ARMA  models 
were  initially  theorized  by  the  work  of  Wold  (1938).  George  Box  and  Gwelym 
Jenkins  are  responsible  for  collating  these  previous  works  and  establishing 
an  approach  to  apply  these  models. 2  The  Box-Jenkins  approach  consist  of 
three  steps  (see  Figure  6-1): 

(1)  Identification  -  the  first  step  in  applying  the  Box-Jenkins 
methodology  was  to  identify  the  degree  of  homogeneity  in  the  data,  i.e., 
how  many  times  the  series  must  be  differenced  to  achieve  stationarity. 

Once  the  degree  of  homogeneity  was  established,  the  data  pattern  was  iden¬ 
tified  as  AR,  MA,  or  ARMA. 


(2)  Estimation  -  After  the  data  pattern  was  correctly  identified, 
parameter  estimates  for  AR,  MA,  or  ARMA  models  were  generated  to  obtain  a 
model  that  best  fit  the  data. 

(3)  Verification  -  Finally,  a  test  run  using  the  estimated  model 
parameters  was  conducted.  The  results  and  diagnostic  checks  were  performed 
on  the  model  parameter  estimates  and  residual  estimates  to  ensure  goodness 
of  fit  and  adequacy  of  the  model.  The  predictive  value  of  the  model  was 
evaluated  and  analyzed  using  historical  data  that  was  not  used  to  develop 
the  original  model.  If  the  model  was  not  adequately  verified,  then  steps 
1-3  were  repeated  until  an  appropriate  model  was  identified. 


Figure  6-1.  Box -Jenkins  Modeling  Approach 


b.  Stationarity.  The  most  crucial  element  in  applying  Box-Jenkins 
models  is  the  principle  of  stationarity.  Stationary  data  are  defined  as 
data  that  are  invariant  with  respect  to  time.  A  stationarity  data  series 
is  characterized  by  a  constant  mean,  variance  and  covariance  throughout  the 
series,  i.e.,  no  change  over  time. 

c.  Data  Transformations.  It  is  uncommon  for  a  data  series  existing  in 
its  natural  form  to  be  stationary.  Thus,  the  data  must  be  transformed  to 
achieve  stationarity.  Three  major  transformations  were  applied  to  the  data 
to  achieve  stationarity:  (1)  differencing  the  series,  (2)  applying  natural 


log  transformations  to  the  series,  and  (3)  applying  a  square  power  trans¬ 
formation  to  the  series.  If  these  techniques  did  not  produce  stationary 
data,  then  differencing  of  the  logged  series  (#2)  or  differencing  of  the 
squared  series  (#3)  was  attempted.  Model  applications  to  differenced 
series  are  referred  to  as  integrated  and  noted  by  the  letter  "I." 

d.  Box- Jenkins  Models.  Once  stationarity  was  achieved,  the  data  were 
modeled  using  the  three  general  classes  of  Box-Jenkins  models:  AR,  MA,  or 
ARMA. 


(1)  Autoregressive  (AR)  Models.  AR  models  follow  the  general  form 

xt  ■  5+4>ixt-l  +  <*>2*t-2  +  ...  +  <£nxt-n  +  et 

where  6  is  drift,  xt  are  the  dependent  observations  of  the  series,  4>n 
are  the  regression  estimates  of  the  model,  and  «t  is  the  error  term.  The 
most  common  models  are  the  AR(1)  model 

xt  =  <5  +  ^lxt-l  +  «t 


and  the  AR(2)  model 


xt  =  5  +  4>]xt-l  +  <A2xt-2  +  et 

Autoregressive  models  differ  from  the  general  regression  equations  in  that 
there  are  no  independent  variables  to  regress  upon.  The  regression  is  per¬ 
formed  on  past  values  of  the  dependent  variable,  thus  the  term 
autoregressive. 6 


(2)  Moving  Average  (MA)  Models.  MA  models  follow  the  general  form 

xt  =  m+  *t  “  *let-l  “  02et-2  -  0n4t-n 

where  m  is  the  mean  of  the  series,  «t  are  the  past  error  terms  and  0n  are 
the  parametric  estimates  of  the  model.  The  most  common  form  of  MA  models 
are  the  MA(1)  model 


Xt  =  M  +  «t  - 


and  the  MA(2)  model 


xt  «  m  +  «t  -  *let-l  -  *2et-2 

Unlike  the  AR  models  which  are  a  linear  function  of  past  observations,  MA 
models  are  a  linear  combination  of  past  errors.  Also,  unlike  the  general 
moving  average  models  where  the  sum  of  parameters  equals  1  (0i  +  03  +  ••• 
0n  =  1),  Box-Jenkins  MA  model  parameters  do  not  necessarily  add  up  to  1. 
Finally,  the  error  terms  of  the  model  are  assumed  to  be  distributed  nor¬ 
mally  with  a  mean  of  zero  (0)  and  a  constant  variance  (  <r2).6 
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(3)  Autoregressive  -  Moving  Average  (ARMA)  Models.  ARMA  models  are 
combination  models  which  are  derived  from  the  following  equality: 

*(B)  xt  =  0(B)  «t 

where  <t>  and  0  are  the  AR  and  MA  parameters,  xt  and  <t  are  the  past  obser¬ 
vations  and  error  terms,  and  B  is  the  backshift  operator  Bxt  =  xt-l«  In 
essence,  this  equality  states  that  a  complex  AR  process  can  be  expressed  as 
a  MA  process  of  infinite  order  and  vice  versa.  The  resultant  of  this 
equality  is  the  general  equation  for  forecasting  Xt: 

xt  s  *lxt-l  +-*.+  *nxt-n  +  5  +  €t  -  *let-l  -•••-  »n‘t-n 

The  combination  of  AR  and  MA  terms  produces  a  model  that  is  more  accurate 
than  the  pure  MA  or  AR  models. 6  A  more  detailed  discussion  of  the  Box- 
Jenkins  modeling  approach  is  contained  in  Appendix  E. 

e.  To  standardize  model  identification,  all  models  are  specified  as 
autoregressive  integrated  moving  averages  (ARIMA)  of  order  p,  d,  q,  where  p 
refers  to  the  order  of  autoregressive  parameters,  d  refers  to  the  number  of 
differencing  transformations,  and  q  refers  to  the  order  of  moving  average 
parameters.  Therefore,  all  models  will  be  referred  to  as  ARIMA  (p,d,q) 
models. 

f.  Box-Jenkins  models  for  the  TWFS  were  developed  using  the  BMDP 
Statistical  Software  package.3  All  figures  depicted  in  Appendix  F  are 
copies  of  the  computer  printouts  from  BMDP  program  applications. 

6-5.  WINTERS  MODEL 

a.  Background.  Historically,  the  fitting  of  systematic  functions  to 
observations  has  typically  relied  on  least-squares  criteria  in  which  all 
the  observations  are  given  equal  weight.  However,  it  is  often  the  case, 
when  data  is  being  observed  as  a  function  of  time,  that  more  weight  should 
be  given  to  the  recent  past,  and  that  observations  taken  a  long  time  ago 
should  be  discounted  in  comparison.  In  1957,  C.  C.  Holt  published  a  paper 
entitled  "Forecasting  Seasonals  and  Trends  by  Exponential ly-weighted  Moving 
Averages."  The  procedure  proposed  therein  addressed  development  of  a  set 
of  weights  proportional  to  powers  of  a  parameter  £,  where  0  was  defined 
to  be  greater  than  zero  but  less  than  unity.  Thus,  the  set  of  weights  were 
1,  0*  02*  etc.  Constraints  were  imposed  whereby  the  sum  of  the  weights 

must  equal  unity,  and  0  must  serve  to  minimize  the  mean  square  error. 

Holt  ultimately  considered  two  parameters,  the  import  of  the  second  being 
to  account  for  a  trend  in  the  data.  In  1960,  P.R.  Winters  extended  Holt's 
method  to  cover  seasonal  effects.  Thus,  the  model  for  which  he  is  respon¬ 
sible  is  c  three-parameter  model. 
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b.  Applications.  Winters  model  owes  its  development  primarily  to  thfe 
fact  that  there  are  many  time  series  that  cannot  be  adequately  modeled  by  a 
polynomial.  Time  series  with  cyclical  or  seasonal  variations  fall  into 
this  category.  For  example,  at  least  a  cubic  equation  (which  has  a  single 
point  of  inflection  between  regions  of  upward  and  downward  concavity)  is 
required  to  capture  the  cyclical  pattern  of  periodic  data.  Furthermore, 
from  an  applications  viewpoint,  many  industrial  time  series  exhibit  sea¬ 
sonal  behavior.  Good  examples  include  the  seasonal  movement  of  specific 
commodities  such  as  POVs. 

c.  General  Form.  The  general  form  of  Winters  model  expresses  an  ob¬ 
servation  xt  at  time  t  as 

xt  *  (ai  +  b2t)  ct  +  «t* 

The  three  parameters  of  the  model  are  ax,  b2,  and  ct,  while  the  term  «t  1S 
taken  to  represent  the  usual  random  error  component.  The  parameter  ax  is 
called  the  permanent  component,  and  is  analogous  to  a  y-intercept. 
Similarly,  the  parameter  b2,  or  trend  factor,  corresponds  to  the  slope  of  a 
simple  linear  equation.  The  third  parameter,  ct»  represents  a  set  of  sea¬ 
sonal  factors  for  each  cycle.  The  seasonal  factors  induce  fluctuations 
above  and  below  the  line  segments  that  are  fitted  to  each  cycle.  The 
Winters  model  as  described  herein  is  a  multiplicative  seasonal  model,  so 
named  because  the  seasonal  parameter  ct  is  applied  multiplicatively,  not 
additively.  Multiplicative  seasonal  models  are  most  appropriate  for  time 
series  in  which  the  amplitude  (or  excursion)  of  the  seasonal  pattern  is 
proportional  to  the  average  level  of  the  series.  This  pattern  was  evident 
in  the  TWFS  data. 

d.  Specific  Form.  The  specific  form  of  Winters  prediction  equation  is 


Wr)-[Si(r)+^)r]fr+,(r+f-i) 

where,  conventionally,  carats  are  used  to  denote  estimates.  The  equation 
gives  the  forecast  at  time  T  for  an  observation  at  time  T  +  r  .  Quantities 
in  parentheses  indicate  the  times  of  computation  of  the  estimates.  Thus 
in  order  to  forecast  period  T  +t,  the  seasonal  factor  which  was  computed 
one  season  (L  periods)  ago  in  period  T  +r-i  must  be  used. 

e.  Parameters .  As  mentioned  earlier,  the  three  parameters  of  the 
Winters  model  are  the  permanent  component,  the  trend  component,  and  the 
seasonal  factor.  Estimates  of  these  parameters  for  the  period  T  are 
weighted  and  combined  with  estimates  from  previous  periods.  The  manner  in 
which  the  current  estimate  of  a  parameter  is  apportioned  with  respect  to  a 
previous  value  is  such  that  the  mean  square  error  is  minimized  over  the 
entire  time  series.  Smoothing  constants  (or  weights)  are  used  to  apportion 
present  and  past  estimates.  For  example,  if  the  smallest  mean  square  error 
were  produced  by  a  weight  of  0.80  for  the  current  estimate  of  a  parameter 
and  0.20,  for  the  previous  estimate  of  the  parameter,  then  this  would  mean 
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simply  that  the  current  estimate  is  four  times  as  important  in  the 
parameter  updating  process  as  the  previous  estimate. 

(1)  Permanent  Component.  The  estimate  of  the  permanent  component  is 
updated  by 


a|(r)~aMr-L)  +(I_a)[5i(r“  0+ W- 1)] 

where  0<  «  <  1.  Note  that  the  value  of  xy  is  divided  by  cy  (T  -  L),  which 
is  the  estimate  of  the  seasonal  factor  for  period  T  computed  one  season  (L 
periods)  ago.  This  is  done  in  order  to  eliminate  seasonal  fluctuations 
from  xt»  i.e.,  to  deseasonalize  the  current  observation.  The 
deseasonalized  observation  is  then  combined  with  the  contribution  of  the 
permanent  component  and  trend  for  the  previous  period  T-l.  This  shifts  the 
origin  of  time  to  the  end  of  the  current  period.  The  adjustment  for  sea¬ 
sonality  can  best  be  understood  by  considering  the  case  when  cy  (T-L)  is 
greater  than  1.  This  occurs  when  the  value  in  period  T-L  is  greater  than 
average  in  its  seasonality.  Dividing  xy  by  this  number  greater  than  1 

gives  a  value  that  is  smaller  than  the  original  value  by  a  percentage  just 
equal  to  the  amount  that  the  seasonality  of  period  T-L  was  higher  than  the 
average.  Of  course,  the  opposite  adjustment  occurs  when  the  seasonality 
number  is  less  than  unity.  It  should  be  noted  that  the  reason  for  using 
the  seasonal  factor  from  the  previous  season  (L  periods  ago)  is  that  the 
seasonal  factor  for  the  current  season  cannot  be  computed  until  the 
permanent  component  itself  is  calculated. 

(2)  Trend  Component.  The  estimate  of  the  trend  component  is  updated 
by 


far)  -  /8[a,(r)  -  a,(r- 1>] + (i  -  p)b2(T-  » 


where  0<  0  <  1.  This  equation  is  exactly  as  Holt's  equation  for  smoothing 
the  trend.  The  estimate  of  the  trend  component  is  simply  the  smoothed  dif¬ 
ference  between  two  successive  estimates. of  the  permanent  component.  The 
procedure  of  determining  the  trend  component  is  similar  to  evaluating  the 
slope  of  a  line  segment,  where  the  endpoints  of  the  line  segment  correspond 
to  the  beginning  and  end  of  the  period  T. 

(3)  Seasonal  Factor.  The  estimate  of  the  seasonal  factor  is  updated 
by 


*AT)mywn  +(,-*)W-z.) 


where  0  <7  <  1.  This  equation  specifies  the  seasonal  index  as  the  ratio 
of  the  current  value  of  the  series,  xy,  to  the  current  smoothed  value  for 
the  series,  <ri(T).  If  xy  is  larger  than  ai(T),  the  ratio  will  be  greater 
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than  1,  while  if  it  is  smaller  than  ai(T),  the  ratio  will  be  less  than  1. 

It  is  important  to  understand  that  al(T)  is  a  smoothed  average  value  of  the 
series  that  does  not  include  seasonality.  The  values  of  xj,  however,  do 
contain  seasonality  (as  well  as  randomness).  Notice  that  equation  smoothes 
(weights)  the  current  observed  seasonal  variation  (xj/ai(T))  with  the 
estimate  of  the  seasonal  factor  for  period  T  computed  L  periods  ago.  That 
was  the  last  opportunity  to  observe  this  portion  of  the  seasonal  pattern. 

f.  Smoothing  Constants  and  Model  Initialization.  The  method  of  ini¬ 
tializing  model  parameters  and  of  solving  for  the  smoothing  constants  is 
not  central  to  understanding  the  Winters  method.  A  short  explanation  of 
these  procedures  has  been  included  as  Appendix  G. 

6-6.  POV  FORECASTING  MODEL  DEVELOPMENT 

a.  POV  data  is  fairly  constant  on  an  annual  basis;  however,  the  monthly 
data  is  highly  seasonal.  Stationary  data  was  achieved  after  differencing 
the  series  at  a  lag  of  12  months. 

b.  The  appropriate  model  for  POV  data  was  a  nonseasonal  ARIMA  (2,0,0) 
model  and  seasonal  ARIMA  (0,1,1)  model: 


(1  -  .02216  -  .3554  B2)  (1  -  Bl2)  Xt  =  (1  -  .8581  b12) 

(.19)  (3.05)  (16.77) 

X2  (3,20)  =  16.35 

The  results  indicate  that  all  of  the  estimated  coefficients  of  the  model, 
except  are  significant  at  the  5  percent  significance  level.  Verifica¬ 
tion  of  the  model's  adequacy  using  the  Box-Pierce  test  indicates  that  the 
autocorrelations  of  the  estimated  model  residuals  are  not  significantly 
different  from  zero  at  the  30  percent  significance  level.  (Details  of  the 
Box-Pierce  test  and  model  adequacy  are  contained  in  paragraph  E-6,  Appendix 
E).  Using  this  model  to  forecast  FY  82  POV  shipments  resulted  in  a 
forecast  of  501.37  K/MTON  versus  an  actual  shipment  of  510.25  K/MTON,  which 
corresponds  to  an  underforecast  of  1.7  percent  (see  Figure  6-2).  The  MSE 
of  the  POV  forecast  for  FY  82  was  8.54. 

c.  Variations  of  the  aggregate  POV  data  model  were  used  to  forecast  POV 
mode  requirements.  The  POV  container  shipments  during  FY  82  were  under¬ 
forecast  by  20.1  percent.  The  POV  breakbulk  shipment  forecast  was  over¬ 
forecast  by  8.8  percent.  Specific  details  pertaining  to  Box-Jenkins 
modeling  of  POV  data  and  all  other  commodities  are  contained  in  Appendix  F. 

d.  When  the  Winters  model  was  used  to  forecast  POVs  for  1982,  an  under¬ 
forecast  of  5.3  percent  occurred.  The  aggregate  forecast  was  483.12 
K/MTON,  corresponding  to  an  MSE  of  12.92.  Separation  of  this  commodity  by 
mode  of  shipment  produced  smaller  errors  of  -3.3  percent  (breakbulk)  and 
3.0  percent  (containerized).  However,  based  on  increased  mean  square 
errors,  the 


6-10 


overall  fits  were  not  as  good.  The  specific  Winters  parameters  for  each 
commodity  and  month  have  been  included  at  Appendix  H. 


Figure  6-2.  POV  Cargo  Forecast  -  FY  82 


6-7.  GENERAL  CARGO  FORECASTING  MOOEL  DEVELOPMENT 

a.  General  cargo  accounts  for  the  largest  amount  of  cargo  shipped  under 
any  single  commodity  code.  General  cargo  accounts  for  approximately  50 
percent  of  all  cargo  shipped  by  the  US  Army  during  any  individual  year. 

This  commodity  code  not  only  includes  operational  and  facility  maintenance 
items,  but  also  force  modernization  program  data.  The  inclusion  of  this 
force  modernization  data  within  the  time  series  data  hinders  the  develop¬ 
ment  of  accurate  forecast  models  for  the  general  cargo  commodity.  If  new 
TAC  codes  are  used  to  isolate  force  modernization  equipment,  the  accuracy 
of  general  cargo  forecasts  will  improve.  To  demonstrate  this  point,  the 
general  cargo  data  series  was  separated  into  two  components:  (1)  A205  TAC 
data— primarily  DARCOM  items  and  (2)  non-A205  TAC  data— other  items  shipped 
as  general  cargo.  The  results  of  this  effort  will  be  discussed  later. 


b.  Annual  general  cargo  data  fluctuates  within  a  narrow  range  of  ±10 
percent;  however,  the  data  is  seasonal.  Stationarity  was  achieved  after 
two  successive  differencing  operations.  First  the  series  was  differenced 
at  a  lag  of  12  months.  This  operation  did  not  produce  stationary  data, 
thus  the  seasonally  adjusted  series  was  differenced  at  a  lag  of  1  month. 

c.  The  appropriate  model  for  general  cargo  data  was  a  mixed  ARIMA 
(1,1,1)  nonseasonal  model  and  a  seasonal  ARIMA  (0,1,1)  model: 

(1  -  .35(79  B)  (1  -  .12,  (1  -  •)  *t  ■  (1  -  U  -  .ffijf*)  ‘t 

x2  (3,20)  =  12.8 


The  results  indicate  that  all  of  the  estimated  coefficients  are  signficant 
of  the  5  percent  significance  level.  The  Q  statistic  for  20  degrees  of 
freedom  is  12.8.  Thus,  the  null  hypothesis,  estimated  model  residuals  are 
uncorrelated,  should  not  be  rejected  at  the  10  percent  significance  level. 
The  general  cargo  seasonal  model  forecasted  that  1856.61  K/MTON  would  be 
transported  during  FY  82  versus  an  actual  shipment  of  2021.55  K/MTON, 
equating  to  an  underforecast  of  4.6  percent  error  (see  Figure  6-3).  (Note: 
another  model  developed  for  general  cargo  underforecast  the  FY  82  lift  by 
4.2  percent.  The  goodness  of  fit  of  this  model  was  not  as  good  as  the 
ARIMA  (1,1,1)  x  (0,1,1)  model,  but  the  forecast  accuracy  of  both  models 
should  be  tracked  in  the  future.)  The  MSE  of  the  general  cargo  forecast 
for  FY  82  was  318.54. 

d.  Box-Jenkins  modeling  of  the  general  cargo  mode  shipments  produced  a 
9.8  percent  underforecast  of  container  data  and  an  8.5  percent  underfore¬ 
cast  of  breakbulk  data.  The  non-A205  (non-DARCOM)  forecast  for  general 
cargo  was  underforecast  by  2.2  percent. 

e.  Application  of  the  Winters  model  to  general  cargo  resulted  in  a 
forecast  of  1870.38  K/MTON.  The  prediction  was  about  7.5  percent  below  the 
actual  shipment,  while  the  MSE  was  279.89. 

f.  Several  excursions  were  conducted  on  general  cargo  in  order  to 
assess  this  important  category  accurately.  The  Winters  model  was  accurate 
to  within  1.4  percent  for  non-A205  (non-DARCOM)  shipments  (breakbulk  and 
container  combimed)  and  to  within  1.6  percent  for  non-A205  (container 
only).  Overall,  the  Winters  model  did  not  fit  the  breakbulk  time  series 
well  (28.6  percent  overforecast) .  The  Winters  general  container  forecast 
was  within  4.9  percent  of  the  actual  data. 
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Figure  6-3.  General  Cargo  Forecast  -  FY  82 


6-8.  HHG  FORECASTING  NOOEL  DEVELOPMENT 

a.  Transportation  arrangements  for  DOD-sponsored  household  goods  (HHG) 
are  made  by  MTMC,  using  primarily  two  methods:  freight  forward  (Code  4), 
and  those  that  are  MSC  processed  (Code  5).  Freight  forward  shipments  are 
shipments  of  HHG  that  are  contracted  with  private  carriers  from  door  to 
door;  shipments  of  HHG  that  are  shipped  using  the  assets  of  MSC  are  Code  5 
shipments.  Additionally,  some  HHG  are  returned  to  CONUS  via  MAC  when  cargo 
space  is  available. 

b.  Data  used  in  the  study  only  consisted  of  Code  5  shipments.  Accord¬ 
ing  to  other  HHG  shipment  data  that  was  gathered  from  MTMC,  Code  5  ship¬ 
ments  comprise  less  than  4  percent  of  the  total  HHG  tonnaged  shipped  during 
a  given  fiscal  year.  However,  any  aberration  in  the  commercial  shipping 
process  will  cause  HHG  shipments  by  MSC  to  increase  significantly.  For 
example,  assuming  that  Code  5  shipments  comprise  4  percent  of  total  HHG 
transported,  a  1  percent  change  in  Code  4  shipments  due  to  a  commercial 


shipping  strike  will  increase  Code  5  shipments  by  25  percent.  This  fact 
complicates  the  development  of  accurate  models  to  forecast  the  amount  of 
HHG  shipped  by  MSC.  It  also  lends  credence  to  the  argument  that  MTMC 
should  be  involved  in  the  HHG  forecasting  process.  An  accurate  forecasting 
model  for  all  HHG  shipments  cannot  be  developed  until  additional 
probability  and  economic  analysis  of  commericial  shipping  is  integrated 
into  the  forecast  model. 

c.  Only  data  from  October  1978  through  September  1981  was  used  to  de¬ 
velop  the  forecasting  models  due  to  the  volatility  of  HHG  during  the  period 
under  study.  As  noted  in  Figure  F-26,  Appendix  F,  HHG  shipments  in  FY  78 
were  approximately  three  times  greater  than  any  other  year  due  to  commer¬ 
cial  shipping  aberrations.  Annual  HHG  data  during  FY  79  through  FY  81  ex¬ 
hibited  a  slight  trend  downward,  but  the  monthly  data  was  seasonal.  To 
achieve  stationarity,  the  data  were  differenced  at  a  lag  of  12  months. 

d.  The  appropriate  model  for  HHG  data  was  a  nonseasonal  ARIMA  (0,0,1) 
model  and  a  seasonal  ARIMA  (1,1,0)  model  with  a  series  mean: 


(1  +  .3198812)  (1-B12)  Xt  «  1.163  +(1-.3397B^) 
(-2.33)  „  (-3.67)  (1.88) 

x2  (2,20)  =  6.26 


All  of  the  estimated  coefficients,  except  04,  are  significant  at  the  5 
percent  significance  level.  04  is  significant  at  the  7  percent  signifi¬ 
cance  level.  Verification  of  the  model's  adequacy  using  the  Box-Pierce 
test  with  20  degrees  of  freedom  indicates  that  the  null  hypothesis  (uncor¬ 
related  residuals)  should  not  be  rejected  at  the  .5  percent  significance 
level.  This  model  forecasted  FY  82  HHG  shipments  to  be  80.19  K/MTON  versus 
an  actual  lift  of  87.68  K/MTON  (see  Figure  6-4).  The  MSE  of  the  forecast 
was  5.24  and  the  annual  forecast  error  was  -8.5  percent.  Mode  forecast  of 
HHG  data  resulted  in  a  23.4  percent  overforecast  of  container  data  and  a 
16.9  percent  underforecast  of  breakbulk  data. 

e.  The  Winters  fit  resulted  in  a  forecast  of  76.079  K/MTON,  amounting 
to  an  overall  underforecast  of  13.2  percent.  The  MSE  was  3.361. 

Separation  of  household  goods  into  breakbulk  and  containerized  shipments 
only  served  to  worsen  the  forecast.  The  breakbulk  fit  produced  an 
overforecast  of  43  percent,  while  the  containerized  model  resulted  in  a 
36.5  percent  underforecast. 
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Figure  6-4.  HHG  Cargo  Forecast  Data  -  FY  82 


6-9.  COAL  FORECASTING  MODEL  DEVELOPMENT 

a.  Coal  shipments  are  concentrated  on  the  shipping  routes  between  East 
Coast/Gulf  Coast  to  Europe.  Less  than  .01  percent  of  all  coal  shipped  is 
inbound  for  ports  other  than  Europe.  Coal  forecasts  are  prepared  by  the 
Energy  Center,  Rheinau  and  forwarded  to  Defense  Logistics  Agency  (DLA- 
DFSC)«  DLA  is  responsible  for  soliciting  the  contracts  for  the  coal  re¬ 
quirements  and  approving  the  negotiations.  Once  contracts  are  approved, 
the  Energy  Center  at  Rheinau  and  the  MSC  coordinate  the  shipment  scheduling 
of  coal. 

b.  During  October  1977  to  September  1981  (timeframe  used  to  construct 
the  statistical  forecasting  model),  several  factors  affecting  coal  ship¬ 
ments  occurred. 

(1)  In  recent  years,  contract  negotiations  for  coal  procurements  have 
been  delayed  due  to  several  reasons,  which  in  turn  have  delayed  coal  ship¬ 
ments  during  a  given  fiscal  year.  Thus,  the  data  obscure  the  fact  that 
some  coal  shipments  which  are  ordered  to  be  delivered  in  one  fiscal  year 
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are  delayed  and  actually  shipped  in  another  fiscal  year.  Over  time,  this 
accordion  effect  will  disappear.  However,  accurate  forecasting  for  any 
given  year  is  hindered  because  the  current  data  are  clouded  by  events  ex¬ 
ternal  to  the  forecasting  system. 

(2)  Ouring  the  conservation  efforts  of  the  past  few  years,  coal 
burners  in  Europe  have  transferred  their  source  of  fuel  from  anthracite  to 
bituminous.  Although  the  effort  is  made  for  cost  and  conservation  pur¬ 
poses,  the  true  impact  of  the  change  cannot  be  determined  from  the  data. 

The  total  tonnage  of  coal  shipped  has  not  changed  appreciably  since  any 
reduction  of  anthracite  has  been  mostly  offset  by  the  increase  in 
bituminous. 

In  sum,  mathematical  forecasting  for  coal  shipments  should  be  modified 
based  upon  knowledge  of  the  energy  program  and  external  events  that  impact 
on  the  program. 

c.  Coal  data  during  the  period  under  study  did  not  exhibit  a  recurring 
pattern  or  seasonal  trend.  However,  it  is  evident  that  some  adjustments  to 
the  raw  data  are  necessary  to  improve  the  accuracy  of  the  forecasts. 
Stationary  data  was  achieved  after  differencing  the  raw  series  at  a  lag  of 
1  month. 

d.  The  appropriate  model  for  coal  data  was  an  ARIMA  (1,1,1)  model: 


(1  -  .34618)  (1  -  8)  Xt  =  (1  -  .859B)  et 
(9.46)  (2.06) 

x2  (2,20)  =  10.87 


The  results  indicate  that  all  of  the  estimated  coefficients  are  significant 
at  the  5  percent  significance  level.  The  estimated  ARIMA  (1,1,1)  model  was 
verified  using  the  Box-Pierce  test  and  a  visual  examination  of  the 
estimated  residuals  of  the  model.  The  Q  statistic  test  is  10.87,  which 
indicates  that  the  null  hypothesis,  uncorrelated  residuals,  should  not  be 
rejected  at  the  6  percent  significance  level.  The  ARIMA  (1,1,1)  model 
forecasted  the  FY  82  coal  tonnage  shipped  to  be  373.53  K/MTON  as  compared 
to  the  historical  FY  82  lift  of  375.94  K/MTON  (-.6  percent)  and  the  actual 
FY  82  DLA  coal  contract  of  440.67  K/MTON  (15.2  percent  error).  The 
forecast  MSE  was  192.60.  Figure  6-5  illustrates  the  monthly  forecast. 

e.  When  the  Winters  model  was  used  to  predict  the  historical  FY  82 
lift,  a  forecast  of  351.93  K/MTON  resulted.  Although  the  average  forecast 
error  was  only  6.4  percent,  the  MSE  was  fairly  high  (210.18),  indicating  a 
relatively  poor  fit. 
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Figure  6-5.  Coal  Cargo  Forecast  -  FY  82 


6-10.  AMMUNITION  FORECASTING  MODEL  DEVELOPMENT 

a.  Ammunition  requests  fall  into  two  general  categories:  (1)  ammu¬ 
nition  to  increase  war  reserve  stocks  and  (2)  ammunition  to  be  used  for 
training  purposes.  The  war  reserve  buildup  program  is  an  annually  funded 
program  that  is  reviewed  every  6  months  to  address  the  manufacturing, 
transportation,  handling,  and  stockpiling  constraints  of  ammunition 
buildup.  Ammunition  dedicated  for  training  purposes  is  strictly  based  upon 
demand  and  number  of  combat  units  in  a  particular  theater.  Additionally, 
safety  testing  of  ammunition  complicates  the  accurate  forecasting  of 
ammunition  requirements.  If  manufactured  ammunition  fails  safety  tests, 
then  ammunition  shipment  plans  must  be  rearranged. 
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b.  The  ammunition  movement  planning  process  consists  of  two  principal 
agencies:  OCSLOG  and  Armament,  Munitions  and  Chemical  Command  (AMCCOM). 
DCSLOG  functions  as  a  screening  filter  for  the  war  reserve  buildup  program 
and  training  ammunition  requests  from  theater  commands.  Requests  are 
reviewed  bi annually  at  Committee  on  Ammunition  Logistic  Support  (CALS) 
meetings,  where  ammunition  allocations  are  determined  for  each  command 
based  upon  priority  and  manufacturing  constraints.  The  approved 
allocations  are  then  used  by  AMCCOM  to  direct  movement  of  ammunition  to  the 
selected  commands.  At  the  subsequent  CALS  meeting  (every  6  months),  the 
status  of  ammunition  reserves  and  manufacturing  capability  is  reviewed  and 
adjustments  are  made  to  achieve  the  stated  goals.  As  with  coal,  annual 
ammunition  forecasts  should  be  based  upon  mathematical  models  and  inside 
information  regarding  the  status  of  ammunition  programs. 

c.  Ammunition  data  is  relatively  constant  on  an  annual  basis  and  the 
monthly  data  exhibits  some  seasonality.  However,  attempts  to  capture  the 
seasonal  aspects  of  the  data  were  futile  since  stationarity  could  not  be 
achieved.  Like  coal,  the  ammunition  data  should  be  monitored  carefully  and 
adjusted  with  caution.  Stationary  data  was  achieved  after  differencing  the 
original  series  at  a  lag  of  1  month. 

d.  The  appropriate  model  for  ammunition  data  was  an  ARIMA  (2,1,4) 
model : 


(1+.8418B  +  .5399  B2)  (1-B)  Xt  ■  (1+.3454B4  -  .368985  -  .3875B6  -  .2561B?) 
(-7.33)  (-4.59)  (-3.11)  (3.86)  (4.57)  (2.43) 

x2  (6,20)  =  17.61 


All  of  the  estimated  coefficients  are  significant  at  the  5  percent  signifi¬ 
cance  level.  The  Box-Pierce  test  for  the  model  indicates  that  the  null 
hypothesis,  estimated  residuals  are  uncorrelated,  should  not  be  rejected  at 
the  40  percent  significance  level.  The  model  forecasted  that  ammunition 
shipments  during  FY  82  would  be  189.62  K/MTON  versus  an  actual  lift  of 
179.31  K/MTON  of  ammunition  during  FY  82  (see  Figure  6-6).  The  MSE  of  the 
ammunition  forecast  was  119.91  and  the  annual  error  was  5.7  percent. 

e.  When  the  Winters  model  was  used  to  forecast  ammunition,  a  forecast 
of  158.080  K/MTON  resulted,  which  corresponds  to  an  error  of  -11.8  percent. 
The  fit  was  quite  good,  except  for  the  existence  of  a  heavy  ammunition 
shipment  (outlier)  during  July  of  FY  82.  MSE  was  107.59. 
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Figure  6-6.  Ammunition  Cargo  Forecast  -  FY  82 


6-11.  SPECIAL  FORECASTING  MOOEL  DEVELOPMENT 

a.  The  majority  of  items  shipped  under  commodity  code  special  comprise 
program  items  that  relate  to  the  Army  Modernization  Improvement  Memorandum 
(AMIM).  They  include  items  such  as  Ml  tanks,  5-ton  trucks,  and  track  ve¬ 
hicles.  However,  there  are  some  program  items  which  are  small  enough  to  be 
shipped  in  the  general  cargo  category  such  as  force  modernization  program 
equipment  shipped  by  Communications  and  Electronics  Command  (CECOM). 

Radios  and  other  electronic  items  that  are  shipped  by  CECOM  that  are  small 
in  size  are  shipped  as  general  cargo.  Items  shipped  under  commodity  code 
special  are  primarily  AMIM  program  items;  however,  not  all  AMIM  items  are 
shipped  as  special  cargo.  For  purposes  of  this  study,  AMIM  items  could  not 
be  isolated  from  the  special  cargo  category  due  to  the  current  TAC 
structure.  However,  with  improved  data  base  management  and  expanded  TAC, 
more  accurate  forecasts  of  the  special  cargo  commodity  can  be  made. 


6-1 


b.  Special  cargo  exhibited  a  trend  that  could  possibly  be  interpreted 
as  program  funding  of  the  defense  budget  and  some  aspects  of  seasonality. 

To  achieve  stationary  data,  the  original  series  was  differenced  at  a  lag  of 
12  months. 

c.  The  appropriate  model  for  special  cargo  was  an  ARIMA  (2,0,0)  non- 
seasonal  model  and  an  ARIMA  (0,1,1)  model: 


(1  -  .1496B  -  .3875B2)  (1  -  b!2)  lnXt  =  (1  -  .801B12  .  .150314)  <t 
(1.21)  (3.16)  (10.61)  (1.65) 

x2  (4,20)  =  13.15 


All  of  the  estimated  coefficients  except  4>i  and  are  significant  at  the 
5  percent  significance  level.  The  model  was  verified  using  the  Box-Pierce 
test  and  a  visual  check  of  the  estimated  residuals.  The  Q  statistic  for 
this  model  with  20  degrees  of  freedom  is  13.15.  This  diagnostic  indicates 
that  the  residuals  are  not  significant  and  that  the  null  hypothesis  should 
not  be  rejected  at  the  16  percent  significance  level.  This  model  fore¬ 
casted  that  559.42  K/MT0N  of  special  cargo  would  be  shipped  during  FY  82 
versus  an  actual  amount  of  189.83  K/MT0N  (see  Figure  6-7).  The  annual 
forecast  error  was  14.2  percent  and  the  forecast  MSE  was  159.23.  It  should 
be  noted  that  the  amount  of  special  cargo  tonnage  shipped  during  FY  82  is 
much  lower  than  any  of  the  preceding  years.  In  fact,  the  average  tonnage 
of  special  cargo  shipped  over  the  past  5  years  was  669.17  K/MT0N.  In  all 
cases,  forecasting  methodologies  are  employed  based  upon  the  assumption 
that  the  behavior  of  future  observations  will  not  differ  greatly  from  the 
behavior  of  past  observations.  When  the  future  is  radically  different  than 
the  past,  as  in  the  case  of  special  cargo  tonnage  for  FY  82,  forecasts 
cannot  predict  the  future  with  any  degree  of  accuracy. 

d.  Similarly,  the  radical  change  in  special  cargo  shipments  during  FY 
82  caused  the  Winters  model  to  overestimate  the  actual  lift  by  28.1  per¬ 
cent.  The  forecast  MSE  was  404.018;  the  Winters  model  forecasted  a  lift  of 
627.52  K/MT0N. 
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Figure  6-7.  Special  Cargo  Forecast  -  FY  82 


6-12.  CARGO  TRAILER/CONEX  FORECASTING  MODEL  DEVELOPMENT 

a.  CONEX  shipments,  as  annotated  in  the  data  (commodity  70),  were  in 
fact  said  to  be  roll  on/roll  off  (RORO)  trailer  shipments.  Empty  CONEX 
data  was  negligible,  therefore,  it  was  combined  with  trailer  data  and  the 
resulting  series  was  modeled  as  cargo  trailers.  CONEX  data  demonstrated  an 
upward  annual  linear  trend  and  the  monthly  data  exhibited  seasonality. 
Stationarity  was  achieved  after  differencing  the  original  series  at  a  lag 
of  12  months. 

b.  The  appropriate  model  for  special  cargo  data  was  an  ARIMA  (0,0,1) 
nonseasonal  model  and  an  ARIMA  (0,1,1)  seasonal  model  with  a  series  mean: 


(1  .  b12)  Xt  =  .4057  +  (1  +  .2912B)  (1 
0  (4.09)  (-2.32) 

x2  (2,20)  =  12.74 


.8264B12)  et 
(14.97) 
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The  results  indicate  that  all  of  the  estimated  coefficients  are  significant 
at  the  5  percent  significance  level.  The  Q  statistic  for  the  estimated 
model  residuals  at  20  degrees  of  freedom  is  12.74  which  is  significant  at 
the  16  percent  significance  level.  Therefore,  the  null  hypothesis,  uncor¬ 
related  residuals,  should  not  be  rejected  at  the  16  percent  significance 
level.  The  model  forecasted  CONEX/cargo  trailer  shipments  for  FY  82  to  be 
76.43  K/MTON  versus  an  actual  lift  for  FY  82  of  74.47  K/MTON  (see  Figure  6- 
8).  The  MSE  of  the  CONEX  forecast  was  .81,  and  the  annual  forecast  error 
was  2.6  percent. 

c.  When  the  Winters  model  was  used  to  estimate  CONEX/cargo  trailer 
shipments,  the  forecast,  77.03  K/MTON,  exceeded  the  actual  lift  by  2.56 
K/MTON,  or  about  3.4  percent.  The  MSE  was  .78  K/MTON,  and  the  fit  was 
quite  good. 
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Figure  6-8.  CONEX  Cargo  Forecast  -  FY  82 
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6-13.  FREEZE  FORECASTING  MODEL  DEVELOPMENT 

a.  AAFES  shipments  of  freeze  goods  account  for  approximately  80  percent 
of  freeze  cargo  shipped  during  a  given  year.  Freeze  shipments  increased  24 
percent  during  FY  82;  however,  AAFES  freeze  shipments  increased  30  percent 
during  the  same  period.  Thus,  like  coal  and  ammunition,  program 
information  from  AAFES  must  be  integrated  into  the  overall  forecast  of 
freeze  data. 

b.  Freeze  data  depicts  a  trend  and  pattern  very  similar  to  CONEX  data. 
The  trend  during  the  period  under  study  was  an  upward  linear  function,  and 
the  monthly  data  exhibited  seasonality.  Stationarity  was  achieved  after 
differencing  the  seasonally  adjusted  series  at  a  lag  of  1  month. 

c.  The  appropriate  model  for  freeze  data  was  an  ARIMA  (2,1,2)  non- 
seasonal  model  and  a  (0,1,1)  seasonal  model: 


(1+.7740B  +  .6789B2)  (1-812)  (1-B)  Xt  =  (1-.4729B  +  .5038BU)  (1-.7753B12)  et 
(-6.76)  (-5.96)  (4.32)  (-5.76)  (12.04) 

x2  (5,20)  =  14.66 

The  results  indicate  that  all  of  the  coefficients  of  the  model  are  signifi¬ 
cant  at  the  5  percent  significance  level.  The  Q  statistic  for  20  degrees 
of  freedom  is  14.66.  Therefore,  the  null  hypothesis,  residuals  are  uncor¬ 
related,  should  not  be  rejected  at  the  20  percent  significance  level. 

The  freeze  seasonal  model  forecasted  that  31.86  K/MTON  would  be  transported 
during  FY  82  versus  an  actual  shipment  of  36.42  K/MTON  (see  Figure  6-9). 

The  MSE  of  the  freeze  forecast  was  .49,  and  the  annual  forecast  error  was  - 
12.5  percent. 

d.  When  the  Winters  model  was  used  to  forecast  freeze  shipments,  a  12.8 
percent  underforecast  resulted.  The  model  forecasted  a  lift  of  31.758 
K/MTON,  and  the  MSE  of  the  forecast  was  0.458,  indicative  of  a  relatively 
good  fit. 
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Figure  6-9.  Freeze  Cargo  Forecast  -  FY  82 


6-14.  CHILL  FORECASTING  MODEL  DEVELOPMENT 

a.  As  stated  previously,  raw  billing  data  from  MSC  were  compiled  into 
monthly  time  series  data  (cargo  shipped  monthly).  The  sorting  process  used 
to  obtain  the  monthly  time  series  data  was  done  under  the  assumption  that 
debits  and  credits  contained  In  the  raw  data  base  would  be  eliminated  by 
this  process.  Furthermore,  it  was  assumed  that  this  data  manipulation 
would  not  produce  any  negative  numbers  (i.e.,  credit  balance).  This  as¬ 
sumption  follows  from  normal  accounting  practices  when  a  credit  is  recorded 
to  offset  previous  debit.  However,  due  to  problems  with  the  raw  data  base, 
a  negative  number  was  generated  from  this  process  (November  1976— see 
Figure  D-l,  Appendix  D).  The  cause  of  this  negative  number  can  only  be 
resolved  after  a  thorough  review  of  the  MSC  and  Army  finance  and  accounting 
records  which  is  beyond  the  scope  of  this  study.  Therefore,  the  first  2 
months  of  the  data  series  were  eliminated.  Time  series  data  from  December 
1977  through  September  1981  were  used  to  develop  forecasting  models  for  the 
corrmodlty  chill. 

b.  Chill  shipments  were  relatively  constant  on  an  annual  basis,  but 
monthly  data  did  exhibit  seasonality.  Stationarity  was  achieved  after 
transforming  the  original  series  using  natural  logarithms  and  then  dif¬ 
ferencing  the  logarithmic  series  at  a  lag  of  12  months. 


c.  The  appropriate  model  for  chill  data  was  a  nonseasonal  ARIMA  (0,0,2) 
model  and  a  seasonal  ARIMA  (1,1,0)  model: 


(1  +  .6371812) (l  .  b12)  ]n  xt  -  (1  -  .3012B  +  .2716B2)  «t 
(-5.61)  (2.0)  (-1.83) 

x2  (3,20)  =  15.63 


All  of  the  estimated  coefficients  except  #2  are  significant  at  the  5  per¬ 
cent  significance  level.  The  Q  statistic  for  20  degrees  of  freedom  is 
15.63,  which  indicates  that  the  residuals  are  not  significant  and  that  the 
null  hypothesis  should  not  be  rejected  at  the  25  percent  significance 
level.  Also,  all  of  the  estimated  residuals  fall  within  the  95  percent 
confidence  interval  and  appear  to  be  distributed  randomly.  The  FY  82  fore¬ 
cast  for  chill  using  this  model  was  15.32  K/MTON  of  cargo  shipped  versus  an 
actual  lift  of  14.13  K/MTON  (see  Figure  6-10).  The  MSE  of  the  chill  fore¬ 
cast  was  .21,  and  the  annual  forecast  error  was  8.4  percent. 

d.  Application  of  Winters  model  to  the  chill  time  series  resulted  in  a 
forecast  of  13.03  K/MTON.  This  was  an  aggregate  underforecast  of  1.094 
K/MTON,  or  -7.7  percent.  The  fit  was  good,  and  the  MSE  of  the  forecast  was 
.13. 


Figure  6-10.  Chill  Cargo  Forecast  -  FY  82 


6-15.  SUMMARY  OF  FORECASTING  RESULTS.  Table  6-2  portrays  the  forecasting 
results  of  the  Box-Jenkins  and  Winters  models  for  each  commodity  for  FY  82. 


Table  6-2.  FY  82  Forecasting  Results 


S' 

Commodity 

Percent  error 

Present  system 

Box-Jenkins 

Winters 

FY  81 

FY  82 

FY  82 

FY  82 

General 

22.0 

3.5 

4.6 

7.5 

Breakbulk  (4.7) 

168.0 

63.0 

9.8 

28.6 

Container  (44.5) 

7.8 

10.0 

8.5 

4.9 

Nonprogram  (33.5) 

— 

— 

2.2 

1.4 

Special  (16.6) 

15.0 

66.0 

14.0 

28.0 

Coal  (10.4) 

22.0 

13.8 

.6 

6.4 

POV  (13.5) 

.4 

7.3 

1.7 

5.3 

Ammo  (3.9) 

36.0 

3.9 

.3 

11.8 

Household  (2.5) 

19.0 

11.0 

8.5 

13.2 

Cargo  tlr  (2.1) 

30.0 

14.6 

1.6 

2.9 

Reefer  (1.2) 

14.0 

39.2 

7.0 

10.0 

(  )  Percent  of  total  lift 


CHAPTER  7 

SATISFACTION  OF  ELEMENTS  OF  ANALYSIS 


7-1.  GENERAL.  The  course  of  the  study  was,  to  a  significant  extent, 
guided  by  all  the  sponsor's  designated  essential  elements  of  analysis 
(EEA).  However,  identification  of  the  variances  between  forecasts  and 
lift,  conditions  which  contributed  to  the  variances,  and  effective  methods 
for  improving  the  forecast  predominated  the  study  efforts.  All  of  the 
questions  in  the  EEA  were  answered  during  the  study,  albeit  some  to  a 
greater  degree  than  others. 

7-2.  SYNOPSIS  OF  EEA.  Discussion  of  the  EEA  is  contained  throughout  the 
report.  A  capsule  of  the  EEA  and  the  answers  to  them  follow: 

a.  What  are  the  recorded  variances  between  long-range  transportation 
workload  forecasts  and  actual  utilization  of  shipping  cargo?  Tables  2-1 
and  2-2,  Chapter  2,  contain  the  variances  by  commodity  and  shipping  mode 
for  FY  81  and  FY  82. 

b.  Do  systemic  conditions  exist  which  contribute  to  unrealistic  fore¬ 
casting?  If  so,  what  are  they?  Paragraph  2-2,  Chapter  2,  details  the  con¬ 
ditions  observed  in  the  current  system.  Most  significant  systemic  contri¬ 
butions  to  the  errors  are  the  lack  of  an  effective  feedback  system,  some 
agency  reports  based  on  a  limited  knowledge  of  the  forecasted  commodity, 
inadequacies  in  the  TAC  structure,  inadequate  program  input  into  the  fore¬ 
cast,  and  separate  long-  and  short-range  forecasting  systems. 

c.  What  is  the  economic  and  operational  impact  of  long-range  forecasts 
which  are  at  variance  with  actual  utilization?  Table  2-5,  Chapter  2,  shows 
the  MSC  industrial  fund  FY  82  losses  which  were  directly  attributable  to 
errors  in  Army  forecasts.  Operationally,  the  forecasting  variance  resulted 
in  underutilization  of  MSC  cargo  carrying  capacity.  The  MSC-controlled 
fleet  was  approximately  40  percent  utilized  in  FY  1982. 

d.  Do  short-range  forecasts  impact  on  long-range  forecasting?  If  so 

how?  With  the  exception  of  CONUS  outbound  household  goods  and  privately 
owned  vehicles,  long-range  and  short-range  forecasts  are  prepared  by  the 
same  personnel.  In  some  cases,  the  long-range  forecast  impacts  the  short- 
range,  as  opposed  to  vice  versa,  as  in  those  cases,  the  short-range  fore¬ 
casts  are  derived  from  the  long-range. 

e.  What  is  the  impact  of  the  current  separation  of  responsibility  for 
long-range  and  short-range  forecasting?  Separation  of  the  responsibility 
for  long-range  and  short-range  forecasting  has  created  a  dual  system  of 
reporting  which  does  not  contribute  to  the  accuracy  of  the  forecasts.  It 
has  produced  long-range  and  short-range  forecasts  which  are  usually  not  in 
agreement.  It  requires  the  maintenance  of  similar  records  at  LCA  and  HQDA, 
and  diffuses  responsibility  f6r  management  of  the  forecasting  system. 
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f.  What  Methodologies  exist  in  other  services  which  could  be  applied  to 
the  resolution  of  the  Army's  problem?  Discussion  of  USAF  and  USN  forecast¬ 
ing  procedures  is  contained  in  Chapter  3.  The  USAF  forecasting  procedures 
at  AFLC  could  be  used  in  conjunction  with  the  procedures  and  models  rec¬ 
ommended  in  this  report  as  the  basis  for  an  improved  Army  forecasting  sys¬ 
tem.  AFR  75-15,  Forecast  of  Air  Force  Transportation  Requirements,  and 
AFLC  Regulation  171-125,  Surface  Transportation  Tonnage  and  Cost  System 
(00278),  contain  a  description  of  the  USAF  system. 

g.  What  DARCOM  activities  affect  major  end  items  of  equipment  planned 
for  overseas  distribution?  All  program  and  project  management  decisions 
affecting  the  distribution  of  major  end  items  to  overseas  locations  impact 
the  requirements  for  surface  cargo  space.  In  particular,  changes  to  pro¬ 
grams  which  occur  subsequent  to  the  formulation  of  the  annual  forecast  of 
space  requirements  are  critical  in  that  the  shipping  vessel  procurement  and 
use  is  developed  from  this  forecast.  If  the  transportation  forecasting 
activity  is  not  cognizant  of  the  changes  to  programs,  then  the  match  of 
shipping  resources  to  available  cargo  is  complicated.  A  major  forecasting 
difficulty  has  been  the  information  interface  between  the  program  decisions 
and  the  forecasting  process.  The  analysis  of  the  impact  of  DARCOM  activi¬ 
ties  was  limited  by  the  masking  of  much  of  the  DARCOM  generated  cargo  with¬ 
in  the  general  and  special  cargo  commodities. 

h.  What  is  the  impact  of  the  TAEDP  on  forecasting  of  requirements?  The 
TAEDP  is  provided  to  DARCOM  Inventory  Control  Points  (ICP)  prior  to  the 
preparation  of  the  preliminary  annual  forecast.  It  is  used  in  conjunction 
with  the  Materiel  Movement  Report  (MMR)  and  judgment  to  forecast  major  end- 
item  movement  requirements.  There  is  limited  confidence  in  the  TAEDP  at 
the  ICPs;  consequently,  there  is  significant  reliance  on  the  MMR  and  fore¬ 
caster  judgment. 

i.  What  is  the  commodity  impact  of  AAFES  forecasting  procedures?  The 

AAFES  forecasting  procedures,  per  se,  do  not  impact  the  commodities.  AAFES 
accounted  for  30  percent  of  all  cargo,  50  percent  of  general  cargo,  and  60 
percent  of  reefer  cargo  shipped  in  FY  1982. 

j.  What  are  the  effective  and  cost  effective  methods  for  improving 
forecasting  accuracy?  .  Chapters  5  and  6  identify  procedures  and  methods 
that  would  improve  the  accuracy  of  the  forecasting  system. 


k.  What,  if  any,  unique  commodities  are  masked  as  general  or  special 
cargo? 

(1)  Much  of  the  CECOM  force  modernization  equipment  is  masked  as  gen¬ 
eral  cargo  and  is  not  currently  identifiable.  This  prevents  accurate 
adjustment  of  general  cargo  forecasts  to  reflect  force  modernization 
programs. 
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(2)  Military  sedans  and  smaller  military  wheeled  vehicles  are  classi¬ 
fied  as  and  are  masked  in  special  cargo.  All  special  cargo  is  forecast  as 
transported  via  MSC  breakbulk  ships.  These  vehicles  can  be  transported  via 
breakbulk  roll  on/roll  off  (RO/RO)  or  container  vessels,  and  could  be  iden¬ 
tified  as  a  separate  commodity  or  in  a  single  commodity  with  POVs. 


CHAPTER  8 


CONCLUSIONS 


8-1.  GENERAL.  The  investigation  of  the  current  Army  transportation  work¬ 
load  forecasting  system  and  the  methodologies  used  within  the  process  indi¬ 
cates  that  it  can  be  improved.  Significant  improvements  and  efficiencies 
can  be  gained  from  refining  the  Army  forecasting  system,  using  accurate 
forecasting  tools,  and  amending  the  reporting  requirements  of  the  transpor¬ 
tation  operating  authorities.  Annual  cost  avoidance  potential  should  range 
betwen  $30M  and  $100M. 

8-2.  FORECASTING  SYSTEM 

a.  The  Army  transportation  workload  forecasting  system  can  be  improved 
by  using  improved  forecasting  methodologies  and  revising  the  forecasting 
process. 

b.  A  sufficient  data  base  has  been  created  in  this  study  which  can  be 
used  as  a  basis  with  forecasting  techniques  for  making  accurate  forecasts 
and  for  incorporation  into  a  data  base  used  by  Army  forecasters. 

c.  The  Box-Jenkins  and  Winters  models  can  provide  accurate  forecasts  of 
the  Army  over-ocean  surface  cargo  transportation  requirements. 

d.  When  used  in  conjunction  with  program  information,  the  forecasts 
developed  using  the  Box-Jenkins  and  Winters  models  should  provide  the  tran¬ 
sportation  operating  agencies  (TOA)  with  accurate  forecasts  of  Army  over¬ 
ocean  surface  cargo  transportation  requirements. 

e.  The  forecasting  of  Army  transportation  requirements  by  the  TOA, 
using  a  forecasting  model  and  program  input  from  the  Army,  is  considered  to 
be  the  most  efficient,  accurate,  and  cost  effective  system. 

f.  Forecasting  by  a  single  Army  agency  would  achieve  similar  accuracy 
as  forecasting  by  the  TOA  and  would  also  result  in  significant  cost 
savings.  If  forecasting  by  the  TOA  is  not  achieveable  in  the  near  term, 
forecasting  by  a  single  Army  agency  is  the  most  efficient  alternative.  If 
forecasting  by  the  TOA  is  ultimately  developed,  then  the  transition  from 
forecasting  by  a  single  Army  activity  to  the  TOA  would  require  minimal  ef¬ 
fort. 

g.  If  the  present  forecasting  structure  is  to  be  retained,  accuracy 
would  be  improved  by  combining  long-range  and  short-range  forecasting  sys¬ 
tems,  using  an  accurate  multiyear  data  base,  applying  uniform  forecasting 
methods,  structuring  TACs  to  identify  each  shipper  and  program,  upgrading 
MECHTRAM,  and  developing  an  effective  feedback  program  to  inform  shippers 
on  the  cargo  shipped. 


h.  HQ  DARCOM  is  the  most  appropriate  location  for  the  forecasting  ac¬ 
tivity  if  it  is  decided  that  a  single  activity  within  the  Army  will  perform 
the  forecasting  function.  Predominant  reasons  for  this  conclusion  are: 

(1)  DARCOM  is  a  major  Army  shipper. 

(2)  HQ  DARCOM  is  cognizant  of  and  has  a  command  relationship  with  the 
program  activities  which  have  the  most  significant  impact  on  irregular 
shipping  requirements  (e.g.,  force  modernization  or  ammunition  program 
decisions).  Consequently,  coordination  problems  would  be  minimized. 

(3)  Non-DARCOM  shipments  are  accurately  predictable  and,  for  the  most 
part,  can  be  forecasted  using  a  forecasting  model. 

(4)  HQ  DARCOM  has  the  capability,  with  augmentation,  to  support  the 
forecasting  system. 

i.  HQDA  and  LCA  also  have  the  capability  to  support  the  system,  and 
could,  if  augmented  with  the  appropriate  skills,  perform  the  forecasting 
function. 

j.  TAC  structure  revision  would  allow  identification  of  program 
specific  cargo  within  the  general  and  special  commodity  categories  and 
facilitate  forecasting  of  these  commodities. 

k.  Interdependence  of  the  transportation  budgets  and  workload  forecasts 
would  be  enhanced  if  they  used  the  same  data  sources. 

l.  Accurate  and  timely  program  data  for  force  modernization, 
ammunition,  fuels,  and  personnel  related  programs  would  improve  forecast 
accuracy. 

m.  The  automated  systems  supporting  the  USAF  surface  tonnage  and  cost 
system,  and  the  forecast  of  Air  Force  transportation  requirements  process 
is  an  appropriate  paradigm  for  the  Army's  ADP  system  to  support  transpor¬ 
tation  workload  forecasting  when  used  with  the  forecasting  models  developed 
in  this  study.  It  is  also  an  .appropriate  model  for  the  cost  and  perfor¬ 
mance  aspects  of  the  system. 

n.  There  are  no  indications  that  separation  of  budget  and  forecasting 
responsibilities  contribute  to  forecasting  error. 

8-3.  EXTERNAL  FACTORS 

a.  The  requirement  to  forecast  the  mode  of  shipment  as  breakbulk  or 
container  is  a  significant  contributor  to  forecasting  error.  Resolution  of 
this  problem  with  MSC  and  MTMC,  resulting  in  elimination  of  this 
forecasting  requirement,  would  eliminate  at  least  one-third  of  the 
forecasting  error. 
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b.  Short-range  forecasts  appear  to  be  of  minimal  value  to  the  TOAs  and 
do  not  justify  the  amount  of  effort  required  for  their  preparation. 
Submission  of  annual  reports  and  periodic  change  reports  should  suffice. 

c.  MSC  billing  system  appears  to  contain  significant  inaccuracies  and 
delays  in  its  reports.  Stringent  audits  of  MSC  bills,  establishment  of  a 
history  of  accuracy  and  timeliness  of  the  MSC  billing  system,  and  develop¬ 
ment  of  procedures  with  MSC  to  resolve  system  defects  which  are  identified 
should  result  in  more  accurate  costs  and  an  accurate  data  base  for  forecast 
formulation. 

d.  The  amount  of  household  goods  processed  by  MSC  is  a  small  proportion 
of  the  total  household  goods  shipped.  MTMC  decides  whether  household  goods 
will  be  shipped  by  a  freight  forwarder  or  through  MSC.  As  a  small  shift 
from  freight  forward  shipment  to  MSC  would  have  a  major  impact  on  forecast 
accuracy,  MTMC  involvement  in  the  forecast  formulation  should  reduce  the 
potential  for  error  in  the  Army  household  goods  forecast. 


CHAPTER  9 

ACTIONS  TO  IMPROVE  THE  TRANSPORTATION  FORECASTING  SYSTEM 


9-1.  GENERAL.  This  chapter  outlines  those  actions  identified  in  this 
study  that,  if  taken,  should  improve  the  Army's  transportation  workload 
forecasting  process.  Regardless  of  the  system  chosen,  the  actions 
concerning  TAC  revision,  resolution  of  mode,  household  goods  and  short- 
range  forecasting  requirements,  refinement  of  forecasting  procedures,  and 
budget  and  forecast  interdependence,  if  executed,  would  improve  forecast 
accuracy.  It  is  assumed  during  the  implementation  of  suggested  actions: 

(1)  that  the  current  forecasting  system  will  be  maintained  until  the 
selected  system  is  totally  in  place,  and  (2}  that  a  timetable  and 
responsibilities  for  the  completion  of  necessary  actions  will  be 
established. 

9-2.  SELECT  FORECASTING  SYSTEM  STRUCTURE.  Selection,  from  the  proposed 
forecasting  system  structures,  of  that  system  or  combination  of  systems 
which  most  satisfactorily  meets  the  requirements  of  an  improved  forecasting 
system  is  a  necessary  prerequisite  to  the  implementation  of  other  changes. 
Forecasting  at  a  single  location,  using  historic  and  program  information, 
is  considered  the  most  efficient  and  effective. 

9-3.  RESOLVE  TAC  STRUCTURE.  Restructuring  of  the  TAC  to  include  TAC 
identification  of  the  major  program  under  which  the  commodity  is  shipped 
(i.e.,  Force  Modernization)  is  required  to  isolate  force  modernization.  If 
either  decentralized  forecasting  system  option  is  selected,  then  including 
identification  of  the  program  and  the  shipping  command/ agency  in  the  TAC 
structure  is  suggested. 

9-4.  RESOLVE  MODE  FORECASTING  REQUIREMENT.  Resolution  with  the  surface 
transport  TOA  to  eliminate  the  requirement  to  forecast  whether  cargo  is  to 
be  shipped  breakbulk  or  container  and  that  the  forecast  state  only  the 
commodity  and  route  is  necessary  to  eliminate  a  significant  source  of 
forecasting  error. 

9-5.  RESOLVE  SHORT-RANGE  FORECASTING  REQUIREMENT.  Resolving  with  MSC  and 
MTMC  to  eliminate  the  monthly  short-range  requirements  forecast  and 
replacing  it  with  periodic  reports  of  changes  to  the  annual  forecast  is 
suggested. 

9-6.  RESOLVE  HOUSEHOLD  GOODS  FORECASTING  PROCEDURE.  Establishment  of  a 
process  whereby  MTMC  will  participate  in  the  development  of  the  household 
goods  forecast  to  determine  what  proportion  of  household  goods  will  be 
shipped  as  Code  4  and  Code  5  should  reduce  the  potential  for  forecast 
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9-7.  ESTABLISH  REPORTING  REQUIREMENTS.  Establishing  reporting  require¬ 
ments,  based  on  the  system  selected  and  the  resolution  of  forecasting 
requirements  with  the  TOA,  will  ensure  understanding  of  revised  reporting 
responsibilities.  The  type,  content,  and  frequency  of  reports  should  be 
established  in  coordination  with  the  forecasting  agencies/commands,  budget 
agencies,  and  the  TOA. 

9-8.  DEFINE  LINKAGE  BETWEEN  FORECAST  AND  BUDGET.  Establishing  mechanisms 
to  ensure  that  the  forecasting  of  transportation  requirements  and  budget 
processes  are  coordinated  and  have  access  to  the  same  source  data  will 
assist  both  forecast  and  budget  preparation. 

9-9.  DEVELOP  ADP  SYSTEM  REQUIREMENTS.  This  action  is  also  a  function  of 
the  forecasting  system  selected.  If  the  present  system  or  variations 
thereof  are  selected,  then  modification  of  the  MECHTRAM  system  to  provide  a 
more  extensive  data  base,  evaluate  the  forecasts  made  by  the 
commands/agencies,  produce  information  of  greater  utility  to  the  budget 
process,  and  provide  timely  cost  and  performance  data  will  improve  the 
forecasting  process.  If  a  single  Army  agency  is  to  perform  the 
forecasting,  then  replacing  or  overhauling  the  MECHTRAM  system  is 
suggested.  In  either  of  these  cases,  including  in  the  system  those 
capabilities  now  embodied  in  the  AFLC  system  and  the  forecasting  models 
developed  in  this  study  will  facilitate  system  improvement.  If  the 
forecasts  are  to  be  prepared  by  the  TOA,  the  modification  of  the  MECHTRAM 
system  to  accept  program  input  and  produce  consolidated  program  reports, 
accurately  process  MSC  bills,  and  provide  comprehensive  cost  and 
performance  reports  on  the  surface  cargo  transportation  process  will 
provide  necessary  forecast  and  budget  information. 

9-10.  UPDATE  FORECASTING  MODELS.  Depending  on  the  time  between  completion 
of  the  models  by  the  study  agency  and  their  use  by  the  forecasters,  it  may 
be  necessary  to  refine  the  forecasting  models.  The  models  should  be 
checked  by  comparing  the  most  recent  lift  data  available  to  forecasts 
developed  by  the  model  for  the  time  periods  in  question. 

9-11.  REVISE  REGULATIONS.  Changes  to  AR  55-23  and  AR  55-30,  reflecting 
changes  to  the  system,  will  be  necessary  to  ensure  understanding  of  the 
revised  process.  Additionally,  change  recommendations  to  JCS  Pub  15  should 
reflect  the  changes  to  forecasting  procedures  which  may  be  developed  with 
the  TOA. 
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APPENDIX  B 


STUDY  DIRECTIVE 


DAL0-TSP-C11 


DEPARTMENT  OF  THE  ARMY 

OFFICX  OF  THE  DEFVTY  chief  of  staff  for  logistics 
WASHINGTON,  D.C.  20JI0 


11  May  1983 


SUBJECT:  Transportation  Workload  Forecasting  (TWF)  Study 


Director 

U.  S.  Army  Concepts  Analysis  Agency 
8120  Woodmont  Avenue 
Bethesda,  Maryland  208 14 


1.  Purpose  of  Study  Directive.  This  directive  tasks  the  Concepts  Analysis 
Agency  to  conduct  the  subject  study. 

2.  Study  Title.  Transportation  Workload  Forecasting  (TWF)  Study. 

3*  Background.  Current  forecasting  procedures  of  Army  cargo  and  mail 
workload  requirements  directed  by  AB  55-30  prescribe  input  from  seventeen  major 
commands/agencies/activities,  world-wide.  These  consolidated  requirements  are 
submitted  by  HQDA  to  the  Military  Sealift  Command  (MSC)  and  the  Military 
Airlift  Command  (MAC)  in  accordance  with  Joint  Chiefs  of  Staff  Publication  15. 
The  MSC  provides  the  Military  Traffic  Management  Command  (MTMC)  a  copy.  MTMC, 
MSC  and  MAC  utilize  this  data  to  generate  their  industrial  fund  budgets. 

History  reveals  significant  variances  in  forecasted  requirements  versus  actual 
lift,  which  result  in  distorted  budgets  by  both  the  shipper  service  and 
MSC/MAC/MTMC .  Transportation  funds  to  pay  these  overseas  movements  are  cen¬ 
trally  budgeted  at  HQDA,  with  the  0.  S.  Army  Finance  and  Accounting  Center 
(USAFAC)  the  designated  office. 

4.  Study  Proponent  and  Proponent's  Study  Director.  HQDA  ODCSLOG  is  the  study 
proponent.  Director  of  Transportation,  Energy  and  Troop  Support,  ODCSLOG 
(DALO-TSP) ,  will  be  the  Proponent's  study  representative. 

5.  'Study  Agency.  U.  S.  Army  Concepts  Analysis  Agency  (CAA). 

6.  Terms  of  Reference. 


a.  Statement  of  the  Problem.  Current  Army  transportation  workload 
forecasting  procedures  result  in  unrealistic  forecasts  of  Army  lift  require¬ 
ments  to  MAC,  MSC,  and  MTMC. 

b.  Purpose .  To  develop  procedures  to  improve  US  Army  transportation 
workload  forecasting. 

c.  Scope.  This  study  will  focus  on  the  long-range  surface  transportation 
workload  forecast  process  and  its  impact  on  the  Army  budget  and  on  transpor¬ 
tation  costs.  Short-range  forecasting  will  be  examined  to  the  extent  that  it 
impacts  on  or  influences  the  long-range  forecast. 
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d.  Objectives.  Determination  of  the  nature  and  extent  of  the  TWF  problem, 
exploration  of  alternative  solutions  to  the  problem,  evaluation  of  the  alter¬ 
native  solutions  in  terms  of  cost  and  feasibility,  and  development  of  an  out¬ 
line  plan  to  manage  the  TWF  system  improvements. 

e.  Tasks. 

(1)  Identify  the  nature  and  extent  of  U.  S.  Army  transportation 
workload  forecasting  variances  from  actual  lift. 

(2)  Determine  the  impact  of  long-range  transportation  workload 
forecasts  on  the  Army  budget,  on  the  industrial  funds  of  Military  Sealift 
Command  (MSC) ,  the  Military  Airlift  Command  and  Military  Traffic  Management 
Command  (MTJK) ,  and  on  rates  established  for  Army  second  destination  cargo 
movements. 

(3)  Examine  feasible  and  cost  effective  methods  for  increasing 
accuracy  in  forecasting. 

(4)  Examine  and  evaluate  alternative  locations  for  forecasting 
responsibility. 

(5)  Recommend  an  operationally  and  cost  effective  transportation 
workload  forecasting  methodology. 

f.  Timeframe.  Current. 

g.  Assumptions.  Transportation  workload  forecasting  requirements  will 
remain  unchanged  for  the  duration  of  the  study. 

h.  Essential  Elements  of  Analysis  (EEA) . 

(1)  What  are  the  recorded  variances  between  long-range  transportation 
workload  forecasting  and  actual  utilization  of  cargo  shipping? 

(2)  Do  systemic  conditions -exist  which  contribute  to  unrealistic 
forecasting?  If  so,  what  are  they? 

(3)  What  is  the  economic  and  operational  impact  of  long-range  fore¬ 
casts  which  are  at  variance  with  actual  utilization? 

(4)  Do  short-range  forecasts  impact  on  long-range  forecasting?  If  so 

how? 

(5)  What  is  the  impact  of  the  current  separation  of  responsibility 
for  long-range  and  short-range  forecasting? 
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(6)  What  methodologies  exist  in  the  other  services  which  could  be 
applied  to  the  resolution  of  the  Army  problem? 

(7)  What  DARCOM  activities  affect  major  items  of  equipment  planned 
for  oversea  distribution? 

(8)  What  is  the  impact  of  the  Total  Army  Equipment  Distribution 
Program  (TAEDP)  on  forecasting  of  requirements. 

(9)  What  is  the  commodity  impact  of  AAFES  forecasting  procedures? 

(10)  What  are  the  feasible  and  cost  effective  methods  for  improving 
forecasting  accuracy? 

(11)  What,  if  any,  unique  commodities  are  masked  as  general  or  special 

cargo? 

7.  Responsibilities . 

a.  The  0DCSL0G  will: 

(1)  Provide  support  as  required  for  its  areas  of  responsibility  and 
interest. 

(2)  Prepare  an  evaluation  of  study  results  IAW  AR  5-5. 

(3)  Establish  and  convene  a  Study  Advisory  Group  (SAG)  under  provi¬ 
sions  of  AR  5-5. 

b.  CAA  will: 

(1)  Establish  a  study  team. 

(2)  Establish  direct  communications  with  0DCSL0G,  D/TRETS,  and  other 
agencies  as  required  for  the  conduct  of  the  study. 

(3)  Provide  periodic  in-process  reviews  (IPR)  and  final  study  docu¬ 
mentation  to  the  study  sponsor. 

c.  USAMSSA  will:  Provide  ADP  support  as  requested. 

8.  Literature  Search. 

a.  Department  of  the  Army,  Office  of  the  Comptroller  of  the  Army,  Report 
on  the  Army  Transportation  Study,  May  1971. 
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b.  Evaluation  of  Second  Destination  Transportation  Funding,  U.  S.  Army 
Logistics  Evaluation  Agency,  29  December  1978. 

e.  Defense  Logistics  Agency  studies  and  reports. 

d.  MSC  reports. 

e.  Army  Inventory  Management  Agency  studies. 

f.  USAF  and  DSN  transportation  workload  forecasting  methodologies. 

g.  OSD  transportation  workload  forecasting  studies  and  reports. 

9*  References . 

a.  JCS  Pub  15,  dated  2  June  1975* 

b.  AR  55-23,  dated  17  March  1978. 

c.  AR  55-30,  dated  15  August  1982. 

d.  AR  55-133,  dated  18  February  1977. 

e.  AR  59-8,  dated  20  August  1982. 

f.  MECHTRAM  Users  Manual,  dated  June  1978. 

g.  AR  11-18,  October  1975. 

h.  DA  PAM,  May  1976. 

i.  AR  11-28,  December  1975. 

10.  Administration. 

a.  Support 

(1)  Funding  for  temporary  duty  (TDY)  and  travel  associated  with  the 
study  will  be  provided  by  each  participating  agency. 

(2)  Headquarters  or  agencies  represented  in  the  Study  Advisory  Group 
will  provide  crwn  TOY,  per  diem,  and  travel  funds. 

b.  Milestone  Schedule 

First  IPR  May  1983 
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Second  IPS 
Third  IPR 


September  1983 
December  1983 


o.  Control  Procedures 

(1)  Periodic  IPRs  will  be  provided  to  the  study  sponsor  by  the  study 

team. 

(2)  Docuaentation  required  by  AR  5-5,  including  DD  Form  1498  and  DD 
Form  1473,  and  a  final  report  to  include  an  Executive  Summary  will  be  sub¬ 
mitted  by  CAA. 

d.  Coordination 

(1)  Direct  coordination  between  CAA,  DALO-TSP,  and  forecasting  activ¬ 
ities  is  authorized.  For  purposes  of  any  possible  data  collection,  coor¬ 
dination  between  CAA  and  the  submitting  activity  is  directed.  Information 
copies  of  all  data  inputs  should  be  provided  to  HQDA,  DAL0-TSP-C1 1 . 

(2)  This  study  directive  has  been  coordinated  with  CAA  in  accordance 
with  AR  5-5. 

FOR  THE  DEPUTT  CHIEF  OF  STAFF  FOR  LOGISTICS: 


JIMMY  DA  ROSS 
Brigadier  General,  GS 
Director]  of  Transportation, 
S Energy  and  Troop  Support 
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APPENDIX  D 

COMMODITY  TIME  SERIES  DATA 


D-l.  GENERAL.  The  purpose  of  this  appendix  Is  to  provide  a  list  of  the 
aggregate  time  series  data  that  was  used  to  build  the  forecasting  models 
and  test  their  accuracy. 

D-2.  COMMODITY  DATA 

a.  Each  figure  contains  three  columns  of  data:  (1)  column  1  Identifies 
the  year  and  month,  (2)  column  2  depicts  the  amount  of  measurement  tons 
that  were  lifted  during  the  month,  and  (3)  column  3  lists  the  number  of 

•  records  or  shipments  that  were  aggregated  to  develop  the  lift  figure  in 
column  2. 

b.  As  discussed  in  paragraph  6-2,  aggregate  time  series  data  was 
developed  for  each  commodity  over  all  routes  as  well  as  each  commodity  for 
a  particular  route.  The  data  contained  in  this  appendix  is  the  commodity 
data  for  all  routes  combined.  The  specific  route  data  (424  routes)  is 
available  upon  request. 
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Figure  0-1.  Chill  Cargo  Data 
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Figure  0-2.  Freeze  Cargo  Data 
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Figure  D-3.  Coal  Cargo  Data 
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Figure  0-4.  POV  Cargo  Data 
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Figure  D-5.  Ammunition  Cargo  Data 
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Figure  D-6.  General  Cargo  Data 
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Figure  D-7.  HH6  Cargo  Data 
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Figure  D-8.  CONEX  Cargo  Data 


APPENDIX  E 


BOX-JENKINS  METHODOLOGY 


E-l.  GENERAL 

a.  Box-Jenklns  models  are  a  flexible  class  of  stable  linear  statistical 
models  that  are  used  to  model  stochastic  time  series  data.  The  models  are 
first  fitted  to  historical  time  series  data  using  parametric  estimates. 

The  fitted  model  is  then  used  to  provide  forecasts  of  future  observations. 


b.  This  appendix  will  discuss: 

(1)  Stochastic  properties  of  time  series  models  and  the  principle  of 
stationarity. 

(2)  Basic  model  autocorrelation  and  partial  autocorrelation 
functions. 

(3)  Seasonal  models. 

(4)  The  backshift  operator. 

(5)  Model  diagnostic  checks. 

E-2.  STOCHASTIC  PROPERTIES  AND  STATIONARITY 

a.  The  assumption  of  a  historical  stochastic  time  series  differentiates 
Box-Jenklns  models  from  other  time  series  models.  Some  time  series  models 
are  simple  extrapolation  formulas  which  fail  to  account  for  the  stochastic 
properties  of  time  series.  Box-Jenklns  models  assume  that  the  historical 
data  evolves  stochastically  and  attempt  to  model  the  stochastic  properties 
of  the  time  series  for  forecasting  purposes.5 

b.  A  stochastic  process  assumes  that  each  observation  xi,  xg,  ...  xt  is 
randomly  drawn  from  a  probability  distribution.  Using  Box-Jenkins  tech¬ 
niques,  an  attempt  is  made  to  duplicate  the  stochastic  process  in  hopes  of 
understanding  the  probability  distributions  of  future  observed  values  and 
providing  accurate  predictions  of  the  future.  A  perfect  duplication  of  the 
stochastic  process  is  practically  impossible.  This  would  require  a  joint 
distribution  of  all  possible  combinations  of  time  series  values  xj,  X2  ... 
xt-  If  the  time  series  is  large,  one  can  quickly  see  that  the  resulting 
probability  distribution  function  would  be  immense.  In  view  of  this  im¬ 
possible  task,  the  next  best  thing  is1  to  model  the  characteristics  of  the 
stochastic  process.  If  the  characteristics  can  be  identified,  then  the 
randomness  of  the  series  can  be  approximated  and  used  to  predict  future 
values.6 
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c.  Stochastic  properties  of  time  series  data  affect  model  development 
through  the  principle  of  stationarity.  Stationary  data  are  data  whose 
stochastic  processes  remain  constant  over  time.  If  the  stochastic  process 
is  constant  over  time,  i.e.,  stationarity  exists,  then  coefficients  similar 
to  those  used  in  regression  analysis  can  be  estimated  to  model  the  data. 

One  of  the  critical  assumptions  of  regression  analysis  is  the  assumption  of 
a  constant  linear  relationship  between  two  variables.  If  the  relationship 
remains  linear  over  time,  then  the  process  is  stationary,  and  the  model  is 
adequate.  If  not,  then  the  process  is  nonstationary  and  the  model  is  not 
correct. 

d.  The  principal  of  stationarity  assumes  that  the  stochastic  properties 
of  the  data  are  invariant  with  respect  to  time.  The  data  are  invariant  in 
that  the  stochastic  processes  of  the  data  are  constant  throughout  all  in¬ 
tervals  of  time.  The  stochastic  processes  are  assumed  to  be  in 
equilibrium,  and  any  variance  from  a  constant  mean  is  assumed  to  be  the 
same  at  any  point  in  time. 

e.  As  mentioned  earlier,  a  stochastic  process  assumes  that  each  obser¬ 
vation  is  randomly  drawn.  Future  predictions  of  the  stochastic  process 
would  be  based  on  the  conditional  probability  distribution  for  the  series. 
Thus,  prediction  of  xj+i  would  be  made  given  the  probability  distribution 
of  the  series  xi  ...  xj.  In  mathematical  terms,  this  conditional  proba¬ 
bility  would  be  expressed  as: 


P  (xj  +  l| xi,  X2  ...  xt) 

If  the  time  series  is  stationary,  the  joint  distribution  and  conditional 
distribution  will  remain  constant  throughout  time.  Thus, 

P(xt»  •••»  xt+n)  s  P(xt+m*  •••»  xt+n+m) 

P(xt)  =  P(xt+m) 

Finally,  if  the  series  is  stationary,  the  mean,  variance,  and  covariance  of 
the  data  will  remain  constant  over  time.6 

f.  The  autocorrelation  function  is  a  statistical  tool  that  is  used  in 
Box-Oenkins  models  to  describe  the  stochastic  process  of  a  time  series  and 
provides  an  understanding  of  the  probability  distributions  of  the  time 
series.  The  autocorrelation  function  \  is  defined  as: 
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If  a  stochastic  process  consists  of  independently  distributed  random 
variables  with  a  mean  of  zero,  then  the  autocorrelation  function  for  pq  =  1 
and  (>|(50  for  all  lags  k>0.  This  process  indicates  the  existence  of 
white  noise,  which  means  that  the  forecast  of  xj+]  =  0  for  all  1>0.  When 
white  noise  is  achieved,  random  errors  of  the  modeling  process  have  been 
eliminated  since  the  predicted  value  of  the  error  term  equals  zero. 

g.  As  in  other  statistical  models,  Box-Oenkins  uses  a  sample  of  time 
series  to  make  predictions  about  the  population.  Thus,  the  sample 
autocorrelation  function  is  defined  as: 


A 

p,. 


T-k 

(xt  •  ;><Vk  - *> 

- ' - 7T” 

1  (x.  -  x)Z 

t=l  1 


This  statistic  is  used  on  Box-Jenkins  models  to  compute  the  autocorrelation 
function  for  different  values  of  lag  k.  If  the  data  are  generated  by  a 
stationary  process,  the  autocorrelation  function  estimates  should  fall 
to  zero  quickly  as  k  increases.  A  failure  of  to  drop  off  quickly  to 
zero  indicates  the  existence  of  nonstationary  data.  To  test  whether 
successive  coefficients  of  the  autocorrelation  function  are  equal  to  zero, 
i.e.,  generated  by  white  noise,  one  would  employ  the  Bartlett  test. 

Bartlett  determined  that  if  the  series  was  generated  by  white  noise,  the 
sample  autocorrelation  coefficients  are  approximately  distributed  with  a 
normal  distribution  of  mean  of  zero  and  standard  deviation  l/VT,  where  T 
equals  the  number  of  observations  in  the  series.  Thus,  if  the  sample 
coefficients  fall  within  the  confidence  interval,  the  sample  coefficients 
are  assumed  to  be  zero.6 

h.  Finally,  if  the  data  are  not  generated  by  a  stationary  process,  a 
technique  known  as  differencing  will  help  achieve  stationary  data.  Dif¬ 
ferencing  is  defined  as: 


"t  =  *t  -  *t-l 


where  the  new  series  wt  is  then  analyzed  for  stationarity.  If  the  data  do 
not  exhibit  stationarity  with  first  order  differencing,  the  series  w^  can 
be  differenced  again.  Additionally,  if  the  variance  of  the  data  does  not 
remain  constant  over  time,  logarithm  transformations  of  the  data  and  pos¬ 
sible  differencing  of  the  transformed  data  should  be  attempted  to  achieve 
stationarity. 

E-3.  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION  FUNCTIONS 

a.  Autoregressive  (AR)  models  attempt  to  describe  the  process  xt  with  a 
weighted  sum  of  past  values  of  the  series  xt_n  and  a  random  disturbance 
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xt  =  *lxt-l  +  ^2xt-2  +  *•*  +  <*nxt-n  +  5  +  t  t 


where  are  the  AR  parameters  5  is  the  mean  of  the  series. 

The  autocorrelation  function  for  the  AR  models  is  used  to  help  determine 
the  number  of  lags  in  the  model: 


An  AR(1)  process  depicts  a  function  that  declines  geometrically  from  p0  = 
1.  An  AR(2)  process  can  be  portrayed  as  an  oscillating  or  sinusoidal  func¬ 
tion  that  dampens  geometrically.  The  partial  autocorrelation  functions  for 
AR  models  closely  resemble  the  pattern  for  MA  model  autocorrelation  func¬ 
tion  and  vice  versa.  Thus,  to  confirm  the  existence  of  an  AR(1)  model,  one 
would  expect  to  find  a  partial  autocorrelation  function  with  a  significant 
coefficient  at  lag  1  and  zeros  for  all  coefficients  with  lags  k  >1.  Simi¬ 
larly,  the  AR(2)  partial  autocorrelation  function  would  depict  two  signifi¬ 
cant  coefficients  and  then  zeros  for  all  lags  k>  2. 2 

b.  Moving  average  (MA)  models  attempt  to  describe  the  process  xt  by  a 
weighted  sum  of  current  and  lagged  disturbance  terms,  ‘t-n* 


xt  s  M+  *t  "  Vt-1  '  *2‘t-2  "  •“  '  Vt-n 


where  m  is  the  mean  of  the  series  and  0n  are  the  MA  parameters. 


The  autocorrelation  function  for  MA  models  will  depict  how  many  disturbance 
terms  should  be  used  in  the  model: 


■*> 
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Autocorrelation  functions  for  MA  models  will  portray  significant  coeffi¬ 
cients  which  correspond  to  the  number  of  disturbance  terms  to  include  in 
the  model.  Successive  coefficients  should  decline  to  zero  rapidly.  The 
partial  autocorrelation  functions  for  MA  models  will  depict  coefficients 
that  dampen  exponentially  to  zero  or  dampen  in  an  oscillating  manner  to 
zero. 

c.  Autoregressive-moving  average  (ARMA)  models  are  mixed  models  which 
attempt  to  describe  the  process  xt  as  a  function  of  past  values,  lagged 
random  disturbances,  and  a  current  disturbance  term. 


‘t  *  Vt-1  *  Vt-2  +  —  ♦  Vt-r.  +  *  +  «t  '  Vt-1  '  Vt-2  - 


•  •  «  ”  9  (, 

n  t-n 


The  autocorrelation  and  partial  autocorrelation  functions  for  ARMA  models 
can  be  depicted  in  many  combinations.  For  more  information,  reference 
Box-Oenkins*  "Time  Series  Analysis"*  for  a  complete  description  of  possible 
function  portrayals. 

E-4.  BACKSHIFT  OPERATOR 


a.  The  "backshift  operator"  is  the  use  of  a  mathematical  concept  to 
simplify  model  building.  The  backshift  operator  B  is  defined  as: 

Bxt  =  xt-l 


In  more  general  terms,  the  backshift  operator  is  defined  as  BnXt,  where  n 
relates  to  the  number  of  past  values  used. 4  The  following  will  explain  the 
use  of  the  backshift  operator: 
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xt-l  =  BXt 

xt-2  *  B2xt 

xt-3  =  B3Xt 
•  • 

•  • 

•  • 

xt-n  =  Bnx^. 

b.  All  ARMA  models  can  be  expressed  using  the  backshift  operator.  For 
example,  an  AR(1)  model  is  depicted  as  follows: 

xt  s  *lxt-l  +  «t 

since  BXt  *  xt-l  the  equation  can  be  written  as: 
xt  s  *lBxt  +  «t 
xt  “  ♦lBxt  =  et 
(1  -  «/»iB)xt  =  «t 

Similarly,  an  MA  model  is  depicted  as  follows: 
xt  *  «t  -  #l«t-l 

since  B  *  <t-l»  the  equation  can  be  written  as: 

xt  *  «t  -  *lB«t 

xt  =  (1  -  *lB)«t 

c.  A  more  complex  model  such  as  an  AR(2)  is  rewritten  as  follows: 
xt  *  *lxt-l  +  *2xt-2  +  «t 

since  BXt  *  xt-l  and  xt-2  =  Bxt-l  *  B2xt*  the  ecluat'*on  can  now  be  written 
as: 

xt  *  *lBxt  +  *2B2xt  +  *t 
xt  "  ♦lBxt  -  *2B2xt  3  «t 
(1  -  <*>iB  -  <>2b2)  xt  =  *t 
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Similarly,  an  MA(2)  model  is  depicted  as: 

*t  =  (1  -  *lB  -  *2B2)€t 

The  use  of  the  backshift  operator  is  especially  useful  in  expressing 
seasonal  models.  This  will  be  explained  in  the  next  section. 

E-5.  SEASONAL  NOOELS 

a.  Seasonal  Box-denkins  models  are  a  particular  class  of  models  which 
incorporate  multiplicative  properties.  Seasonal  data  are  identified 
through  autocorrelation  analysis  of  the  time  series  data.  If  seasonality 
exists,  autocorrelation  coefficients  will  peak  at  lags  12,  24,  and  36  and 
should  be  correlated  with  each  other.  The  existence  of  seasonality  in  the 
data  violates  the  principle  of  stationarity  and  must  be  removed  through 
differencing.  In  this  case,  the  resulting  series  wt  is  a  seasonally 
adjusted  series: 


«t  *  xt  -  xt-l2 


If  the  series  wt  is  stationary,  it  is  modeled,  and  if  it  is  nonstationary, 
the  series  wt  is  differenced. 

b.  Once  stationary  data  is  achieved,  two  distinct  models  are  combined. 
One  model  is  a  seasonal  model  which  captures  the  seasonal  correlation  be¬ 
tween  observations  a  year  apart.  The  second  model  is  a  nonseasonal  model 
which  explains  the  dependency  of  observations  within  a  given  year.  The 
following  equations  illustrate  this  process. 

(1)  A  seasonal  AR  model  is  described  as  follows: 

*t  a  *12xt-12  +  «t 

xt  =  ^12Bxt.ii  +  «t 


■r. 
r  v; 


»  *  w  . 


i . .  • . 
V.v 


& 


AV. 


titter 


xt  3  *12B12xt  +  «t 

xt  -  *12B12xt  =  «t 

(1  -  ngB^xt  s  at 


(E-l) 
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(2)  Likewise,  a  seasonal  MA  model  is  constructed  as  follows: 

xt  ■  «t  -  *12«t-12 
H  =  «t  -  ^12B«t-ll 


xt  8  «t  -  0B12«t 

xt  •  U  -  «12B12)«t  (E-2) 

c.  The  two  seasonal  models  above  now  explain  the  relationship  between 
the  observations  separated  by  a  year  (i.e.,  October  1978,  October  1979). 
However,  a  nonseasonal  relationship  exists  between  the  months  within  a  year 
(i.e.,  October  1978,  November  1978).  Thus,  the  error  term  at  in  the  sea¬ 
sonal  model  would  not  be  totally  independent.  A  relationship  exists  be¬ 
tween  the  seasonal  model  error  terms  that  must  also  be  explained  in  order 
to  completely  analyze  the  data.  Nonseasonal  MA  or  AR  models  are  used  to 
explain  this  relationship. 


(1  =  «t  (E-3)  AR  nonseasonal  model 

at  =  (1  -  <?iB)«t  (E-4)  MA  nonseasonal  model 

Combining  seasonal  (E-l)  and  nonseasonal  (E-3)  components  provides  the 
following: 

(1  -  4>iB)(l  -  <*i2B12)xt  -  «t 
where  seasonal  (E-2)  and  nonseasonal  (E-4)  components  are  AR. 
xt  -  (1  -  *lB)(l  -  *12B12)«t 

where  seasonal  and  nonseasonal  components  are  MA.  Models  may  also  be  mixed 
where  the  seasonal  and  nonseasonal  components  are  opposite. 

(1  -  <t> 12B*2)  xt  a  (1  -  0lB}«t  seasonal  AR 

nonseasonal  MA 

(1  -  $iB)  xt  s  (1  -  012B*2)«t  seasonal  MA 

nonseasonal  AR 
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E-6.  MODEL  DIAGNOSTICS 

a.  The  next-.to-last  step  in  model  building  is  diagnostic  checks  on  the 
fit  of  the  model.  There  are  two  basic  diagnostic  checks  that  are  performed 
once  the  Box-Jenkins  model  is  estimated. 

b.  The  first  test  involves  a  visual  comparison  of  the  autocorrelation 
function  of  the  model  to  the  autocorrelation  function  of  the  original 
series.  This  particular  check  is  a  very  subjective  test.  If  the  auto¬ 
correlation  functions  are  similar,  then  the  model  is  assumed  to  be  valid. 

If  the  autocorrelation  functions  are  different,  then  the  model  adequacy  is 
suspect. 

c.  The  second  test  involves  a  quantitative  analysis  of  the  residuals  of 
the  model.  Box-Jenkins  model  building  assumes  that  the  error  terms  are 
normally  distributed  with  a  mean  of  zero  and  a  variance  1/T.  If  the  re¬ 
siduals  of  the  model  are  characterized  by  those  properties,  the  residuals 
closely  resemble  the  properties  of  white  noise.  Statistical  results  by 

G.  E.  P.  Box  and  D.  A.  Pierce  have  developed  the  Box-Pierce  statistic  Q: 


K 


k=l 


where  T  is  the  number  of  observations,  k  is  the  number  of  lags,  and  P\^  is 
the  estimated  residual  at  each  lag  k.  The  statistic  is  approximately  dis¬ 
tributed  as  a  chi-square  distribution  with  k  -  p  -  q  degrees  of  freedom  (p 
is  the  AR  order  and  q  is  the  MA  order).  Although  any  lag  greater  than  5  is 
sufficient  for  this  test,  the  normal  rule  of  thumb  is  to  include  enough 
lags  to  have  at  least  20  degrees  of  freedom. 6 

d.  The  null  hypothesis  H0  for  the  Box-Pierce  statistic  is  that  the  re¬ 
siduals  are  not  correlated  with  each  other,  have  a  mean  of  zero  and  a  vari¬ 
ance  1/T,  i.e.,  the  properties  of  the  residuals  resemble  white  noise.  Once 
calculated,  the  Q  statistic  is  compared  to  values  of  a  chi-square  distribu¬ 
tion  for  a  given  number  of  degrees  of  freedom.  If  the  Q  statistic  is  less 
than  the  value  in  the  chi-square  table,  the  null  hypothesis  is  not  rejected 
for  a  given  significance  level.  If  the  Q  statistic  is  greater  than  the 
chi-square  value,  the  null  hypothesis  is  rejected  for  a  given  significance 
level.  For  example,  if  the  Q  statistic  is  11.1  for  20  degrees  of  freedom, 
then  the  null  hypothesis  would  not  be  rejected  at  the  95  percent  confidence 
level,  since  11.1  is  less  than  12.4. 
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BOX-JENKINS  COMMODITY  ANALYSIS 


F-l.  GENERAL.  The  purpose  of  this  appendix  is  to  detail  the  steps  of 
analysis  that  were  performed  during  the  Box-Jenkins  model  development  for 
each  commodity.  The  analysis  discussion  will  consist  of:  (1)  a  visual 
examination  of  the  raw  data  patterns  and  moving  average  patterns,  (2)  the 
achievement  of  stationarity  through  differencing  or  transformations,  (3)  an 
identification  of  model  type  through  analysis  of  the  autocorrelation  and 
partial  autocorrelation  functions,  (4)  the  estimation  of  model  parameters, 
and  (5)  the  verification  of  model  fit  for  e;.ch  commodity.  The  analysis 
discussion  of  commodities  appears  in  this  appendix  as  follows: 

a.  POV,  pp  F-l  through  F-14. 

b.  General  cargo,  pp  F-15  through  F-29. 

c.  HHG,  pp  F-30  through  F-42. 

d.  Coal,  pp  F-43  through  F-65. 

e.  Ammunition,  pp  F-66  through  F-79. 

f.  Special,  pp  F-80  through  F-93. 

g.  Cargo  trail er/CONEX,  pp  F-94  through  F-106. 

h.  Freeze,  pp  F-107  through  F-121. 

i.  Chill,  pp  F-122  through  F-134. 

F-2.  POV 

a.  The  raw  data  for  POVs  is  depicted  in  Figure  F-l.  In  its  raw  form, 
the  series  exhibits  a  seasonal  trend  with  the  peaks  occurring  in  the  summer 
months  (June- July)  and  the  troughs  occurring  in  the  winter  months  (January- 
February).  This  pattern  is  more  pronounced  in  the  3-month  and  6-month 
moving  average  schematics  (Figures  F-2  and  F-3).  However,  the  mean  of  the 
series  is  fairly  constant  over  time  as  displayed  in  Figure  F-4. 

b.  Autocorrelation  function  analysis  of  the  raw  series  confirms  the 
seasonality  of  the  series  as  noted  by  the  large  autocorrelation  coeffi¬ 
cients  at  lags  6,  12,  18,  24,  30,  and  36  (Figures  F-5  and  F-6).  The  strong 
seasonal  trend  violates  the  principle  of  data  stationarity;  therefore,  the 
series  was  differenced  at  a  lag  of  12  months  to  eliminate  the  seasonal 
trend.  Figure  F-7  depicts  the  seasonally  adjusted  series. 
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c.  The  autocorrelation  function  of  the  seasonally  differenced  series  is 
depicted  in  Figure  F-8.  Aside  from  the  high  autocorrelation  coefficients 
at  lags  12  and  24,  the  series  exhibits  random  behavior  and  is  considered  to 
be  stationary.  The  significant  autocorrelation  coefficients  at  lags  12  and 
24  suggest  an  MA  seasonal  model.  Examination  of  the  partial  autocorrela¬ 
tion  function  (Figure  F-9)  identifies  a  possible  AR  nonseasonal  model. 
Again,  without  considering  the  autocorrelation  coefficients  at  lags  12  and 
24,  the  autocorrelation  pattern  is  similar  to  an  AR  process.  In  sum,  the 
specific  model  applicable  to  the  POV  seasonally  adjusted  series  is  a  non¬ 
seasonal  ARIMA  (2,0,0)  model  and  a  seasonal  ARIMA  (0,1,1)  model: 


(1  -  *1  B  -  *2  b2)  U  -  Bl2)  Xt  *  (1  -  *12  Bl2)«t 

d.  The  results  of  this  model  (Figure  F-10)  are  as  follows: 


(1  -  .0221B  -  .3554  02)  (1  -  Bl2)  Xt  -  (1  -  .8581  Bl2)n 
(.19)  (3.05)  (16.77) 

x*  (3,20)  =  16.35 


The  results  indicate  that  all  the  estimated  parameters  of  the  model,  except 
</>  i  are  significant  at  the  5  percent  significance  level.  Verification  of  the 
model's  adequacy  using  the  Box-Pierce  test  indicates  that  the  autocorrela¬ 
tions  of  the  estimated  model  residuals  are  not  significantly  different  from 
zero  at  the  32  percent  significance  level  (Figure  F-ll).  Thus,  the  null 
hypothesis,  residuals  are  uncorrelated,  should  not  be  rejected  at  the  32 
percent  significance  level. 

e.  However,  the  significant  autocorrelation  coefficent  at  lag  24  of  the 
estimated  autocorrelation  function  suggests  that  the  model  is  not  correctly 
specified.  Attempts  to  adjust  the  seasonal  model  to  an  ARIMA  (0,1,2)  re¬ 
sulted  in  a  nonstationary  model  since  #i>1.0.  Adding  another  seasonal 
model  to  the  original  model  resulted  in  the  following: 


(1  -  <t>\  B  -*2  B2)  (1  -  Bl2)  Xt  -  (1  -  *12  Bl2)  (1  -  024  B24)‘t 


This  model  eliminated  the  significant  coefficient  at  lag  24,  however  this 
model  is  more  restrictive  and  the  MSE  of  the  FY  82  forecast  was  larger  than 
the  (2,0,0)  x  (0,1,1)  seasonal  model.  These  model  adjustments  should  be 
tracked  during  future  forecasts  to  determine  if  they  are  warranted. 

f.  Using  this  model  to  forecast  FY  82  resulted  in  a  forecast  of  501.37 
K/MTON  versus  an  actual  shipment  of  510.252  K/MTON.  The  MSE  of  the  POV  fore¬ 
cast  for  FY  82  was  8.54  (Figure  F-12),  and  the  annual  forecast  error  is 
-1.7  percent. 
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Figure  F-6.  PACF  -  Raw  Data  -  POV 
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Figure.  F-8.  ACF  -  Lagged  12  Months 
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Figure  F-10.  Model  Parameters  -  POV 
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Figure  F-ll.  Estimated  Residuals  -  POV 
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Figure  F-12.  Forecast  FY  82  -  POV 
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F-3.  GENERAL  CARGO 

a.  The  raw  data  for  general  cargo  is  depicted  in  Figure  F-13.  The  raw 
data  appears  to  follow  a  seasonal  pattern  which  falls  in  the  early  part  of 
the  fiscal  year  and  rises  at  the  end  of  the  fiscal  year.  This  seasonal 
pattern  is  more  clearly  manifested  in  the  3-month  and  6-month  moving  ave¬ 
rage  diagrams  (Figures  F-14  and  F-15).  However,  it  should  be  noted  that 
the  seasonal  pattern  is  not  repeated  regularly  and  that  a  general  trend 
exists  in  the  data.  This  fact  is  verified  by  the  12-month  moving  average 
diagram  (Figure  F-16)  which  depicts  a  possible  5-year  cycle  in  the  data. 
Currently,  there  is  not  enough  data  to  confirm  the  existence  of  this  cycle 
throughout  time;  thus,  this  hypothesis  requires  continued  analysis. 

b.  Autocorrelation  function  analysis  of  the  data  confirms  the  existence 
of  seasonality  in  the  data  (Figures  F-17  and  F-18)  as  noted  by  the  large 
autocorrelation  coefficients  at  lags  6,  12,  18,  24,  30,  and  36.  To  elimi¬ 
nate  seasonality  in  the  data,  the  raw  series  was  differenced  at  a  lag  of  12 
months.  The  plot  of  the  new  series  is  shown  in  Figure  F-19.  It  is  readily 
obvious  that  this  series  is  not  stationary,  as  noted  by  the  downward  trend 
of  observations  over  time.  Nonstationarity  in  the  data  is  more  clearly 
evidenced  by  the  autocorrelation  function  (Figure  F-20  )  which  fails  to 
fall  quickly  to  zero.  In  an  effort  to  achieve  stationarity,  the  seasonally 
adjusted  series  was  differenced  at  a  lag  of  1  month.  The  resulting  auto¬ 
correlation  function  (Figure  F-21)  indicates  that  a  weak  stationary  series 
has  been  generated. 

c.  Analysis  of  the  autocorrelation  function  identifies  a  possible  mixe*' 
ARMA  nonseasonal  model  and  a  possible  MA  seasonal  model.  The  nonseasonal 
mixed  ARMA  model  is  suggested  based  upon  the  gradual  decline  (e*''«pt  lag  3) 
of  the  autocorrelation  and  partial  autocorrelation  coefficients  (Figures 
F-21  and  F-22).  The  autocorrelation  function  suggests  a  MA  seasonal 
process  as  manifested  by  the  gradual  decay  of  autocorrelations  at  lags  12 
and  24.  Thus,  the  data  suggest  a  multiplicative  seasonal  with  a  mixed 
ARIMA  nonseasonal  (1,1,1)  model  and  a  seasonal  ARIMA  (0,1,1)  model: 


(1  -*i  B)  (1  r  Bl2)  (l  -  B)  Xt  =  (l  -  el  B)  (1  -012  t 


d.  The  results  of  the  model  (Figure  F-23)  are  as  follows: 


(1  -  .3579  B)  (1  -  B12)  (1  -  B)  Xt  -  (1  -  .7962  B)  (1  -  .7753 

(2.5)  (8.31)  (13.43) 

X2  (3,20)  =  12.8 


F-15 
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The  results  indicate  that  all  of  the  estimated  parameters  are  signficant  of 
the  5  percent  significance  level.  The  Q  statistic  for  20  degrees  of  free¬ 
dom  is  12.8  (Figure  F-24).  Thus,  the  null  hypothesis,  estimated  model  re¬ 
siduals  are  uncorrelated,  should  not  be  rejected  at  the  13  percent  sig¬ 
nificance  level. 

e.  The  general  cargo  seasonal  model  forecasted  that  1856.61  K/MTON  would 
be  transported  during  FY  82  versus  an  actual  shipment  of  2021.55  K/MTON  The 
MSE  of  the  general  cargo  forecast  for  FY  82  was  318.54  (Figure  F-25),  and 
the  annual  forecast  error  was  -8.2  percent.  As  stated  in  paragraph  6-7, 
another  model  for  general  cargo  forecasted  FY  82  shipments  with  a  -4.6 
percent  error.  Both  models  should  be  tracked  into  the  future  for 
forecasting  accuracy  and  model  fit. 
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Figure  F-15.  Six-month  Moving  Average  -  General 
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Figure  F-16.  Twelve-month  Moving  Average  -  General 
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Figure  F-18.  PACF  -  Raw.  Data  r  General 
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Figure  F-20.  ACF  -  Lagged  12  Months 
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Figure  F-24.  Estimated  Residuals  -  General 
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Figure  F-25.  Forecast  FY  82  -  General 
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F-4.  HH6 

a.  The  raw  series  for  HH6  shipments  is  depicted  in  Figure  F-26.  The 
raw  data  depicts  the  volatility  of  HHG  shipments  when  MSC  is  required  to 
transport  normal  Code  4  shipments.  During  FY  78  HHG  shipments  by  MSC  were 
more  than  three-times  as  large  as  any  other  fiscal  year  in  the  data  set. 

Due  to  this  distortion  in  the  data,  only  data  from  October  1978  through 
September  1981  was  used  to  develop  the  models.  The  HHG  time  series  has 
only  36  observations,  much  less  than  is  normally  required  to  achieve  satis¬ 
factory  results  from  time  series  analysis. 

b.  Aside  from  the  distortion  during  FY  78,  the  raw  data  depict  a  sea¬ 
sonal  trend  with  peaks  during  the  summmer  months  and  troughs  during  the 
winter  months.  Evidence  of  seasonality  in  the  data  is  also  seen  in  the  6- 
month  moving  average  (Figure  F-27).  However,  the  12-month  moving  average 
(Figure  F-28)  suggests  that  the  mean  of  the  series  is  fairly  constant  over 
the  latter  years  and  appears  to  exhibit  a  slight  downtrend. 

c.  Autocorrelation  function  analysis  (Figures  F-29  and  F-30)  confirms 
the  existence  of  seasonality  in  the  data  as  manifested  by  the  significant 
autocorrelation  coefficient  at  lags  12  and  24.  The  series  was  differenced 
at  a  lag  of  12  months  to  eliminate  the  seasonality  and  achieve  stationary 
data.  The  plot  of  the  seasonally  adjusted  series  is  depicted  in  Figure 
F-31.  Autocorrelation  and  partial  autocorrelation  function  analysis  of  the 
data  (Figures  F-32  and  F-33)  indicate  that  stationarity  has  been  achieved. 

d.  Analysis  of  Figures  F-32  and  F-33  suggest  that  the  data  be  modeled 
with  a  seasonal  and  nonseasonal  model.  In  contrast  to  all  of  the  other 
models,  the  seasonal  model  for  HHG  was  AR  in  nature.  This  was  discovered 
after  comparing  the  results  of  AR  and  MA  seasonal  models.  The  nonseasonal 
model  was  identified  as  an  MA  process.  The  resulting  model  was  a 
nonseasonal  ARIMA  (0,0,1)  model  and  a  seasonal  ARIMA  (1,1,0)  model: 

(1-  <*>12  B12)  (1-B12)  Xt  *  m  +  (1-  04  B4)«t 

The  nonseasonal  model  was  built  with  only  one  parameter  rather  than  four 
since  the  first  three  coefficients  lacked  any  statistical  significance. 
Also,  the  mean  of  the  seasonally  adjusted  series  was  significant  at  the  6 
percent  significance  level. 

e.  The  results  of  the  model  (Figure  F-34)  were  as  follows: 

(1  +  .3198B12)  (1-b12)  Xt  ■  -1.163  +(l-.3397B4)e t 
(-2.33)  (-3.67)  (1.88) 

X2  (2,20)  =  6.26 


All  of  the  estimated  parameters  except  04,  are  significant  at  the  5  per¬ 
cent  significance  level.  04  is  significant  at  the  7  percent  signifi¬ 
cance  level.  Verification  of  the  model's  adequacy,  using  the  Box-Pierce 
test  with  20  degrees  of  freedom,  indicates  that  the  null  hypothesis  (uncor¬ 
related  residuals)  should  not  be  rejected  at  the  .5  percent  significance 
level  (Figure  F-35). 

f.  This  model  forecasted  FY  82  HHG  shipments  to  be  80.2  K/MTON  versus 
an  actual  lift  of  87.7  K/MTON  (Figure  F-36).  The  MSE  of  the  HHG  forecast 
was  5.24  percent,  and  the  annual  forecast  error  was  -8.5  percent. 
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Figure  F-27.  Six-month  Moving  Average  -  HHG 


Figure  F-28.  Twelve-month  Moving  Average  -  HHG 
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Figure  F-29 ..  ACF  -  Raw  Data  -  HHG 
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Figure  F-30.  PACF  -  Raw  Data  -  HHG 
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Figure  F-32.  ACF  -  Lagged  12  Months 
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Figure  F-33.  PCAF  -  Lagged  12  Months 
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Figure  F-35.  Estimated  Residuals  -  HHG 
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Figure  F-36.  Forecast  FY  82  -  HHG 
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F-5.  COAL 

a.  The  raw  series  for  commodity  coal  is  depicted  in  Figure  F-37.  The 
raw  data  do  not  depict  a  recurring  pattern  or  seasonal  trend.  Examination 
of  the  3-month  and  6-month  moving  average  diagrams  (Figures  F-38  and  F-39) 
affirms  the  absence  of  a  predictable  pattern.  The  12-month  moving  average 
diagram  (Figure  F-40)  suggests  the  possibility  of  a  5-year  cycle;  however, 
there  is  not  enough  data  to  verify  this  hypothesis. 

b.  Autocorrelation  function  analysis  of  the  data  confirms  the  absence 
of  seasonality  in  the  data  (see  Figures  F-41  and  F-42).  However,  the  auto¬ 
correlation  function  of  the  raw  series  indicates  that  a  trend  does  exist  in 
the  data  and  that  stationarity  is  not  present.  To  achieve  stationarity, 
the  raw  series  was  differenced  at  a  lag  of  one  month.  The  plot  of  the  dif¬ 
ferenced  series  is  illustrated  in  Figure  F-43.  Figures  F-44  and  F-45  de¬ 
pict  the  autocorrelation  and  partial  autocorrelation  functions  of  the  dif¬ 
ferenced  series.  These  functions  as  well  as  a  visual  examination  of  the 
plot  of  the  differenced  series  indicate  the  existence  of  stationary  data. 

c.  The  autocorrelation  and  partial  autocorrelation  functions  of  the 
differenced  series  suggest  that  the  data  could  be  modeled  using  an  ARIMA 
(0,1,1)  process.  Only  the  first  lag  of  the  autocorrelation  function  is 
significant,  and  the  remaining  coefficients  are  distributed  randomly  within 
the  95  percent  confidence  interval.  The  partial  autocorrelations  decay 
exponentially  to  zero  as  the  lags  increase.  Besides  an  ARIMA  (0,1,1) 
model,  the  data  could  also  be  modeled  using  an  ARIMA  (1,1,1)  model.  Thus, 
the  appropriate  models  are: 


(1  -  B)  Xt  =  (1  -  <?i)e  t 
or 

(i  -0i  b)  (i  -  b)  xt  =  (i  - n)u 


d.  Test  of  both  models  indicated  that  the  ARIMA  (1,1,1)  model  was  the 
most  appropriate  model.  Thus,  the  resulting  forecast  model  (Figure  6-46) 
is: 

(1  -  .34618)  (1  -  B)  Xt  =  (1  -  .859B)«t 
(9.46)  (2.06) 

x2  (2,20)  =  10.87 

e.  All  of  the  estimated  parameters  were  significant  at  the  5  percent 
significance  level.  The  model  was  verified  using  the  Box-Pierce  test  and  a 
visual  examination  of  the  estimated  residuals  of  the  model.  The  Q  statis¬ 
tic  test  is  10.87,  which  indicates  that  the  null  hypothesis,  uncorrelated 
residuals,  should  not  be  rejected  at  the  5  percent  significance  level 
(Figure  F-47).  Also,  all  of  the  estimated  residuals  fall  within  the  95 
percent  confidence  .interval  and  do  not  depict  any  noticeable  trend  or 
pattern. 
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f.  The  ARIMA  (1,1,1)  model  forecasted  the  FY  82  coal  tonnage  shipped  to 
be  406.51  K/MTON  (Figure  F-48)  as  compared  to  the  historical  FY  82  lift  of 
375.94  K/MTON  and  the  actual  FY  82  OLA  coal  contract  of  440.67  K/MTON.  The 
MSE  of  the  coal  forecast  was  192.60,  and  the  annual  forecast  error  was  -.6 
percent  (actual  FY  82)  and  15.2  percent  (DLA). 

g.  As  discussed  earlier,  several  external  factors  have  affected  coal 
shipments  during  October  1977  to  September  1981.  A  close  examination  of 
the  raw  data  indicates  that  contract  negotiation  delays  have  seriously  af¬ 
fected  coal  shipments  and  that  a  detailed  adjustment  of  coal  data  is  needed 
to  develop  an  accurate  forecasting  model.  The  need  for  data  adjustment  is 
clearly  identified  by  the  large  outlier  observations  contained  in  the  dif¬ 
ferenced  series  (Figure  F-43). 

h.  In  particular,  the  data  point  for  April  1980  is  a  large  outlier  that 
affects  the  model-bui Iding  process.  The  observation  for  April  1980  was 
adjusted  to  the  average  shipment  for  April  during  the  5-year  period.  The 
autocorrelation  and  partial  autocorrelation  functions  (Figures  F-49  and 
F-50)  of  the  adjusted  series  again  depict  a  trend  in  the  data.  Also,  it 
should  be  noted  that  almost  all  of  the  autocorrelation  coefficients  for  the 
adjusted  data  are  greater  than  the  autocorrelation  coefficients  for  the 
original  series  (Figure  F-41).  Thus,  the  data  are  very  sensitive  to  data 
manipulation,  and  further  data  adjustments  should  be  done  with  extreme 
caution. 

i.  The  differenced  adjusted  series  is  illustrated  in  Figure  F-51.  The 
autocorrelation  and  partial  autocorrelation  functions  (Figures  F-52  and 
F-53)  indicate  the  existence  of  weak  stationarity  in  the  data.  Again,  it 
should  be  noted  that  although  the  autocorrelation  function  of  the  differ¬ 
enced  adjusted  series  indicates  a  possible  ARIMA  (1,1,1)  model,  the 
function  varies  from  the  differenced  unadjusted  series.  In  particular, 
lags  10,  17,  and  24  are  more  significant  than  before.  This  would  suggest 
that  additional  parameters  need  to  be  added  to  the  original  model.  To 
retain  some  parsimony  in  the  model  only  one  parameter  was  added.  The  form 
of  the  model  is: 

(1  -  *1  B)  (1  -  B)  Xt  =  (1  -  n  B  -  8 10  BlO)«t 

j.  The  results  of  the  model  (Figure  F-54)  are  as  follows: 

(1  -  .21688)  (1  -  B)  Xt  =  (1  -  . 7693B  +  .2705B10)«t 
(1.46)  „  (9.69)  (-3.8) 

X2  (3,20)  =  10.51 

All  of  the  estimated  parameters  are  significant  at  the  5  percent  signifi¬ 
cance  level  except  The  Q  statistic  indicates  a  good  fit  for  the  model, 
thus  the  null  hypothesis  should  not  be  rejected  at  the  5  percent  signifi¬ 
cance  level  (Figure  F-55).  It  should  be  noted  that  the  "goodness  of  fit" 
test  accounts  for  the  significant  residual  coefficient  at  lag  17.  Attempts 
to  reduce  the  significance  of  this  coefficient  added  further  complication 
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and  restriction  to  the  original  model.  Also,  until  all  of  the  data  are 
adjusted  correctly,  additional  modifications  to  the  original  model  should 
not  be  attempted. 

k.  The  adjusted  model  forecasted  that  398.80  K/MTON  of  coal  would  be 
shipped  during  FY  82  (Figure  6-56)  which  corresponds  to  an  annual  forecast 
error  of  6.1  percent. 
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Figure  F-41.  ACF  -  Raw  Data  -  Coal 
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Figure  F-42.  PACF  -  Raw. Data  -  CocM 
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Figure  F-44.  .ACF  -  Laaged  1  Month 
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Figure  F-45.  PACF  -  Lagged  1  Month 
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Figure  F-47.  Estimated  Residuals  -  Coal 
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Figure  F-50.  PACF. -'Raw  Data  -  Coal  Adjusted 
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Figure  F-53.  PACF  -  Lagged  1  Month 
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Figure  F-54.  Model  Parameters  -  Coal  Adjusted 
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Figure  F-56.  Forecast  FY  82  -  Coal  Adjusted 
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F-6.  AMMUNITION 

a.  The  raw  series  for  ammunition  is  depicted  in  Figure  F-57.  The  raw 
series  depicts  a  recurring  pattern  with  the  troughs  occurring  in  the  first 
two  quarters  of  each  fiscal  year  and  the  peaks  in  the  latter  two  quarters. 
This  pattern  is  more  identifiable  in  the  3-month  and  6-month  moving  average 
diagrams  (Figures  F-58  and  F-59).  However,  a  visual  examination  of  the  12- 
month  moving  average  (Figure  F-60)  indicates  that  the  amount  of  ammunition 
shipped  from  year  to  year  has  remained  relatively  constant  with  a  slight 
trend  downwards. 

b.  Autocorrelation  function  analysis  of  the  model  data  confirms  the 
evidence  of  a  trend  in  the  data  and  suggests  the  possibility  of  seasonal 
data  (Figure  F-61).  The  seasonality  in  the  data  was  eliminated  by  dif¬ 
ferencing  the  series  at  a  lag  of  12  months.  The  autocorrelation  function 
of  the  resulting  series  is  depicted  in  Figure  F-62.  Stationarity  in  the 
data  has  not  been  achieved  as  evidenced  by  an  increasing  trend  in  the  auto¬ 
correlation  function  for  lags  1-4.  Attempts  to  achieve  stationarity 
through  differencing  or  other  transformations  of  the  seasonally  adjusted 
series  were  futile.  Therefore,  the  original  series  was  differenced  at  a 
lag  of  one  month  and  the  plot  is  depicted  in  Figure  F-63. 

c.  The  autocorrelation  and  partial  autocorrelation  of  the  differenced 
series  are  illustrated  in  Figures  F-64  and  F-65.  Analysis  of  the 
autocorrelation  function  indicates  that  an  ARIMA  (2,1,0)  model  is 
appropriate  to  model  the  data.  This  model  was  tested  and  the  resulting  fit 
of  the  model  was  not  good  due  to  significant  autocorrelation  coefficients 
at  lags  4  through  7.  To  improve  the  fit  of  the  model,  MA  parameters  were 
added.  The  resulting  model  was  an  ARIMA  (2,1,4)  model  of  the  form: 

(1-  01  B  -  <t> 2  B2)  (1  -  B)  Xt  =  (1  -  04  B4  -  05  B$  -  96  B6  -  97  B7)«t 

d.  The  results  of  the  model  (Figure  F-66)  are  as  follows: 

(1+.8418B  +  .5399  B2)  (1-B)  Xt  =  (l+.3454B4-.3689B5-.3875B6-.25617)t t 
(-7.33)  (-4.59)  (-3.11)  (3.86)  (4.57)  (2.43) 

X2  (6,20)  =  17.61 


All  of  the  estimated  parameters  are  significant  at  the  5  percent  signif¬ 
icance  level.  The  Q  statistic  for  the  model  (Figure  F-67)  indicates  that 
the  null  hypothesis,  estimated  residuals  are  uncorrelated,  should  not  be 
rejected  at  the  40  percent  significance  level. 


F-66 


■V 
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e.  The  model  forecasted  that  ammunition  shipments  during  FY  82  would  be 
189.62  K/MTON  (Figure  F-68).  The  actual  lift  of  ammunition  during  FY  82 
was  179.31.  The  MSE  of  the  ammunition  forecast  was  119.91  and  the  annual 
forecast  error  was  5.7  percent.  It  should  be  noted  that  this  model,  used 
in  conjunction  with  a  slightly  modified  model,  resulted  in  an  annual 
forecast  error  of  .3  percent. 

f.  Problems  with  attaining  stationarity  in  the  data  are  in  part  attrib¬ 
uted  to  outliers  in  the  data  series.  The  amounts  of  tonnage  shipped  during 
February  1980,  November  1980,  and  January  1981  are  far  below  the  average 
amount  of  ammunition  shipped  during  a  given  month.  Like  coal,  the  ammuni¬ 
tion  data  series  should  be  monitered  and  adjusted  with  extreme  care.  An 
attempt  was  made  to  modify  the  data  series,  however,  additional  work  in 
this  area  is  required  before  adequate  modeling  can  proceed. 
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Figure  F-59.  Six-month  Moving  Average  -  Ammunition 
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Figure  F-60.  Twelve-month  Moving  Average  -  Ammuniti 
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Figure  F-61.  ACF  -  Raw  Data  -  Ammunition 
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Figure  F-62.  ACF  -  Lagged  12  Months 
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Figure  F-64.  ACF  -  Lagged  1  Month 
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Figure  F-65.  PACF  -  Lagged  1  Month 
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ESTIMATION  by  bacxcasting  method 
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Figure  F-66.  Model  Parameters  -  Ammunition 
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Figure  F-67.  Estimated  Residuals  -  Ammunition 
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Figure  F-68.  Forecast  FY  82  -  Ammunition 


F-7.  SPECIAL 


a.  The  raw- series  for  the  special  cargo  commodity  is  illustrated  in 
Figure  F-69.  The  raw  series  does  not  depict  a  discernible  pattern  or  trend 
in  the  data.  Also,  examination  of  the  3-month  and  6-month  moving  averages 
does  not  identify  a  visible  pattern  (Figures  F-70  and  F-71).  The  12-month 
moving  average  (Figure  F-72)  suggests  that  a  general  trend  may  exist  in  the 
data  which  could  be  interpreted  as  the  program  funding  of  the  defense  bud¬ 
get.  In  fact,  the  trend  shows  a  gradual  decline  during  the  period  under 
study. 

b.  Autocorrelation  function  analysis  of  the  data  (Figures  F-73  and 
F-74)  suggests  the  possibility  of  seasonal  data.  Therefore,  the  series  was 
differenced  at  a  lag  of  12  months  in  an  effort  to  achieve  stationarity. 

The  variance  of  the  seasonally  differenced  series  exhibited  heteroscedas- 
ticity  (change  in  variance  over  time).  Thus,  the  original  series  was 
transformed  using  natural  logarithms  and  then  differenced  at  a  lag  of  12 
months.  The  differenced  logarithmetic  series  is  found  in  Figure  F-75.  The 
autocorrelation  and  partial  autocorrelation  functions  of  the  data  (Figures 
F-76  and  F-77)  indicates  that  a  random  stationary  series  has  been 
generated. 

c.  Analysis  of  the  autocorrelation  and  partial  autocorrelation  func¬ 
tions  indicates  that  the  data  could  be  modeled  as  a  possible  ARIMA  (2,0,0) 
nonseasonal  process  and  an  ARIMA  (0,1,1)  seasonal  model.  Thus,  the  appro¬ 
priate  model  is: 


(1  -*i  B  -  02B2)  (1  -  B1?)  lnxt  =  (1  -  B 12  Bl2)«t 

d.  The  results  of  this  model  (Figure  F-78)  are  as  follows: 

(1  -  .1496B  -  .387562)  (1  -  Bl2)  lnXt  =  (1  -  .801B12  -  .1503Bl4)*t 

(1.21)  (3.16)  (10.61)  (1.65) 

x2  (4,20)  =  13.15 

All  of  the  estimated  parameters  except  <t>i  and  #14  are  significant  at  the  5 
percent  significance  level.  The  additional  MA  parameter  (  0^4)  was  added 
to  the  model  since  the  autocorrelation  coefficient  at  lag  14  of  the  origi¬ 
nal  hypothesized  seasonal  model  was  highly  significant.  The  model  was  ver¬ 
ified  using  the  Box-Pierce  test  and  a  visual  check  of  the  estimated  re¬ 
siduals.  The  Q  statistic  for  this  model  with  20  degrees  of  freedom  is 
13.15  (Figure  F-79).  This  diagnostic  indicates  that  the  residuals  are  not 
significant  and  that  the  null  hypothesis  should  not  be  rejected  at  the  15 
percent  significance  level.  All  of  the  residuals  are  within  the  95  percent 
confidence  interval. 


CAA-SR-84-2 


e.  This  model  forecasted  that  559.42  K/MTON  of  special  cargo  would  be 
shipped  during  FY  82  (Figure  F-80)  versus  an  actual  amount  of  489.83  K/MTON 
The  annual  forecast  error  was  12.4  percent  and  the  forecast  MSE  was  159.23. 
Although  this  discrepancy  would  cause  one  to  doubt  the  usefulness  of  this 
model,  it  should  be  noted  that  the  amount  of  special  cargo  tonnage  shipped 
during  FY  82  is  much  lower  than  any  of  the  preceding  years.  In  fact,  the 
average  tonnage  of  special  cargo  shipped  over  the  past  5  years  was  669.17 
K/MTON.  In  all  cases,  forecasting  methodologies  are  employed  based  upon 
the  assumption  that  the  behavior  of  future  observations  will  not  differ 
greatly  from  the  behavior  of  past  observations.  When  the  future  is 
radically  different  than  the  past,  as  in  the  case  of  special  cargo  tonnage 
for  FY  82,  forecasts  cannot  predict  the  future  with  any  degree  of  accuracy. 


Figure  F-71.  Six-month  Moving  Average  -  Special 
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Figure  F-73.  ACF  -  Raw  Data  -  Special 
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Figure  F-74.  PACF  -  Raw  Data  -  Special 
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Figure  F-76.  ACF  -  Lagged  12  Months 
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Figure  F-77.  PACF  -  Lagged  12  Months 
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Figure  F-79.  Estimated  Residuals  r  Special 
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Figure  F-80.  Forecast  FY  82  -  Special 
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F-8.  CARGO  TRAILER/CONEX 

a.  The  raw  series  for  cargo  trailer  shipments  is  depicted  in  Figure 
F-81.  The  raw  data  do  exhibit  a  seasonal  pattern  that  trends  upwards  over 
time.  The  increasing  trend  over  time  is  more  clearly  evident  after  review¬ 
ing  the  3-month  and  6-month  moving  average  diagrams  (Figures  F-82  and 
F-83).  According  to  the  12-month  moving  average  diagram  (Figure  F-84),  the 
increasing  trend  appears  to  be  an  increasing  linear  function  over  time. 
However,  it  should  be  noted  that  over  the  past  2  years  (FY  81-FY  83),  the 
trend  has  remained  flat. 

b.  Autocorrelation  function  analysis  (Figures  F-85  and  F-86)  of  the 
data  confirms  the  existence  of  the  seasonality  in  the  data  as  indicated  by 
the  damping  peaks  at  12,  24,  and  36  lags.  To  eliminate  the  seasonality, 
the  series  was  differenced  at  a  lag  of  12  months.  Figure  F-87  depicts  the 
seasonally  adjusted  series.  The  autocorrelation  and  partial  autocorrela¬ 
tion  functions  of  the  seasonally  adjusted  series  are  illustrated  in  Figures 
F-88  and  F-89.  The  autocorrelation  function  indicates  that  the  series  is 
stationary  and  that  a  mean  exists  in  the  series  (t-value  of  mean  is  signif¬ 
icant  at  the  5  percent  significance  level). 

c.  Analysis  of  the  autocorrelation  and  partial  autocorrelation 
functions  suggests  that  the  data  could  be  modeled  with  an  ARIMA  (0,0,1) 
nonseasonal  model  and  an  ARIMA  (0,1,1)  seasonal  model.  Thus,  the  suggested 
model  takes  the  form: 

(1  -  B12)  Xt  -  m  +  (1  -  B)  (1  -  elz  Bl2)et 

d.  The  results  of  the  model  (Figure  F-90)  are  as  follows: 

(1  -  Bl2)  Xt  =  .4057  +  (1  +  .29128)  (1  -  .8264Bl2)*t 
(4.09)  (-2.32)  (14.97) 

x2  (2,20)  =  12.74 

The  results  indicate  that  all  of  the  estimated  parameters  are  significant 
at  the  5  percent  significance  level.  The  Q  statistic  for  the  estimated 
model  residuals  at  20  degrees  of  freedom  is  12.74  which  is  significant  at 
the  13  percent  significance  level  (Figure  F-91).  Therefore,  the  null 
hypothesis,  uncorrelated  residuals,  should  not  be  rejected  at  the  13  per¬ 
cent  significance  level. 

e.  The  model  forecasted  CONEX/cargo  trailer  shipments  for  FY  82  to  be 
76.43  K/MTON  (Figure  F-92)  versus  an  actual  lift  for  FY  82  of  74.47  K/MTON. 
The  MSE  of  the  forecast  was  .81  and  the  annual  forecast  error  was  2.6 
percent. 
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Figure  F-83.  Six-month  Moving  Average  -  CONEX 
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Figure  F-85.  ACF  -  Raw  Data  -  CONEX 
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Figure  F-89.  PACF  -  Lagged  12  Months 
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Figure  F-90.  Model  Parameters  -  CONEX 
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Figure  F-91.  Estimated  Residuals  -  CONEX 
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Figure  F-92.  Forecast  FY  82  -  CONEX 
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F-9.  FREEZE 

a.  The  raw  series  for  commodity  freeze  is  depicted  in  Figure  F-93.  The 
raw  series  does  not  depict  any  visible  pattern;  however,  the  data  exhibit 
an  upward  trend  over  time.  The  3-month  moving  average  (Figure  F-94)  indi¬ 
cates  that  the  data  may  be  seasonal  with  peaks  and  troughs  occurring  at  the 
end  and  beginning  of  the  fiscal  year,  respectively.  The  6-month  moving 
average  (Figure  F-95)  clearly  identifies  an  upward  trend  in  the  data,  but 
the  seasonal  trend  is  muted.  Finally,  the  12-month  moving  average  (Figure 
F-96)  indicates  a  steady  increase  of  tonnage  shipped  throughout  the  data 
period.  In  fact,  the  12-month  moving  average  suggests  that  the  annual  data 
follow  a  linear  trend  upwards. 

b.  Autocorrelation  function  analysis  (Figure  F-97)  of  the  data  affirms 
the  trend  and  seasonality  of  the  freeze  data.  The  trend  is  verified  by  the 
high  number  of  positive  coefficients  during  the  first  12  lags.  The  season¬ 
ality  of  the  data  is  suggested  by  the  significant  coefficients  at  lags  11, 
23,  and  35.  Seasonality  in  the  data  was  eliminated  by  differencing  the 
series  at  a  lag  of  12  months.  The  plot  of  the  seasonally  adjusted  series 
is  shown  in  Figure  F-98.  The  autocorrelation  and  partial  autocorrelation 
function  of  the  seasonally  adjusted  series  (Figures  F-99  and  F-100)  indi¬ 
cate  that  stationarity  in  the  data  has  not  been  achieved.  To  achieve 
stationarity  the  seasonally  adjusted  data  was  differenced  at  a  lag  of  one 
month.  The  resulting  autocorrelation  and  partial  autocorrelation  functions 
are  depicted  in  Figures  F-101  and  F-102. 

c.  The  autocorrelation  function  of  the  differenced  series  suggests  the 

possibility  of  a  mixed  nonseasonal  model  and  a  MA  seasonal  model.  The  in¬ 
itial  model  hypothesized  was  an  ARIMA  (1,1,1)  nonseasonal  model  and  an 
ARIMA  (0,1,1)  seasonal  model.  The  autocorrelation  function  of  the 
estimated  residuals  for  the  model  depicted  significant  coefficients  at 
lags  2,3,5,  and  11.  Next  a  ARIMA  (2,1,2)  nonseasonal  model  was  tried.  The 

coefficient  for  lag  2  was  statistically  insignificant  and  the  autocor¬ 
relation  coefficient  at  lag  11  of  the  estimated  residual  was  still  sig¬ 
nificant.  Finally,  an  ARIMA  (2,1,2)  model  incorporating  the  parameter 

11  as  the  second  MA  parameter  was  attempted.  The  form  of  this  model  was: 


(1-4>1  8  -  <t>2  B2)  (1-812)  (l-B)  Xt  =  (1-  h  B  -  BU)  U“  ^B12)^ 
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d.  The  results  of  the  model  (Figure  F-103)  are  as  follows: 


(1+.7740B+.6789B2)  (1-B12)  (1_b)  Xt  =  (1-.4729B+.5038BH)  (l-.7753Bl2)«t 
(-6.76)  (-5.96)  (4.32)  (-5.76)  (12.04) 

X2  (5,20)  =  14.66 


The  results  indicate  that  all  of  the  estimated  parameters  of  the  model  are 
significant  at  the  5  percent  significance  level  (Figure  F-104).  The  Q 
statistic  for  20  degrees  of  freedom  is  14.66.  Therefore,  the  null  hy¬ 
pothesis,  residuals  are  uncorrelated,  should  not  be  rejected  at  the  22 
percent  significance  level. 

e.  The  freeze  seasonal  model  forecasted  that  31.9  K/MTON  would  be 
transported  during  FY  82  versus  an  actual  shipment  of  36.4  K/MTON  (Figure 
F-105).  The  MSE  of  the  forecast  was  .49,  and  the  annual  forecast  error  was 
-12.5  percent. 
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Figure  F-96.  Twelve-month  Moving  Average  -  Freeze 
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Figure  F-101.  ACF  -  Lagged  12  Months  and  1  Month 
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Figure  F-102.  PACF  -  Lagged  12  Months  and  1  Month 
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Figure  F-104.  Estimated  Residuals 
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F-10.  CHILL 


a.  The  raw  series  for  chill  is  shown  in  Figure  F-106.  The  raw  series 
indicates  the  evidence  of  a  recurring  pattern  and  the  possibility  of 
seasonality  in  the  data.  Visual  examination  of  the  plots  of  the  3-month 
and  6-month  (Figures  F-107  and  F-108)  moving  averages  indicates  that  the 
data  do  follow  a  seasonal  trend.  However,  the  12-month  moving  average 
indicates  that  the  overall  mean  of  the  series  during  the  last  6  1/2  years 
has  been  constant  (Figure  F-109). 

b.  Autocorrelation  function  analysis  of  the  modeled  data  confirms  the 
fact  that  seasonality  does  exist  in  the  data  (Figures  F-110  and  F-lll). 
Therefore,  to  eliminate  seasonality,  the  series  was  differenced  at  a  lag  of 
12  months.  This  process  did  not  produce  stationary  data  since  the 
resulting  data  plot  exhibited  heteroscedasticity  (change  in  variance  over 
time).  The  original  series  was  then  transformed  using  natural  logarithms 
and  the  logarithmic  series  was  then  differenced  as  a  lag  of  12  months.  The 
plot  of  the  logarithmic  differenced  series  is  depicted  in  Figure  F-112. 

The  autocorrelation  and  partial  autocorrelation  functions  of  the  seasonally 
adjusted  series  are  shown  in  Figures  F-113  and  F-114.  The  autocorrelation 
function  suggests  that  the  data  is  stationary. 

c.  Analysis  of  the  autocorrelation  and  partial  autocorrelation  function 
suggests  the  possibility  of  an  ARIMA  (0,0,2)  nonseasonal  model  and  an  ARIMA 
(1,1,0)  seasonal  model.  An  AR  seasonal  process  was  used  since  the  results 
from  the  AR  model  were  better  than  the  MA  seasonal  model.  The  suggested 
model  is  of  the  form: 

(1  -  4>12B^)(1  ■  B^2)  in  =  (1  -  0jB  -  #2^)*  t 


d.  The  results  of  the  model  (Figure  F-115)  are  as  follows: 

(1  +  .637lBl2)(l  -  Bl2)  In  Xt  =  (1  -  .3012B  +  .2716B2)et 
(-5.61)  .  (2.0)  (-1.83) 

X2  (3.20)  =  15.63 


All  of  the  estimated  parameters  except  6 2  are  significant  at  the  5  percent 
significance  level  (Figure  F-116).  The  Q  statistic  for  20  degrees  of  free¬ 
dom  is  15.63,  which  indicates  that  the  residuals  are  not  significant  and 
that  the  null  hypothesis  should  not  be  rejected  at  the  26  percent  signifi¬ 
cance  level.  Also,  all  of  the  estimated  residuals  fall  within  the  95  per¬ 
cent  confidence  interval  and  apprear  to  be  distributed  randomly.  The  FY  82 
forecast  for  chill  using  this  model  was  15.32  K/MTON  of  cargo  shipped 
versus  an  actual  lift  of  14.13  K/MTON.  The  MSE  of  the  chill  forecast  was 
.21,  and  the  annual  forecast  error  was  8.4  percent. 
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igiire  F- 107 .  Three-month  Moving  Average  -  Chill 


Figure  F-109.  Twelve-month  Moving  Average 
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Figure  F-110.  ACF  -  Raw  Data  -  Chill 
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Figure  F-lll.  PACF  -  Raw  Data 


F-128 


autocorrelations 


CAA-SR-84-2 


PARTIAL 

AUTOCORRELATIONS 

1-  12 

-.15 

.20 

.29 

.07 

-.02 

.13 

-.15 

-.21 

.07 

-.24 

.06 

-.20 

ST.E. 

.14 

.14 

•  14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

13-  2H 

-.1C 

.10 

-.05 

.08 

.05  - 

.21 

.16 

.03 

.06 

-.06 

.08 

-.36 

ST.E. 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

25-  36 

-.21 

.17 

.07 

-.07 

-.10 

.02 

-.08 

.06 

-.04 

-.OS 

-.04 

-.03 

ST.E. 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

.14 

PLOT  OF 

PARTIAL 

AUTOCORRELATIONS 

LAfi  CC 

-1.0 

>RR. 

-.8 

-.6 

-.4 

-.2 

.0 

.2 

•—4  — — 

.4 

.6 

.8 

1.0 

—4 

1 

-.147 

♦ 

xxxx  i 

4 

2 

.199 

♦ 

IXXXXX 

4 

3 

.291 

4 

1XXXXXXX 

4 

.072 

4 

IXX 

4 

5 

-.017 

4 

1 

4 

6 

.126 

4 

IXXX 

4 

7 

-.146 

4 

XXXX  I 

4 

a 

-.214 

4 

XXXXX  I 

4 

9 

.U75 

4 

IXX 

4 

10 

-.238 

«  XXXXXX  1 

4 

11 

.059 

4 

IX 

4 

12 

-,2C2 

4 

XXXXX I 

4 

13 

-.097 

4 

XX  1 

4 

14 

.101 

4 

IXXX 

4 

15 

-.052 

4 

XI 

4 

16 

.082 

4 

IXX 

4 

17 

.051 

4 

IX 

4 

18 

-.210 

4 

XXXXX  I 

4 

19 

.165 

4 

IXXXX 

4 

20 

.034 

4 

IX 

4 

21 

.062 

4 

IXX 

4 

4? 

-.063 

4 

XX  I 

4 

23 

.077 

4 

IXX 

4 

24 

-.057 

4 

X  I 

4 

25 

-.214 

4 

XXXXX  I 

4 

I6 

.168 

4 

IXXXX 

4 

27 

.067 

4 

IXX 

4 

28 

-.069 

4 

XX  I 

4 

29 

-.099 

4 

XX  I 

4 

30 

.013 

4 

I 

4 

31 

-.084 

4 

XX  T 

4 

32 

.061 

4 

IXX 

4 

33 

-.036 

4 

X  I 

4 

34 

-.049 

4 

X  I 

4 

35 

-.043 

4 

XI 

4 

36 

-.034 

4 

X  I 

4 

Figure  F-114.  PACF  Lagged  12  Months 
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ESTIMATION  BY  B ACKCAST1N6  METHOD 

RELATIVE  CHANGE  IN  RESIDUAL  SUN  OF  SQUARES  LESS  THAN  .1000-03* 
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Figure  F-115.  Model  Parameters  -  Chill 


F-132 


CAA-SR-84-2 


AUTOCORRELATIONS 


It.!! 

Jl3 

*:?$ 

:!|-:K  :» 

:il 

0.0  .07  -.21 

•14  .14  .14 

:B  -:B 

13-29 

ST.E. 

":iS 

:?!’ 

*:li  ":Si  -:B 

•  04 

•  16 

.11  -.OS  .09 
.16  .16  .16 

•09  -.17 
.16  .16 

25-  36 
ST.E. 

•  CS 

•  17 

s» 

:» 

W  *:??  -:B 

—.36  0.3 

.17  .17 

PLOT  OF  AUTOCORRELATIONS 

•!  «0  "*  •  3  -#6  -«2  »Q  i2  *4  *6  »8  1*0 


L«b 

CORK .  1 

■4--—  —  4--— 

I 

1 

•  038 

IX 

2 

-.039 

XI  < 

3 

•2ca 

IXXXXX< 

4 

— .066 

XXI 

4 

S 

.064 

4 

IXX 

4 

6 

.107 

4 

IXXX 

4 

7 

-.120 

4 

XXXI 

4 

a 

-.004 

4 

I 

9 

10 

.071 

-.208 

4 

4 

IXX 

XXXXXI 

4 

4 

11 

.083 

4 

IXX 

4 

12 

-.035 

, 

4 

XI 

4 

13 

-.198 

4 

XXXXXI 

4 

14 

.042 

4 

IX 

4 

15 

-.124 

4 

XXXI 

4 

16 

-.046 

4 

XI 

17 

-.088 

4 

XXI 

4 

18 

-.166 

4 

XXXXI 

4 

19 

.044 

4 

IX 

4 

20 

.115 

4 

IXXX 

4 

21 

-.084 

4 

XXI 

4 

Si 

4 

4 

IXX 

IXX 

4 

4 

24 

-.174 

4 

XXXXI 

4 

25 

.047 

• 

.  -4  . 

—  -  IX 

4 

26 

.044 

4 

IX 

4 

27 

•  U41 

4 

IX 

4 

28 

.027 

4 

IX 

4 

29 

-.088 

4 

XXI 

4 

30 

•  085 

4 

IXX 

4 

31 

.005 

4 

•  l 

4 

32 

-.005 

4 

I 

4 

33 

-.069 

*  ■ 

XXI 

34 

-.064 

4 

XXI 

4 

35 

-.058 

4 

XI 

4 

36 

-.005 

4 

I 

4 

Figure  F-116.  Estimated  Residuals  -  Chill 
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APPENDIX  6 

WINTERS  MODEL  INITIALIZATION  AND  SMOOTHING  CONSTANTS 


G-I.  GENERAL.  This  section  addresses  initialization  of  model  parameters 
and  development  of  smoothing  constants. 

6-2.  INITIALIZATION.  Chapter  6  discussed  the  procedure  for  updating  the 
parameters  (permanent  component,  trend,  and  seasonal  factor),  given  that 
initial  values  exist.  Upon  option,  initial  estimates  of  the  Winters  model 
parameters  can  be  specified  by  the  user.  Alternatively,  several  heuristic 
algorithms  have  been  devised  to  initialize  parameters  based  on  manipulation 
of  historical  data.  The  initialization  procedure  described  below  is  due  to 
Montgomery, 5  and  is  similar  to  the  one  proposed  by  Winters. 

a.  Trend  Component.  Assuming  that  data  are  available  for  m  seasons, 
then  compute  the  mean  of  all  observations  for  the  first  and  last  of  these 
seasons.  Denote  the  average  observation  for  the  jth  season  by  Xj,  j  = 
l,2,...,m.  Estimate  the  trend  in  the  same  manner  that  would  be  used  to 
compute  a  simple  algebraic  slope.  Since  there  are  m-1  seasons  between 
season  1  and  season  m,  and  since  there  are  L  periods  per  season,  then  the 
initial  estimate  of  the  trend  becomes 


a  xm  -  x, 

b  (0)  *  JO - I 

2  (m  -  1)  L 


(G-l) 


b.  Permanent  Component.  For  initialization  purposes,  it  is  assumed 
that  the  average  observation  iq  for  the  first  season  occurs  timewise  at  the 
middle  of  the  season.  With  this  in  mind,  the  permanent'  component  can  be 
treated  like  a  simple  y-intercept.  Writing  the  equation  in  si ope- intercept 
form  gives 


*1  =  %(0) 


4- V0) 


Since  all  terms  are  known  except  for  the  permanent  component,  equation 
(G-2)  can  be  rewritten  as 


%<0) 


-  -j-fyo) 


G-l 
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c.  Seasonal  Factor.  Since  there  are  m  seasons  and  L  periods  per 
season,  seasonal  factors  are  computed  initially  for  each  of  the  mL  periods. 
Each  factor  is  computed  as  the  ratio  of  the  actual  observation  to  the 
average  seasonally  adjusted  value  for  that  season,  further  adjusted  by  the 
trend.  The  computation  is 


x,  -[(L  +  l)/2  -  jy&^O) 


t  *  1,2,..., mL 


(6-4) 


where  Xj  is  the  average  for  a  season  corresponding  to  the  t  index,  and  j 
is  the  position  of  period  t  within  the  season.  For  example,  if  l<t<L, 
then  i  =  1,  and  if  L  +  l<t<2L,  then  i  =  2.  Equation  (G-4)  produces  m 
estimates  of  the  seasonal  factor  for  each  period.  (In  the  TWF  Study,  m  was 
usually  five,  and  there  were  five  estimates  for  each  month  of  the  year.) 

The  m  estimates  for  each  period  (month)  are  averaged  to  produce  a  single 
estimate  of  the  seasonal  factor  for  each  period  within  the  season. 


c.  =  - 


m  -  1 

Z 

k  =  0 


t  +  kL 


t  *  1,2,. ..,L 


(G-5) 


Finally,  the  seasonal  factors  are  normalized  so  that  they  sum  to  L  (L  =  12 
in  the  study). 


t  =  1,2, . . . ,L  (G-6) 


The  above  procedure  produces  estimates  ^i(O),  ^(O),  and  ^t(O)  assuming 
that  the  origin  of  time  is  immediately  prior  to  period  1.  The  parameters 
may  then  be  updated  by  the  technique  described  in  paragraph  6-5  of  this 
report. 


G-2 
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S-3.  SMOOTHING  CONSTANTS.  Smoothing  constants  are  necessary  in  order  to 
combine  (weight)  previous  estimates  of  parameters  with  their  updated 
values.  Numerical  estimates  of  the  permanent  component,  trend  component, 
and  seasonal  factor  receive  the  weights  a  ,  0  ,  and  7  for  the  current 
interval  T.  These  weighted  estimates  are  combined  additively  with 
complementary  weighted  values  (using  1-a,  1-/3  ,  and  1-7  )  for  the 
previous  time  period  or  season,  as  appropriate.  All  weights  are  varied 
incrementally  so  that  the  parameters  of  the  model  ultimately  provide  the 
best  fit  according  to  some  predetermined  criterion,  i.e.,  mean  square 
error.  Unlike  the  formal  method  of  least  squares,  which  uses  partial 
derivatives  to  develop  a  set  of  simultaneous  linear  equations  (normal 
equations)  that  are  solved  through  matrix  inversion,  the  Winters  method  is 
heuristic  in  nature.  As  such,  the  optimum  set  of  smoothing  constants  is 
determined  by  trial  and  error.  The  coefficients  lie  in  the  interval  (0,1). 
In  order  to  keep  computer  time  requirements  modest,  a  coarse  grid  is  tried 
first.  Values  of  a  ,  0  ,  and  7  are  stepped  across  the  unit  interval  in 
increments  of  0.05  until  all  possible  combinations  of  smoothing  constants 
have  been  examined.  The  set  of  (  <*,/3  ,7  )  producing  the  smallest  mean 
square  error  is  used  in  the  program  as  the  basis  for  a  second,  fine-grained 
search.  A  step  of  0.01  is  used  to  search  in  a  narrow  interval  about  the 
coarse  estimates  of  (a,/3  ,7)  to  yield  refined  values  of  the  smoothing 
constants. 
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THE  WINTERS  METHOD  COMMODITY  ANALYSIS 


For  the  nine  commodities  evaluated  in  the  study,  this  appendix  contains 
three  exhibits  each: 

a.  The  first  exhibit  for  each  commodity  displays  the  top  40  sets  of 
smoothing  constants  obtained  during  the  Winters  model  optimization  process. 
The  sets  of  (a  ,  0  ,  y  )  are  rank-ordered  according  to  the  residual  sums  of 
squares  they  produce.  Since  each  smoothing  constant  is  stepped  across  the 
unit  interval  (0,1)  in  increments  of  0.01,  there  are  1003,  0r  one  million, 
candidate  sets  of  smoothing  constants  evaluated  for  each  commodity.  The 
model  ultimately  defined  for  each  commodity  is  the  one  yielding  the 
smallest  residual  sum  of  squares  over  the  data  interval  defined  up  to,  but 
not  including,  the  forecast  interval.  This  set  of  optimum  smoothing 
constants  is  then  presumed  to  yield  the  best  fit  over  the  forecast 
interval.  It  should  be  noted  that  although  this  selection  criterion 
provides  the  best  fit  over  the  initialization  interval,  it  does  not  always 
guarantee  the  best  fit  over  the  forecast  interval. 

b.  Tne  second  exhibit  for  each  commodity  displays  both  the  raw  and 
fitted  time  series  over  the  entire  model  initialization  interval,  i.e.,  up 
to  but  not  including  the  forecast  interval.  For  most  of  the  nine  commodi¬ 
ties,  60  months  (5  years)  worth  of  data  were  used  in  the  initialization 
phase  to  formulate  the  optimum  model.  Because  transportation  workload 
forecasting  requires  a  lead  time  of  at  least  12  months,  the  parameters 
(permanent  component,  trend,  and  seasonal  factor)  from  the  last  12  months 
(usually  months  49  to  60)  of  the  initialization  phase  were  used  to  produce 
the  12  monthly  estimates  of  the  forecast  phase. 

c.  The  third  exhibit  for  each  commodity  shows  the  actual  forecast 
phase.  Summary  statistics  are  provided  at  the  end  of  each  printout.  Once 
the  Winters  model  smoothing  constants  and  parameters  have  been  developed 
internally,  calculation  of  forecasts  is  straightforward-  For  example, 
applying  the  forecast  formula 


wn  -  [*.(7-)  +  *i(7>]cr+,(r  +  T  -  L) 


to  obtain  an  estimate  for  POVs  shipped  during  month  67  would  require  use  of 
parameters  from  month  55  (T=55;  r=12;  L=12). 

x67  (55)  *  [40.1948  +  (0.0376) (12)]  0.8949 

=  36.37 
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Figure  H-2.  Output  of  the  Initialization  Phase 
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Figure  H-8.  Output  of  the  Initialization  Phase  -  HHG 
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Figure  H-10.  Smoothing  Constant  Optimization  Routine  -  Coal 
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Figure  H-ll.  Output  of  the  Initialization  Phase  -  Coal 
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Figure  H-13.  Smoothing  Constant  Optimization  Routine  -  Ammunition 
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Figure  H-16.  Smoothing  Constant  Optimization  Routine  -  Special 
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Figure  H-17.  Output  of  the  Initialization  Phase  -  Special 
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Figure  H-20.  Output  of  the  Initialization  Phase  -  CONEX 
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Figure  H -22.  Smoothing  Constant  Optimization  Routine  -  Freeze 
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Figure  H-26.  Output  of  the  Initialization  Phase  -  Chill 
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DALO-TSP-C 1 1  (30  Jan  84)  1st  Ind 

SUBJECT:  Transportation  Workload  Forecasting  (TWF)  Study 

DA  Washington,  D.C.  20310  (DALO-TSP-C 1 1 )  2  MAR  1984 

TO:  Director,  US  Army  Concepts  Analysis  Agency,  ATTN:  CSCA-SPP,  8120 
Woodmont  Avenue,  Bechesda,  MD  20814 

Completed  Study  Critique  is  furnished.  There  was  one  editorial  comment,  which 
is  lisced  on  a  separate  page  and  attached  to  the  critique  sheet. 

FOR  THE  DEPUTY  CHIEF  OF  STAFF  FOR  LOGISTICS: 

h^jQQOaJk. 

1  Enel  •BERNARD  J.  CLARK 

as  Colonel,  GS 

Chief,  Performance  Management  Division 
Directorate  of  Transportation,  Energy 
and  Troop  Support 
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DEPARTMENT  OF  THE  ARMY 

US  ARMY  CONCEPTS  ANALYSIS  AGENCY 
8120  WOOOMONT  AVENUE 
SETHESOA.  MARYLAND  20814 


CSCA-SPP 

SUBJECT:  Transportation  Workload  Forecasting  (TWF)  Study 


3  0  JAN  1934 


Deputy  Chief  of  Staff 
for  Logistics 
Department  of  the  Army 
ATTN:  DALO-TSP 
Washington,  DC  20310 


1.  Reference: 

a.  Letter,  DAL0-TSP-C11,  11  May  1983,  subject:  Transoortation  Workload 
Forecasting  (TWF)  Study. 

b.  Letter,  DACS-DMO,  19  October  1983,  subject:  Responsibility  of  Study 
Performing  and  Study  Sponsoring  Organizations. 

2.  The  Deputy  Chief  of  Staff  for  Logistics  requested  that  the  US  Army 
Concepts  Analysis  Agency  study  US  Army  transportation  workload  forecasting 
and  develop  procedures  to  improve  the  system. 

3.  Attached  are  several  copies  of  our  Draft  Study  Report  which  describe 
the  study  approach,  the  current  forecasting  system,  and  several  alternative 
systems  and  methods  that  should  result  in  improved  forecasts  of  over-ocean 
cargo  transportation  requirements.  These  drafts  are  being  provided  IAW 
reference  b  in  order  to  obtain  your  comments  prior  to  publication  of  the 
final  report  of  study.  A  suggested  Study  Critique  Sheet  is  provided  for 
your  use  as  you  desire. 

4.  Request  that  your  comments  be  provided  to  CAA  within  30  days  after 
receipt  of  the  draft  report.  Your  comments,  if  any,  and  our  response 
to  comments,  if  any,  will  be  included  in  the  final  report  if  they  are 
provided  to  CAA  prior  to  our  planned  publication  date. 


2  Incl 
as 


DAVID  C.  HARDISON 
Director 


1-2 


1.  Were  there  any  editorial  comments?  YFS  If  so,  please  list  on 
separate  page  and  attach  to  the  critique  sheet. 

2.  Was  the  work  accomplished  in  a  timely  manner?  yes _ .  If  not, 

please  comment.  _ 


3.  Does  the  work  report  address  adequately  the  issues  planned  for  the 
analysis?  yes _ .  If  not,  please  comment.  _ 


4.  Were  appropriate  analysis  techniques  used?  yes _ .  If  not, 

please  comment.  _ _ 


5.  Are  the  findings  fully  supported  by  good  analysis  based  on  sound 
assumptions?  yes _ .  If  not,  please  explain. _ 


6.  Does  the  report  contain  the  preferred  level  of  details  of  the 
analysis?  ves _ .  If  not,  please  comment.  _ 


7.  Is  the  written  material  fully  satisfactory  in  terms  of  clarity  of 

presentation,  completeness,  and  style?  yes _ •  If  not,  please 

comment.  _  _  _ • 
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PARAGRAPH  6-8a  should  read  as  follows: 

a*  Transportation  arrangements  for  DOD-sponsored  household  goods  (HHG )  are 
made  by  MTMC,  using  primarily  two  methods:  freight  forward  (Code  4),  and  those 
that  arc  MSC  processed  (Code  5).  Freight  forward  shipments  are  shipments  of  HHG 
that  are  Contracted  with  private  carriers  from  door-to-door;  shipments  of  HHG 
that  are  shipped  using  the  assets  of  MSC  are  Code  5  shipments.  Additionally, 
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2.  MODELS,  ROUTINES,  ANO  SIMULATIONS 
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THE  PRINCIPAL  FINDINGS  of  the  work  reported  herein  are  as  follows: 

i 

(1)  The  transportation  workload  forecasting  system  has  produced 
inaccurate  forecasts  resulting  in  inefficient  Military  Sealift  Command 
(MSC)  industrial  fund  operations. 

(2)  Accurate  forecasting  of  cargo  transportation  requirements  can  be 
accomplished  by  forecasting  at  a  single  activity. 

(3)  Either  HQ  Military  Traffic  Management  Command  (MTMC)  or  HQ,  US 
Army  Materiel  Development  and  Readiness  Command  (DARCOM)  is  a  suitable 
location  for  a  single  point  forecasting  activity. 

(4)  The  Box-Jenkins  and  Winters  Forecasting  Models  can  provide 
accurate  forecasts  when  used  in  conjunction  with  program  information. 

(5)  Changes  to  the  allocation  of  transportation  account  codes  and 
requirements  for  forecasting  shipping  mode  are  also  required  to  improve 
forecasting  accuracy. 


THE  MAIN  ASSUMPTION  on  which  the  work  reported  herein  rests  is  that  trans¬ 
portation  workload  forecasting  requirements,  contained  in  JCS  Publication 
15,  would  not  be  changed. 


THE  PRINCIPAL  LIMITATIONS  of  this  work  which  may  affect  the  findings  are  as 
follows: 

(1)  Only  the  forecasting  of  peacetime  over-ocean  surface  cargo  trans¬ 
portation  requirements  was  evaluated. 

(2)  Historical  lift  data  was  extracted  exclusively  from  MSC  records 
and  could  not  be  validated  from  Army  sources. 


THE  SCOPE  OF  THE  STUDY  was  taken  to  include  an  analysis  of  the  Army's  long- 
range  cargo  transportation  requirements  forecasting  process  and  its  impact 
on  budgets  and  transportation  costs. 


THE  STUDY  OBJECTIVE  was  to  develop  cost  effective  systems  and  methods  for 
improving  the  forecasting  of  Army  over-ocean  surface  cargo  transportation 
requirements. 


THE  BASIC  APPROACH  followed  in  this  study  can  be  defined  as:  research  was 
conducted  into  the  nature  and  extent  of  the  forecasting  problem,  to  identify 
its  impact,  and  its  systemic  and  methodological  causes.  Several  alternative 
systems  were  evaluated  based  on  their  relative  costs  and  efficiency.  Then 
a  series  of  mathematical  techniques  was  evaluated  for  suitability  as  fore¬ 
casting  tools.  Two  of  the  techniques,  the  Box-Jenkins  and  Winters  models, 
were  used  to  forecast  the  1982  cargo  transportation  requirements  based  on 
1977  to  1981  MSC  cargo  lift  data. 


REASONS  FOR  PERFORM I NS  THE  STUDY  are  mainly  as  follows:  recent  forecasts 
of  Army  over-ocean  surface  cargo  transportation  requirements  have  been  in¬ 
accurate.  As  a  consequence  MSC  industrial  funds  have  incurred  significant 
losses  and  the  MSC  controlled  fleet  was  not  efficiently  utilized  for  cargo 
transport.  This  study  was  directed  to  develop  methods  to  improve  the  fore¬ 
casts. 


THE  STUDY  SPONSOR  was  the  Deputy  Chief  of  Staff  for  Logistics,  who  also 
established  the  objectives  and  monitored  the  study  activities. 


THE  STUOY  EFFORT  was  directed  by  LTC  James  N.  Keenan,  Strategy,  Concepts 
and  Plans  Directorate. 


COfflENTS  AND  QUESTIONS  may  be  directed  to  CAA,  ATTN:  Assistant  Director 
for  Strategy,  Concepts  and  Plans. 
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requirements  for  forecasting  shipping  mode  are  also  required  to  improve 
forecasting  accuracy. 

THE  MAIN  ASSUMPTION  on  which  the  work  reported  herein  rests  is  that  trans- 
portation  workload  forecasting  requirements,  contained  in  JCS  Publication 
15,  would  not  be  changed. 

TOE  PRINCIPAL  LIMITATIONS  of  this  work  which  may  affect  the  findings  are  as 
follows: 

(1)  Only  the  forecasting  of  peacetime  over-ocean  surface  cargo  trans¬ 
portation  requirements  was  evaluated. 

(2)  Historical  lift  data  was  extracted  exclusively  from  MSC  records 
and  could  not  be  validated  from  Army  sources. 


TOE  SCOPE  OF  THE  STUDY  was  taken  to  include  an  analysis  of  the  Army's  long- 
range  cargo  transportation  requirements  forecasting  process  and  its  impact 
on  budgets  and  transportation  costs. 


THE  STUDY  OBJECTIVE  was  to  develop  cost  effective  systems  and  methods  for 
improving  the  forecasting  of  Army  over-ocean  surface  cargo  transportation 
requirements. 


THE  BASIC  APPROACH  followed  in  this  study  can  be  defined  as:  research  was 
conducted  into  the  nature  and  extent  of  the  forecasting  problem,  to  identify 
its  impact,  and  its  systemic  and  methodological  causes.  Several  alternative 
systems  were  evaluated  based  on  their  relative  costs  and  efficiency.  Then 
a  series  of  mathematical  techniques  was  evaluated  for  suitability  as  fore¬ 
casting  tools.  Two  of  the  techniques,  the  Box-Jenkins  and  Winters  models, 
were  used  to  forecast  the  1982  cargo  transportation  requirements  based  on 
1977  to  1981  MSC  cargo  lift  data. 


REASONS  FOR  PERFORMING  THE  STUDY  are  mainly  as  follows:  recent  forecasts 
of  Army  over-ocean  surface  cargo  transportation  requirements  have  been  in¬ 
accurate.  As  a  consequence  MSC  industrial  funds  have  incurred  significant 
losses  and  the  MSC  controlled  fleet  was  not  efficiently  utilized  for  cargo 
transport.  This  study  was  directed  to  develop  methods  to  improve  the  fore¬ 
casts. 


THE  STUDY  SPONSOR  was  the  Oeputy  Chief  of  Staff  for  Logistics,  who  also 
established  the  objectives  and  monitored  the  study  activities. 


THE  STUDY  EFFORT  was  directed  by  LTC  James  N.  Keenan,  Strategy,  Concepts 
and  Plans  Directorate. 


COWOTS  AND  QUESTIONS  may  be  directed  to  CAA,  ATTN:  Assistant  Director 
for  Strategy,  Concepts  and  Plans. 
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THE  PRINCIPAL  FINDINGS  of  the  work  reported  herein  are  as  follows: 

(1)  The  transportation  workload  forecasting  system  has  produced 
inaccurate  forecasts  resulting  in  inefficient  Military  Sealift  Command 
(MSC)  industrial  fund  operations. 

(2)  Accurate  forecasting  of  cargo  transportation  requirements  can  be 
accomplished  by  forecasting  at  a  single  activity. 

(3)  Either  HQ  Military  Traffic  Management  Command  (MTMC)  or  HQ,  US 
Army  Materiel  Development  and  Readiness  Command  (DARCOM)  is  a  suitable 
location  for  a  single  point  forecasting  activity. 

(4)  The  Box-Jenkins  and  Winters  Forecasting  Models  can  provide 
accurate  forecasts  when  used  in  conjunction  with  program  information. 

(5)  Changes  to  the  allocation  of  transportation  account  codes  and 
requirements  for  forecasting  shipping  mode  are  also  required  to  improve 
forecasting  accuracy. 


THE  MAIN  ASSUMPTION  on  which  the  work  reported  herein  rests  is  that  trans¬ 
portation  workload  forecasting  requirements,  contained  in  JCS  Publication 
15,  would  not  be  changed. 


THE  PRINCIPAL  LIMITATIONS  of  this  work  which  may  affect  the  findings  are  as 
follows: 

(1)  Only  the  forecasting  of  peacetime  over-ocean  surface  cargo  trans¬ 
portation  requirements  was  evaluated. 

(2)  Historical  lift  data  was  extracted  exclusively  from  MSC  records 
and  could  not  be  validated  from  Army  sources. 


THE  SCOPE  OF  THE  STUDY  was  taken  to  include  an  analysis  of  the  Army's  long- 
range  cargo  transportation  requirements  forecasting  process  and  its  impact 
on  budgets  and  transportation  costs. 


THE  STUDY  OBJECTIVE  was  to  develop  cost  effective  systems  and  methods  for 
improving  the  forecasting  of  Army  over-ocean  surface  cargo  transportation 
requirements. 


THE  BASIC  APPROACH  followed  in  this  study  can  be  defined  as:  research  was 
conducted  into  the  nature  and  extent  of  the  forecasting  problem,  to  identify 
its  impact,  and  its  systemic  and  methodological  causes.  Several  alternative 
systems  were  evaluated  based  on  their  relative  costs  and  efficiency.  Then 
a  series  of  mathematical  techniques  was  evaluated  for  suitability  as  fore¬ 
casting  tools.  Two  of  the  techniques,  the  Box-Jenkins  and  Winters  models, 
were  used  to  forecast  the  1982  cargo  transportation  requirements  based  on 
1977  to  1981  MSC  cargo  lift  data. 


REASONS  FOR  PERFORMING  THE  STUDY  are  mainly  as  follows:  recent  forecasts 
of  Army  over-ocean  surface  cargo  transportation  requirements  have  been  in¬ 
accurate.  As  a  consequence  MSC  industrial  funds  have  incurred  significant 
losses  and  the  MSC  controlled  fleet  was  not  efficiently  utilized  for  cargo 
transport.  This  study  was  directed  to  develop  methods  to  improve  the  fore¬ 
casts. 


THE  STUDY  SPONSOR  was  the  Deputy  Chief  of  Staff  for  Logistics,  who  also 
established  the  objectives  and  monitored  the  study  activities. 


THE  STUDY  EFFORT  was  directed  by  LTC  James  N.  Keenan,  Strategy,  Concepts 
and  Plans  Directorate. 


COMMENTS  AND  QUESTIONS  may  be  directed  to  CAA,  ATTN:  Assistant  Director 
for  Strategy,  Concepts  and  Plans. 


